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20.  ABSTRACT  - Continued 

^ehavlor  of  the  compensating  circuits  have  been  derived  and  Indicate  that  the 
error  Is  proportional  to  the  square  of  the  temperature  variation.  For  typical 
probes  with  an  ambient  temperature  Increase  of  40oC,  maximum  errors  In  the 
Indicated  velocity  of  -IX  or  less  can  be  achieved. 

Measurements  In  many  complex  flowflelds,  such  as  In  swirling  flows  or  near 
luff  bodies,  require  the  knowledge  of  the  yaw  sensitivity  of  an  X-wlre  probe 
to  a high  degree  of  accuracy.  The  yaw  relation  Ue^^l  " U ~ b(l-cos°‘ a)  ] 2 
roposed  and  investigated  analytically  and  experimentally.  In  the  derivation 
of  the  response  equation  of  Inclined  cylindrical  sensors  In  steady  and  fluctuat- 
ing velocity  fields,  previous  analyses  have  been  extended  for  the  present  yaw 
relation  to  Include  X-probes  with  dissimilar  sensors  as  well  as  mean  velocities 
at  an  angle  y to  the  line  of  symmetry  of  the  X-wlres. 

For  mean  velocity  measurements,  with  sensors  that  have  yaw  dependence  which 
does  not  closely  follow  the  cosine  law,  neither  the  simple  nor  an  approximate 
weighted  sum  and  difference  technique  applies  when  y^O.  This  makes  analog 
processing  to  obtain  accurate  mean  velocities  extremely  difficult  If  not  Impos- 
sible. However,  for  cases  where  the  mean  streamline  is  aligned  with  the  probe 
axis  of  symmetry,  the  errors  Introduced  were  found  to  be  less  then  SX.  In 
low  level  turbulence,  time  dependent  velocity  components  can  always  be  obtained 
with  comprable  accuracy  using  a weighted  sum  and  difference  technique  once  the 
mean  flow  angle  Is  determined.  However,  the  weighting  must  take  into  account 
deviations  from  the  cosine  law.  For  a typical  sensor  whose  yaw  data  has  a 
standard  deviation  from  the  cosine  law  of  0.5X,  and  when  Y“0,  the  estimates  of 
errors  In  u',  v*  and  uv  are  approximately  3X,  IIX  and  8X,  respectively. 

The  error  In  uv  Increases  substantially  as  y deviates  from  zero.  These 
errors  are  essentially  proportional  to  the  standard  deviation  from  the  cosine 
law,  o,  as  It  Increases  above  its  value  In  this  example.  Even  when  the  sensor 
closely  follows  the  cosine  law,  l.e.,  o ; O.IX,  errors  in  these  quantities  as 
large  as  3X  can  exist.  Discrepancies  were  also  found  between  different  yaw 
relations  In  calculated  values  of  turbulence  velocity  correlations  which  are 
derived  from  the  same  data.  These  discrepancies  depend  on  the  mean  flow  angle, 
the  distribution  of  the  turbulence  energy  between  the  velocity  components  and 
the  accuracy  to  which  the  mean  flow  calibration  data  fits  the  yaw  function. 

The  effect  of  various  parameters  on  the  yaw  sensitivity  of  hot-wires  has  been 
experimentally  determined.  In  all  cases,  the  present  yaw  relation  fits  the 
calibration  data  to  a higher  degree  of  accuracy  than  the  cosine  law  or  the 
Frlehe  and  Schwarz  relation,  where  m ■ 0.5.  The  ratio  of  Improvement  over 
either  relation  Is  smaller  for  the  very  short  sensors  (1/d  r 100)  than  for  the 
long  sensors  (l/d  = 1000).  The  trend  of  the  yaw  parameters  with  either  sensor 
t/d  or  Reynolds  number  docs  not  seem  to  approach  the  cosine  law  but  rather  it 
seems  to  overshoot  It  without  any  observed  asymptotic  values.  Values  of  the 
transverse  sensitivity  coefficient  calculated  from  a yaw  sensitivity  function, 
taking  Into  account  deviations  from  the  cosine  law,  yield  a higher  degree  of 
accuracy  chan  those  obtained  directly  from  the  data.  Based  on  the  Reynolds 
number  dependency  of  the  yaw  coefficients,  one  should  use  different  yaw  cali- 
bration coefficients  for  different  velocities.  If  a single  set  of  parameters 
Is  to  be  used,  the  accuracy  achieved  by  a simple  relation  like  Che  cosine  law 
Is  as  good  as  one  may  expect.  However,  efter  Che  data  has  been  gathered,  Che 
results  must  be  corrected  for  Che  large  errors  Introduced  Into  the  mean  and 
rms  velocity  components.  Each  new  sensor  mounted  on  a probe  must  be  calibrated, 
even  If  the  new  sensor  Is  Identical  to  the  one  being  replaced. 
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ABSTRACT 


I 

A scheme  utilizing  a temperature  probe  immersed  in  the 
working  fluid  to  compensate  for  the  dependence  of  hot-wire 
velocity  calibration  on  ambient  temperature  variations  was 
analyzed.  One  needs  only  to  know  the  properties  of  the 
anemometer  bridge  and  the  velocity  and  temperature  probes 
to  achieve  the  compensation.  Hence,  the  scheme  provides 
means  for  incorporating  the  temperature  compensation 
a priori  to  conducting  the  experiments  without  any  need 
for  temperature  calibration.  Estimates  of  the  errors  intro- 
duced through  the  non-ideal  behavior  of  the  compensating 
circuits  have  been  derived  and  indicate  that  the  error  is 
proportional  to  the  square  of  the  temperature  variation. 

For  typical  probes  with  an  ambient  temperature  increase  of 
40®C, maximum  errors  in  the  indicated  velocity  of  -1%  or 
less  can  be  achieved. 

Measurements  in  many  complex  flowfields,  such  as  in 

swirling  flows  or  near  bluff  bodies , require  the  knowledge  of 

the  yaw  sensitivity  of  an  X-wire  probe  to  a high  degree  of 

accuracy.  The  yaw  relation  U /u.  = [1  - b (l-cos"’a)  ] ^ is 

e 1 

proposed  and  investigated  analytically  and  experimentally. 

In  the  derivation  of  the  response  equations  of  inclined 
cylindrical  sensors  in  steady  and  fluctuating  velocity 
fields,  previous  analyses  have  been  extended  for  the  pres- 
ent yaw  relation  to  include  X-probes  with  dissimilar  sen- 
sors as  well  as  mean  velocities  at  an  angle  y to  the  line  of 
symmetry  of  the  X-wires. 
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For  mean  velocity  measurements,  with  sensors  that  have 
yaw  dependence  which  does  not  closely  follow  the  cosine  law, 
neither  the  simple  nor  an  approximate  weighted  sum  and  dif- 
ference technique  applies  when  This  makes  analog 

processing  to  obtain  accurate  mean  velocities  extremely  dif- 
ficult if  not  impossible.  However,  for  cases  where  the  mean 
streamline  is  aligned  with  the  probe  axis  of  symmetry,  the 
errors  introduced  were  found  to  be  less  then  5%.  In  low 
level  turbulence,  time  dependent  velocity  components  can 
always  be  obtained  with  comprable  accuracy  using  a weighted 
sxim  and  difference  technique  once  the  mean  flow  angle  is 
determined.  However,  the  weighting  must  take  into  account 
deviations  from  the  cosine  law.  For  a typical  sensor  whose 
yaw  data  has  a standard  deviation  from  the  cosine  law  of 
0.5%,  and  when  y=0 , the  estimates  of  errors  in  u' , v'  and 
uv  are  approximately  3%,  11%  and  8%,  respectively.  The  error 
in  uv  increases  substantially  as  y deviates  from  zero.  These 
errors  are  essentially  proportional  to  the  standard  devia- 
tion from  the  cosine  law,  o,  as  it  increases  above  its  value 
in  this  example.  Even  when  the  sensor  closely  follows  the 
cosine  law,  i.e.,  o sO.1%,  errors  in  these  quantities  can 
exist.  Discrepancies  were  also  found  between  different  yaw 
relations  in  calculated  values  of  turbulence  velocity  cor- 
relations which  are  derived  from  the  same  data.  These  dis- 
crepancies depend  on  the  mean  flow  angle,  the  distribution 
of  the  turbulence  energy  between  the  velocity  components 
and  the  accuracy  to  which  the  mean  flow  calibration  data  fits 
the  particular  yaw  function. 
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The  effect  of  various  parameters  on  the  yaw  sensitivity  i 

of  hot-wires  has  been  experimentally  determined.  In  all 
cases,  the  present  yaw  relation  fits  the  calibration  data 

I 

to  a higher  degree  of  accuracy  than  the  cosine  law  or  the 

Friehe  and  Schwarz  relation,  where  m = 0.5.  The  ratio  of  i 

improvement  over  either  relation  is  smaller  for  the  very 
short  sensors  {t/d  * 100)  than  for  the  long  sensors 

(Z/d  * 1000) . The  trend  of  the  yaw  parameters  with  either  i 

sensor  Z/d  or  Reynolds  number  does  not  seem  to  approach  the 

cosine  law  but  rather  it  seems  to  overshoot  it  without  any  ; 

observed  asymptotic  values.  Values  of  the  transverse  ! 

sensitivity  coefficient  calculated  from  a yaw  sensitivity  j 

i 

function,  taking  into  account  deviations  from  the  cosine 
law,  yield  a higher  degree  of  accuracy  than  those  obtained 
directly  from  the  data.  Based  on  the  Reynolds  number  de- 
pendency of  the  yaw  coefficients,  one  should  use  different 
yaw  calibration  coefficients  for  different  velocities.  If 
a single  set  of  parameters  is  to  be  used,  the  accuracy 
achieved  by  a simple  relation  like  the  cosine  law  is  as  good 
as  one  may  expect.  However,  after  the  data  has  been  gath- 
ered, the  results  must  be  corrected  for  the  large  errors 
introduced  into  the  mean  and  rms  velocity  components.  Each 
new  sensor  mounted  on  a probe  must  be  calibrated,  even  if 
the  new  sensor  is  identical  to  the  one  being  replaced. 

The  parameters  of  the  yaw  relation  are  affected  by  the 
presence  of  unsteady  velocities  in  the  flowfield  and  de- 
pend on  the  frequency  of  these  fluctuations.  While  dynamic 
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yaw  calibration  of  the  probes  is  not  necessary,  the  mean 
yaw  sensitivity  of  the  probe  should  be  obtained  in  the  pres- 
ence of  some  velocity  fluctuations,  e.g.,  in  the  presence  of 
some  background  turbulence  of  similar  spectral  content  to 
the  flow  under  consideration. 
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CHAPTER  I 
INTRODUCTION 

Background 

Hot-wire  sensors  are  probably  the  most  commonly  used 
measuring  instrument  in  experimental  fluid  mechanics.  In 
conjunction  with  a constant  temperature  anemometer  unit, 
they  provide  unique  capabilities  which  to  date  are  only 
rivaled  by  the  laser  doppler  velocimeter.  Because  of  vari- 
ous limitations,  such  as  very  high  cost,  ambiguities  rela- 
ted to  the  scattering  particles  and  errors  introduced  by 
continuous  tracking  signal  processing  units,  hot-wire  anemo- 
meters will  most  likely  remain  a primary  instrument  for 
velocity  measurements.  The  hot-wire  does,  however,  have 
many  of  its  own  limitations,  some  of  which  will  always  re- 
main; e.g.,  need  for  the  sensor  and  its  support  to  be  locat- 
ed in  the  flowfield  and  the  nonlinearity  of  the  anemometer 
output.  Numerous  studies  have  been  carried  out  to  develop 
the  instrument  to  its  present  stage.  In  particular,  several 

of  these  studies  have  been  devoted  to  removing  some  of  the 
limitations . 

The  temperature  of  the  working  fluid  in  many  of  the 
filow  facilities  does  not  remain  constant  for  the  entire  du- 
ration of  the  experiment.  These  temperature  changes  may  be 
deliberate,  because  of  some  aspect  of  the  experiment,  or  in- 
cidental, G.g.,  due  to  the  energy  losses  from  the  fan  or 
pump  providing  the  flow  in  the  facility.  In  industrial 
applications  of  hot-wires,  such  temperature  variations  are 


2 


unavoidable.  Since  the  hot-wire  sensor  is  not  only  sensi- 
tive to  the  flow  velocity  but  also  to  the  fluid  temperature, 
density  and  viscosity,  such  variations  lead  in  many  occa- 
sions tc  substantial  errors  in  the  indicated  velocity.  A 
simple  and  inexpensive  method  of  temperature  compensation 
which  would  require  no  extensive  additional  calibration  of 
the  anemometer  would  be  ideal.  For  example,  an  acceptable 
method  would  be  one  which  can  be  accurately  used  once 
basic  properties  of  the  probes  are  known.  The  accuracy  of 
the  method  is  of  course  measured  by  the  errors  introduced 
during  temperature  variations.  Hence,  an  estimate  of  these 
errors  is  an  essential  part  of  the  development  and  documen- 
tation of  the  method. 

In  turbulent  or  complex  flowfields  various  arrays  of 
hot-wire  sensors  are  needed  in  order  to  resolve  the  various 
components  and  correlations  of  the  velocity  field.  Various 
types  have  been  used  over  the  last  three  or  four  decades, 
the  most  common  being  the  X-wire  probe.  In  order  for  any 
hot-wire  array  to  resolve  velocity  components  in  directions 
other  than  the  streamwise  direction  or  to  provide  informa- 
tion on  the  components  of  velocity  in  more  than  one  direc- 
tion, the  dependence  of  the  anemometer  output  on  the  yaw  or 
pitch  of  the  probe  must  be  known.  If  fluctuating  velocities 
or  their  correlations  are  to  be  measured  , this  dependence 
must  be  known  to  a very  high  degree  of  accuracy,  since  the 
derivatives  of  the  output  with  respect  to  the  yaw  or  pitch 
angle  are  required  in  order  to  infer  these  quantities  from 
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the  anemometer  signals. 

In  the  case  of  the  X-probe,  the  dependence  or  sensitivity 
of  the  output  signal  from  both  sensors  must  be  known  as  a 
function  of  the  yaw  angle,  measured  for  example  from  the 
axis  of  symmetry  of  the  X-wires.  When  the  mean  streamlines 
of  the  flow  are  not  parallel  to  this  axis  during  the  entire 
traverse  of  the  probe  in  the  flowfield,  the  mean  velocity 
vector  will  not  be  aligned  with  the  X-wire  axis  of  symmetry 
at  all  times.  In  such  cases  the  yaw  dependence  of  the  probe 
must  be  known  to  the  high  degree  of  accuracy  over  a wider 
range  of  angles.  Many  practical  and  complex  flowfields  are 
included  in  these  cases.  The  one  class  of  flows  which  moti- 
vated this  work  and  is  of  particular  interest  here  can  be 
described  in  general  as  "swirling  flows".  Such  flows  arise 
in  many  problems  as,  for  example,  can  be  seen  from  the  In- 
troduction to  Nagib's  work  [1].  Measurements  in  such  flows 
or  near  bluff  bodies,  in  particular  in  the  presence  of  tur- 
bulence, require  the  information  described  above. 

The  present  study  was  initiated  in  response  to  needs 
which  arose  in  connection  with  the  recent  work  by  Wigeland 
et  [2,3]  and  aims  at  providing  its  results  to 'future 

extensions  of  their  investigations.  To  obtain  the  measure- 
ments presented  in  their  papers  (2,3)  Wigeland,  Ahmed  and 
Nagib  utilized  X-wire  probes.  They  used  well  known  methods 


♦Numbers  in  brackets  refer  to  numbered  references 
in  the  Bibliography. 


4 


to  calibrate  the  probes  and  to  process  the  signal  with  the 
aid  of  analog  circuitry.  Several  such  methods,  which  rely 
on  different  approaches,  have  been  used  by  many  investiga- 
tors and  only  some  of  them  are  documented  in  the  literature. 
In  particular,  some  of  the  approaches  not  documented  in  the 
literature  rely  on  digital  processing  techniques  which  have 
been  recently  introduced  to  experimental  fluid  mechanics 
research,  but  seem  to  be  gaining  considerable  ground  over 
analog  signal  processing  systems.  At  I.I.T.  we  believe  that 
both  analog  and  digital  techniques  should  be  available  to 
the  experimenter  so  that  they  may  be  utilized  in  the  various 
stages  of  diagnostics  and  final  data  collection  of  the  ex- 
periment. Therefore,  in  our  opinion  a most  essential  ingre- 
dient required  in  an  approach  for  dealing  with  the  calibra- 
tion and  utilization  of  X-wire,  as  well  as  multiple  sensor 
probes, is  their  adaptability  to  both  signal  processing  tech- 
niques . 

Almost  all  of  the  methods  used  with  X-wires  rely  at  one 
stage  on  describing  the  yaw  calibration  data  by  some  ana- 
lytical function.  These  vary  from  simple  cosine  functions 
to  high  order  polynomials,  in  a few  cases  the  discrete  cali- 
bration data  are  used  in  the  method.  However,  these  suffer 
from  the  difficulties  involved  in  obtaining  derivatives  of 
the  anemometer  output  with  respect  to  yaw  angle. 

Many  of  these  methods  and  comparisons  between  them  are 
discussed  in  the  following  chapters.  In  the  remainder  of 
this  chapter,  some  of  the  results  in  the  literature  of 


importance  to  the  two  aspects  of  the  present  study,  i.e., 
yaw  and  temperature  dependence  of  hot-wires,  are  briefly 
discussed . 

Relevant  Literature 

Extensive  work  has  been  done  in  the  areas  of  tempera- 
ture and  yaw  dependence  of  hot-wires.  In  regard  to  the  tem- 
perature dependence,  semi-empirical  heat-transfer  laws  have 
been  developed  by  Collis  and  Williams  [4],  Davies  and 
Fisher  [5]  and  Bradbury  and  Castro  [6] . With  the  aid  of 
these  relations,  a number  of  schemes  to  measure  and/or  com- 
pensate for  the  ambient  temperature  variations  are  described 
both  in  the  literature  and  through  commercially  available 
compensators . 

Burchill  and  Jones  (7]  proposed  a scheme  in  which  the 
resistance  in  the  anemometer  bridge  that  controls  the  wire 
temperature  is  manually  adjusted  as  the  fluid  temperature 
changes,  the  adjustment  being  based  on  calibration  curves 
obtained  earlier  at  different  ambient  temperatures.  The 
amount  of  time  wasted  during  an  experiment  is  the  major  ob- 
jection to  this  scheme.  Other  techniques  developed  by 
Chevray  and  Tutu  (8J  and  similarly  by  Ali  (9]  involve  a 
complex  compensating  scheme  which  is  applied  to  the  signal 
from  the  bridge,  after  it  has  been  linearized,  through  the 
linearizer  constants.  Rose  [10]  determined  a method  by 
which  values  read  from  the  anemometer  were  corrected  through 
the  use  of  a correction  factor  which  was  a function  of  tem- 
perature. His  method  was  found  to  give  3%  errors  in 
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indicated  velocity  due  to  only  a 12®F  change  of  temperature. 
Methods  like  those  of  Chevray  and  Tutu  [8]  and  of  Ali  [9] 
are  particularly  suited  when  high  frequency  temperature 
changes  are  present  in  the  flowfield.  Their  methods  as 
well  as  Rose's  can  be  also  adapted  to  digital  signal  process- 
ing techniques. 

A number  of  compensating  units  are  sold  by  the  hot-wire 
anemometer  manufacturers.  However  cost  is  prohobitive  in 
many  cases  and  these  units  can  be  interfaced  only  with  the 
newer  anemometer  models.  Nagib  [1]  and  Tan-atichat  et  al. 

[11]  employed  a simple  compensating  scheme  that  consisted  of 
a temperature  sensing  probe  together  with  a resistive  net- 
work which  was  used  as  the  leg  in  the  bridge  of  a constant 
temperature  anemometer  opposite  the  velocity  sensor.  Using 
this  scheme  they  were  able  to  compensate  to  better  than  1% 
change  in  indicated  velocity  over  a 20“C  change  in  water 
temperature.  However  to  determine  the  values  of  resistances 
in  the  network  required  a long  and  tedious  iterative  pro- 
cedure. Nevertheless,  this  was  by  far  the  simplest  scheme  whicf 
attained  the  desired  accuracy,  as  long  as  the  frequency  re- 
sponse of  the  temperature  probe  was  higher  than  the  frequen- 
cy of  temperature  changes. 

In  regard  to  the  yaw  dependence,  Prandtl  [12] , 

Struminsky  [13]  and  Jones  [14]  determined  that  for  an  in- 
finite heated  cylinder  only  the  velocity  component  normal 
to  it  was  effective  in  the  cooling.  In  this  case  the  yaw 
sensitivity  function  is  given  by 
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f (a)  » cos  a (I-l) 

where  a Is  the  angle  between  the  incident  velocity  and  the 
normal  to  the  wire.  The  effective  velocity,  U ; is  there- 
fore  related  to  the  incident  velocity,  by: 

Ug  - f^(a)  (1-2) 

Hinze  [15] , Webster  [16]  and  Champagne  [17]  have  found 
that  finite  length  hot-wires  produce  deviations  from  the 
cosine  law  and  suggested  the  following  for  the  directional 
sensitivity, 

f^(a)  - (cos^o  + sin^a)**  (1-3) 

where  k is  a small  parameter  accounting  for  the  axial  com- 
ponent of  the  velocity.  All  of  them  report  that  k varies 
with  sensor  t/d  although  considerable  eunount  of  scatter  is 
evident  when  comparing  their  data.  Hinze  [15],  Bruun  [18], 
and  Kjellstrom  et  al. [19]  have  determined  that  k decreases 
with  increasing  velocity,  which  Horvatin  [20]  suggests  may 
be  due  to  the  axial  temperature  profile  of  the  wire  becoming 
more  isothermal  at  the  higher  velocities. 

Other  yaw  relations  have  been  proposed  by  Bruun  [18] , by 
Fujita  and  Kovasznay  [21],  and  by  Friehe  and  Schwarz  [22] 
where  in  their  case, 

f^(a)  - [1  - b(l  - cos‘‘o))^  (1-4) 
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The  Friehe  and  Schwarz  relation  appears  to  fit  most  yaw 
calibration  data  to  a greater  degree  of  accuracy  than  any  of 
the  previous  yaw  relations  examined.  However,  their  data 
122]  indicate  that  hot-films  which  have  small  f/d's  seem  to 
deviate  less  from  the  cosine  law  as  compared  to  hot  wires 
which  have  much  larger  Z/d’s. 

In  an  attempt  to  explain  these  trends,  and  to  examine 
the  conditions  in  water,  Roberts  et  a_l.  [23,24]  generalized 
the  Friehe  and  Schwarz  relation  to 

f^(a)  = [1  -b{l  - cos'"  o)]^  (1-5) 

and  investigated  the  yaw  dependence  of  commercially  available 
hot-film  probes  in  water.  This  generalization  was  based  on 
observed  trends  in  the  data  which  indicated  that  a combina- 
tion of  the  behaviors  depicted  in  Figures  1 and  2 is  more 
suitable  than  either  one  separately.  It  is  quite  clear  from 
these  two  figures  that  the  behavior  in  the  two  cases  is  sig- 
nificantly different  at  large  yaw  angles.  It  should  also  be 
noted  that  the  behavior  of  Figure  1 is  that  of  the  Friehe 
and  Schwarz  relation  while  Equation  (1-5)  combines  the 
trends  shown  in  both  figures.  Roberts  et  al^.  [23,24]  deter- 
mined that  the  accuracy  in  fitting  the  calibration  data 
using  this  relation,  which  is  called  the  present  relation 
in  the  following  chapters,  was  superior  to  the  Friehe  and 
Schwarz  relation.  However,  the  trends  they  sought  to  ex- 
plain still  remained  unresolved  after  their  studies. 

Champagne  et  al.  [25]  and  later  Friehe  and  Schwarz, 


each  using  their  own  yaw  sensitivity  relation,  have  deter- 
mined corrections  to  low  intensity  turbulence  measurements 
using  an  ideal  linearized  X-probe  for  the  case  whore  the 
mean  velocity  is  along  the  axis  of  symmetry  of  the  X-wires. 
These  corrections  have  not  been  determined  yet  in  complex 
flows  where  the  mean  velocity  cannot  be  aligned  with  the 
X-probe,  nor  has  a comparison  been  made  between  the  various 
turbulence  correlations  as  would  be  obtained  from  the  same 
data  using  different  yaw  relations. 

Morrison  et  al.  [26)  show  that  dynamic  calibration  of 
X-probes  yields  significantly  different  results  than  those 
obtained  by  static  calibration.  Broun  [17]  has  demonstrated 
that  part  of  this  difference  is  caused  by  poor  approximation 
of  the  static  calibration.  Bruun's  data  however,  reveals 
that  a difference  does  still  exist  which  he  attributes  to 
experimental  error. 

From  the  literature  cited  above  it  seems  that  there  is 
considerable  disagreement  and  uncertainty  in  the  interpre- 
tation of  signals  from  and  use  of  X-wire  probes.  This  study 
aims  to  clarify  some  of  these  problems. 


Objectives 

The  objective  of  this  work  can  be  divided  into  two 
groups;  one  relating  to  the  temperature  dependence  of  hot- 
wires , and  the  other  to  their  yaw  dependence. 

In  regard  to  the  temperature  dependence,  the  objective 
is  to  analyze  the  compensating  scheme  proposed  by  Nagib  [1] 
to  obtain,  whenever  possible,  exact  solutions  of  the  values 
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of  the  resistances  in  the  compensating  network,  and  to  de- 
rive estimates  of  the  errors  introduced  by  the  technique 
during  typical  temperature  variations.  It  is  considered 
desirable  to  be  able  to  perform  all  of  these  knowing  only 
the  temperature  dependence  of  the  velocity  sensor  and  com- 
pensating probes  and  their  resistances,  which  are  usually 
specified  by  the  manufacturer.  If  these  objectives  are 
accomplished  the  desired  compensation  can  be  achieved  and 
the  expected  errors  estimated  a priori  to  making  any  flow 
measurements  and  without  the  need  for  temperature  calibra- 
tion. 

In  regard  to  the  yaw  dependence  of  hot-wires,  the  first 
objective  is  to  do  a careful  study  of  the  effect  of  various 
parameters  on  the  yaw  sensitivity  as  described  by  the  cosine 
law,  the  Friehe  and  Schwarz  relation  and  the  present  rela- 
tion. Theee  parameters  include  Reynolds  number,  probe 
misalignment  and  length  to  diameter  ratio  of  the  sensor,  l/d. 
Two  possible  methods  of  varying  the  sensor  t/d  were  con- 
sidered during  the  planning  stages.  The  first  is  based  on 
varying  the  active  portion  of  the  sensor  while  keeping  the 
ratio  between  the  prong  spacing  and  the  active  length  con- 
stant. Another  method  is  to  keep  the  prong  spacing  constant 
while  changing  the  sensor  t/d.  The  latter  was  selected  for 
several  reasons,  although  the  former  appeared  attractive  for 
other  reasons.  While  the  length  of  the  active  sensor  is 
easily  changed  during  fabrication  the  prong  spacing  is 
usually  not  readily  altered.  The  separation  of  the  prongs 
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is  often  dictated  by  the  size  of  the  experimental  setup  and 
requirements  of  minimum  probe  interference.  On  the  other 
hand  the  length  of  the  active  sensor  may  be  varied  to  suit 
the  flow  conditions,  including  the  scales  of  the  velocity 
gradients  and  turbulence  field. 

The  second  objective  is  to  extend  the  analysis  of 
Champagne  [17]  and  Friehe  and  Schwarz  [22]  to  cases  when  an 
X-probe,  with  sensors  having  different  yaw  sensitivities, 
is  used  in  flows  where  the  mean  velocity  is  not  along  the 
line  of  symmetry  of  the  X-wires.  In  conjunction  with  this, 
a comparison  of  the  differences  between  the  values  of  the 
turbulence  velocity  correlations  computed  from  the  cosine  law, 
Friehe  and  Schwarz  and  the  present  relation  can  be  made. 

The  third  objective  is  to  examine  the  effect  of  free- 
stream  fluctuations  on  the  yaw  sensitivity  of  hot-wires. 

Since  hot  wires  are  calibrated  for  yaw  sensitivity  under 
static  conditions  with  very  low  turbulence  levels,  i.e. 
calibration  conditions,  operating  them  under  conditions 
where  larger  fluctuations  are  present  may  have  an  effect  on 
the  various  sensitivities. 
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CHAPTER  II 

ANALYSIS  OF  TEMPERATURE  COMPENSATING  CIRCUITS 
FOR  USE  WITH  HOT-WIRES  OR  HOT-FILMS 

A number  of  temperature  compensating  schemes  have  been 
proposed  by  the  hot-film  manufacturers.  Several  of  these 
schemes  depend  on  a temperature  sensing  probe  which  can  be 
used  in  the  anemometer  bridge  shown  in  Figure  3 in  place  of 
the  bridge  resistance  R^.  The  ideal  probe  for  this  purpose 
should  have  the  following  properties; 

1.  Its  resistance  R^  should  be  given  by 

^ = r^MR^  = MRjj 

where  M is  the  bridge  ratio  and  r is  the  overheat  ratio. 

o 

2.  The  dependence  of  R.J,  on  temperature  should  be  iden- 
tical to  that  of  the  velocity  probe  operating  resistance  Rjj. 

3.  The  size  of  the  probe  should  be  sufficiently  large 
so  as  not  to  be  heated  by  the  bridge  current. 

The  latter  requirement  assures  that  the  probe  resistance  is 
independent  of  the  fluid  velocity. 

Such  a probe  is  of  course  almost  impossible  to  manu- 
facture and  even  if  it  could  be  selected  from  a large 
statistical  sample  it  could  only  be  used  with  the  one 
matched  velocity  probe.  Fortunately  most  of  the  hot -film 
and  hot-v/ire  probes  furnished  by  a manufacturer  are  quite 
similar  in  properties  with  their  resistances  falling  within 
a limited  range.  A temperature  compensating  circuit 
similar  to  those  outlined  in  Figure  3 can  therefore  be 
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used,  at  least  over  some  range  of  ambient  temperature 
variation.  A complete  analysis  for  the  determination  of 
the  series  and  parallel  resistances  of  compensating 
circuits  1,  2 and  3 are  presented. 

Analysis 

Based  on  the  bridge  balance  condition  of  the  anemometer 
circuit  in  Figure  3, 


R1/R2  = Rh/R3  = M (II-l) 

where  M is  the  bridge  ratio  and  R^  is  the  equivalent 
resistance  of  the  temperature  compensating  network 

t 

Considering  only  the  first  coefficient  of  the 
temperature  dependence  of  R^  and  R,^  (i.e.,  ignoring  second 
and  higher  order  terms  of  temperature  changes)  we  write 


dR^ 

dT  “ ^o“t  (II ”2) 

and 

af^  = «co“c 

Or  that 

Rt  = 

and 

Rc  - '"•=> 

where  R^^  and  R^^  are  the  resistances  of  the  temperature 
compensating  probe  and  the  cold  resistance  of  the  velocity 
probe,  respectively,  when  the  6unbient  temperature,  T,  is 
equal  to  the  probe  reference  temperature  T^. 
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Assuming  that  is  independent  of  temperature,  vrtiere 


2 — and  a>0 

r Oc  r 


(II-6) 


we  define  a probe  resistance  ratio  as  follows 


(II-7) 


The  resistance  ratio  6 is  dependent  on  temperature  and  its 
value  at  the  original  reference  temperature  T^,  (usually 


room  temperature)  is 


Rm 

° Rc  t*=t 


(II-8) 


We  denote  the  operating  resistance  of  the  velocity  sensor 
by  Rjj,  i.e., 

r^(T)>l  (II-9) 

where  r^  is  the  overheat  ratio.  Since  is  related  to  the 
temperature  compensating  circuit  resistance  through  the 
bridge  relation,  it  is  also  a function  of  the  ambient 
temperature.  First,  we  consider  compensating  circuit  1, 
where 


RpRip 

From  the  bridge  balance  condition,  we  write 


p , m[r  + 

^ Rp+R^ 


(II-IO) 


(II-ll) 


Using  Equations  (II -5)  and  (II-7)  and  defining  non- 
dimensional  series  and  parallel  resistances,  and 
respectively,  where 
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and  Sj.>0  (II  ■‘12) 


and  Pjj.>0  (11-13) 

we  obtain  the  nondimensional  bridge  balance  condition, 


The  governing  equation  for  the  velocity  sensor  (King ' s 
Law)  can  be  written  as 

= A + BU"  (11-15) 

where  U is  the  fluid  velocity  and  Ip  is  the  current  through 

the  probe;  A,  B and  n are  nearly  constant.  Based  on  this 

relation  and  since  for  constant  temperature  operation  the 

anemometer  servo -amplifier  maintains  the  probe  temperature 

and  hence  its  resistance  R constant,  one  may  conclude  that 

n 

two  possible  modes  may  be  useful  for  keeping  the  output 
voltage  constant,  i.e.,  independent  of  ambient  temperature. 
One  can  either  maintain  a constant  overheat  ratio  or 
maintain  a constant  resistance  difference. 

Another  possible  mode  exists  if  we  rewrite  King's  Law 
in  terms  of  the  bridge  output  voltage  and  determine  the 
condition  necessary  to  keep  the  output  independent  of 
ambient  temperature. 

Case  A;  Constant  Overheat  Ratio.  In  this  case  we  set 

r_  ■ constant 
o 


^”^C 


and 


S = 

®r  " R„ 


2p 
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Since 


M = = -O-g 

R3  R3 


(11-16) 


and  the  anemometer  bridge  ratio  M is  constant,  then 


Rc(T) 

; • rtttt  “ 


and 


dT 


= 0 


so  that  using  Equations  (11-2,3,4)  we  obtain 


dR, 

dr" 


Co  C o 
M 


(11-17) 


Case  B:  Constant  Overheat  Resistance  Difference . In 
this  case  v/e  set 

Rj^  - R^  * constant 

so  that 


Using  Equations  (II-l)  and  (II-4) , we  obtain 


dR3 


^Co“c 

M 


(11-18) 


Case  C:  Con stant  Bridge  Output  Voltage . In  this  case 
referring  to  Figure  3,  and  using  the  bridge  balance 
condition,  the  non-linear  governing  Equation  (11-15)  may  be 
rewritten  as 


(Rl+Rjj)  ^ (A+BU") 


(11-19) 


(note  that  if  the  ambient  temperature  remains  constant, 
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this  reduces  to  the  form  = A'  + B’u")  . Assuming  that 
the  coefficients  A,  B,  and  n are  insensitive  to  temperature 
variations  relative  to  the  changes  in  and  with 
ambient  temperature  (the  validity  of  this  assumption  is 
examined  in  Appendix  B) , we  set 


(11-20) 


Using  Equations  (II-l),  (II-4)  and  (11-19),  we  obtain 


RcoV^'^^ 

M 


(11-21) 


where 


T (T) 


2 + 


^i-1 


(11-22) 


Let  us  now  consider  one  final  case.  In  this  case  we 
determine  the  condition  necessary  to  )teep  the  output  voltage 
constant  as  the  ambient  temperature  varies  if  the 
anemometer  output  is  linearized. 

Case  D:  Constant  Linearized  Output  Voltage.  From 
Equation  (11-19)  we  can  write 


Kj^(T)U 


(11-23) 


vhere  is  the  linearized  output  voltage  and  the  lineari- 
zation constant  is  a function  of  T given  by 


Kj  (T) 


B 


(Rl+R,,) 


1/n 


(11-24) 


For  the  output  to  be  independent  of  temperature 
variations,  we  set 
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that  is. 


(11-25) 


dK, 

= 0 (11-26) 


Using  Equations  (II-l) , (II-4) , (11-24)  and  (11-26)  we 
obtain 

dRi  (T) 

ar  M 


Comparing  this  with  Equation  (11-21) , we  find  that 
temperature  compensation  does  not  depend  upon  whether  the 
anemometer  output  voltage  is  linearized  or  not. 

We  now  combine  all  of  the  above  cases  into  a single 
equation  which  specifies  the  functional  dependence  of  the 
compensating  circuit  with  temperature.  This  equation  can 
be  written  as 


dRs 

aT“ 


^Co“c 

Mn 


(11-27) 


where 


for  Case  A 


= \ 1 for  Case  B 


(11-28) 


for  Cases  C and  D 


HtT 

The  solution  of  Equation  (11-27)  along  with  the  non- 


dimensional  bridge  relation  for  each  circuit  will  yield 
values  of  the  series  and  parallel  resistances  for  each 
circuit. 


Compensating  Circuit  1.  Substitution  of 
Equation  (II -10)  into  Equation  (11-27)  and  using 
Equations  (11-12)  and  (11-13)  leads  to  the  solution  for 
the  parallel  resistance  directly 

Pj.  = (/o^HtTB^-I)  (11-29) 


or 


Rp  = lR^(/arMn  Pq-I)  (11-30) 

o 

then  substitution  of  Equation  (11-29)  into  Equation  (11-14) 
leads  to 


= 


rr. 


6M  7a  M'n'  B 
r o 


T=T 


(11-31) 


o 

Here,  the  values  of  Rg  and  are  evaluated  at  the  initial 
reference  temperature  T^. 

Compensating  Circuit  2.  Referring  to  Figure  3,  we 
obtain 


Use  of  Equations  (II-l) , (11-12)  and  (11-13) 
non-dimensional  bridge  balance  condition, 


(11-32) 
leads  to  the 


Sr  + 1 - Prfr^^^r"^^^  ° (11-33) 

L o 3 

Substitution  of  Equation  (11-32)  into  Equation  (11-27) 
yields  a non-dimensional  form 


“ Pr^''®r”'^  ^o  ^ (11-34) 

The  solution  of  Equations  (11-33)  and  (11-34)  for  the 
series  and  parallel  resistances  results  in  the  resistances 
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eM(/a^Mn  e^-1 


jt=t 


(11-35) 


and 


gM 


✓a  Mn  g - 1 


(11-36) 


T»T 


Operating  charts  for  the  selection  of  the  series  and 
parallel  resistances  for  circuits  1 and  2 are  also  plotted 
in  Figures  4 and  5 , respectively . These  charts,  or  the 
equations  given  above,  can  be  used  by  the  reader  to  design 
or  adjust  either  of  the  tvro  compensating  circuits.  One 
needs  only  to  know  the  properties  of  the  anemometer  bridge 
and  the  velocity  and  temperature  probes,  i.e.,  their 
resistances  and  their  temperature  dependence  coefficients, 
in  order  to  achieve  the  compensation,  without  any  need  for 
temperature  calibration. 

Compensating  Circuit  Circuit  3 is  slightly 
different  than  the  previous  two  circuits  presented  in  that 
there  are  four  unknown  resistances  to  determine.  Since  we 
are  working  with  two  relationships,  it  would  seem  that  at 
this  point  two  of  the  resistances  can  be  selected 
arbitrarily.  It  will  be  shown  later  that  some  constraints 
on  the  other  two  resistances  exist.  However,  initially  we 
will  proceed  as  if  there  are  no  constraints  in  order  to 
generate  a family  of  solutions. 
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From  Figure  3,  the  equivalent  resistance  is 


3 ^^P2''’^S2^  * 

Using  the  previous  methods,  the  non-dimensional  bridge 
relation  becomes 

<^r2*=r2>  * 

- 1;' <®rl*''rl^l”’r2Sr2  * ''^1^2'S^1*»  1 * » 

Substitution  of  Equation  (11-37)  into  Equation  (11-27) 
yields  the  non-dimensional  form 


(11-37) 


(Sri+Pri^l) 

(11-39) 

the  resistances  will  be  arbitrary  we 

define 

®r2  = ®'=rl  ® 

(11-40) 

'“r2  = '’■''rl 

(11-41) 

b'  are  constants  to  be  determined. 

define  the  following  quantities, 

= Srl 

(II-42a) 

^2  = '’rl‘Srl*l> 

(II-42b) 

o 

(II-42C) 

o 

(II-42d) 

Substitution  of  Equations  (11-40)  and  (11-41)  into 
Equation  (11-38)  along  with  Equation  (11-42)  leads  to 
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a'Zj^  + Zj 

a'Z,  + Z. 
3 4 


(11-43) 


Substitution  of  Equation  (11-40)  and  Equation  (11-41) 


into  Equation  (11-39)  leads  to 


a'Zj^  + Z2 

ZT 


(11-44) 


where 


“is  = '’rl  - '■rl 


Solving  for  a'  and  b'  in  Equations  (11-43)  and  (11-44), 


we  obtain 


S - "4 


(11-45) 


Z^(Zg-Z^)  -»•  Z2Z3 

Z5Z3 


(11-46) 


Since  a'>0  we  find  from  Equation  (11-45)  that 


Thir.  leads  to 


0 1 ^ /cXj.Mn  6q  - 1 


(11-47) 


The  condition  on  the  coefficient  b'  is  that  b'>0  which 


implies  that 


and  hence 


^1  " 


Z5  > 0 


JjlH— 

/o^Mn  - 1 


(11-48) 


f ^ 
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Having  selected  and  based  on  Equations  (II -47) 
and  (11-48),  P^2  ^£-2  determined  from 

Equations  (11-40)  and  (11-41). 

Putting  the  solution  in  dimensional  form  vje  have  the 
following : 


(1)  Select  Rgj^  such  that 


0 i ■'si 


i [*t(®  •'V"  6,  (11-49) 


(2)  Select  Rpj^  such  that 

R„.  + R 


SI 


^1  ' 

/a  Mn  - ij 

^ r o 

1 * i 

o 

(3)  Determine  Rg2 

from 

P &E 

r ^1 

rVo’'“r”^  ^o 

^S2 

BM 

SI  ^ 

(11-50) 


(11-51) 


T=T. 


and  (4)  Rp^  is  given  by 


'P2 


BMlRp^/a^Mn  - (Rgi+Rpi+V^ 


(11-52) 


T»T, 


Selection  of  A Compensating  Probe 

Ideally,  it  would  be  desirable  to  have  one  temperature 
compensating  probe  for  all  velocity  probes.  For  typical 
probes  this  is  not  always  possible.  However,  Icnowing  the 
range  of  resistances  of  the  velocity  probes  and  their 
temperature  coefficients,  a temperature  probe  may  be 
carefully  selected  for  use  in  most  of  the  cases.  Based  on 
the  solution  of  compensating  circuits  1,  2 and  3,  and 


t 


requiring  that  a real  solution  exist > the  follovdng 
inequality  involving  the  parameters  of  both  probes  must 
hold: 
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Co 

MoCnTr’ 


Co  or  o 


(11-53) 


With  a knowledge  of  the  properties  of  available  velocity 
sensors,  operating  overheat  ratios,  anemometer  bridge  ratio 
and  operation  mode  of  the  compensating  circuit,  a suitable 
temperature  compensator  may  either  be  fabricated  or 
purchased  from  a manufacturer  based  on  this  inequality. 

Besides  this  restriction,  the  compensating  probe  must 
not  be  influenced  by  the  electrical  current  passing  through 
it.  This  results  in  temperature  probes  sensitive  to 
velocity,  which  is  not  only  an  undesirable  feature  but 
also  affects  the  calibration  curve  of  the  velocity  probe  as 
discussed  by  Tan-atichat  et  al. [II] . The  compensating 
circuit  also  has  a frequency  response  to  temperature 
variations  v^ich  depends  on  the  thermal  time  constant  of 


the  temperature  probe  (up  to  several  cycles  per  second 
can  be  obtained  using  commercially  available  probes) . 

Error  Estimates 

We  do  not  expect  "perfect"  compensation  from  the 
proposed  circuit  configurations.  This  is  true  since 
Equation  (11-27)  prescribes  the  dependence  of  on  ambient 
temperature  for  ideal  compensation,  while  the  introduction 
of  a proposed  circuit  and  a temperature  probe  forces  a set 
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function  of  temperature  on  . Therefore,  the  error  will 
be  related  to  the  difference  between  the  two. 

Case  C:  Non-Linear  Output  Voltage.  As  before. 


% - 


(A+BU*^) 


Then,  with  the  assumption  that  A,  B and  n are  not  functions 
of  temperature, 

3E, 

'b 


, 1 
e7  3t“  ~ 


-1 


R<^'i 

-1 

r 9R-3 

r R _ R - R_ 

1— ^ 

2-  C'+ 

RhJ 

3T 

1 ^ 

(11-54) 


If  Rj  behaves  ideally  then  the  terms  in  the  brackets  vgould 
vanish  and  Equation  (11-27)  vould  arise.  Since  this  is 
not  necessarily  true,  we  integrate  Equation  (11-54)  from 
the  initial  reference  temperature  T^  to  some  ambient 


temperature  T to  obtain 


Eb(T)  * Ej^(TQ)e 


1/2 


(11-55) 


v4iere 


M 


I = 


3R. 

dr" 


R^  RiRq 

■ 


“c^Co 


1+ 


MR. 


1- 


MR3  J 


dT 

(11-56) 


I is  a function  of  R^ (T)  which  is  a function  of 


Since 

the  particular  compensating  circuit  used,  the  magnitude  of 
the  change  in  the  output  voltage  with  temperature  will  be 
a function  of  the  circuit  configuration  employed. 


d 


A Sample  Estimate.  We  now  consider  a typical  hot-wire 
and  compensator  combination  to  obtain  an  indication  of  the 
quality  of  the  compensation  using  the  proposed  scheme. 
Assume  we  are  using  a tungsten  hot-wire  whose  resistance 
at  room  temperature  is  lOfJ  and  whose  temperature 
coefficient,  o^,  is  0.0048/®C.  From  Equation  (11-53)  we 
determine  that  a compensator  made  from  nickel  whose  resist- 
ance at  room  temperature  is  150  and  whose  temperature 
coefficient,  a^,  is  0.0067/®C  is  suitable. 

Using  the  analysis  described,  the  values  of  the 


J 


resistances  for  circuit  1 are 


I 


» 


\ 


I 

1 


X 
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Rp  = 49.68fi 
Rg  = 6.48n, 

and  for  compensating  circuit  2 are 

Rp  = 77.6in 

Rg  = 8.44  ft. 

Using  these  values  we  can  determine  R2(T)  and  dR^/dT  and 
thus  evaluate  Equations  (11-55)  and  (11-59)  for  different 
eunbient  temperature  changes  for  circuits  1 and  2.  The 
results  are  shovm  in  Table  1. 

The  estimates  presented  in  Table  1 indicate  that  the 
proposed  compensating  circuits  will  tend  to  slightly  under - 
compensate  for  eunbient  temperature  changes.  The  use  of 
compensating  circuit  1 with  a temperature  change  of  40®C 
produces  a 2%  change  in  the  indicated  velocity,  while 
utilizing  circuit  2 would  reduce  this  error  to  1%  for  the 
same  temperature  change.  Over  the  range  of  temperature 
changes  exaunined  the  error  is  approximately  proportional 
to  the  square  of  the  temperature  difference. 

Next,  we  ask  the  following  question:  by  changing  from 
circuit  1 to  circuit  2,  the  error  was  reduced  by  a factor 
of  2;  can  we  then  minimize  the  error  by  going  to  circuit  3? 

At  the  reference  temperature,  the  compensating  circuit 
must  satisfy  constraint  equations  as  derived  previously. 

It  must  satisfy  the  bridge  relation,  i.e., 

- 0 , (11-60) 


«3“ 


^o«C 


IT 
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Table  1.  Sample  Error  Estimates  for  Typical  Velocity  and 
Temperature  Probes  in  Circuits  1 and  2. 


T - T^(“C) 


0 

10 

20 

30 

40 

Rj,(T)  (fi) 

10.00 

10.48 

10.96 

11.44 

11.92 

R^(T)  (fl) 

15.00 

16.01 

17.01 

18.02 

19.02 

Circuit 

1 

R3(T)  (fi) 

18.00 

18.59 

19.15 

19.70 

20.23 

dR^ 

35^ 

fn  1 

0.0593 

0.0575 

0.0557 

0.0541 

0.0525 

E^(T)-E^(To) 

(%) 

n 

-0.03 

-0.11 

-0.25 

-0.44 

U 

(%) 

n* 

-0.13* 

-0.52* 

-1.17* 

-2.05* 

u 

Circuit 

2 

R3 (T)  (0) 

18.00 

18.59 

19.17 

19.73 

20.28 

dR3 

af“ 

h] 

0.0593 

0.0581 

0.0570 

0.0559 

0.0548 

EbfTI-Eb'V 

(%) 

0 

-0.01 

-0.05 

-0.12 

-0.25 

0 #ffl  1 

(%) 

0* 

-0.06* 

-0.23* 

-0.53* 

-0.94* 

n - 0.43 
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and  the  temperature  relation,  Equation  (11-21),  i.e.. 


dT 


M 


0 , 


(11-61) 


T=T. 


for  the  constant  output  voltage  modes.  The  error  involved 
is  related  to  I which  is  defined  by  Equation  (11-56) . 

Since  in  circuits  1 and  2 there  are  only  two  adjustable 
resistances,  I is  fixed,  i.e.,  the  error  is  unique. 
Circuit  3,  however,  has  four  undetermined  resistances  so 
that  I can  be  minimized  subject  to  the  two  constraints. 

For  a prescribed  temperature  difference, 
which  one  wishes  to  compensate  over. 


I ^^^S1'^S2'^P1'^P2^ 


(11-62) 


where 


f(Rsi-Rs2'«Pl''>P2>  “ 


and  where  I*  is  a )cnown  function  of  temperature,  namely, 
MdR,  r R,  R.sT  r R,  ■] 


I*  s 


MdR,  r Rf,  R,R„1  r I 

MR,  “IT?  *^Co“c  MR,  J 

• *5  W J W J J 


(11-64) 


F^J 


A necessary  condition  that  I has  a minimum  value  is  that 


dl  - II—  dR 


wr:  ""'si  aLo  *^^82  ^*S»1  ?R^  ^^2  “ ° 


'SI  ” S2 

subject  to  two  constraints 


(11-65) 


and 
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where  (j)^  and  defined  in  Equations  (11-60)  and 

(11-61) . Using  the  method  of  Lagrange  multipliers,  we 
obtain  six  equations  in  six  unknowns,  ^S2'  *Vl'  ^2' 

and  Xj#  vdiere  and  Xj  are  the  Lagrange  multipliers. 
The  six  equations  are: 


3 


SI 


[T,  X,3R-  X23^R, 

J.  ^ 


SI 


SI 


(II-66a) 


. fT,  X,3R-  X53^R. 

^ I ^ I*  dT  4*  ^ ^ «f  ^ ^ 

^ It 

o 


(II-66b) 


fT,  X,3R-  R, 

I It 


(II-66c) 


kii 

o 


(II-66d) 


r R^i 
n - o C 
.3  M 


T-T 


(II-66e) 


3R3 

W 


M 


T-T. 


(II-66f) 


In  the  above  equations  R^  is  given  by  Equation  (11-37) . 
Solution  of  the  above  system  of  equations  will  lead  to 
values  of  resistances  which  make  circuit  3 the  optimized 
circuit  configuration.  The  solution  of  this  system  of 
equations  has  not  yet  been  attempted.  However,  in  principle 


t 


31 


f they  can  be  solved  numerically. 

The  performance  of  the  constant  output  voltage  mode 
was  based  upon  the  assumption  that  the  coefficients  in 
T King's  Law  were  independent  of  temperature.  From  the 

experimental  results  in  Appendix  B,  we  see  that  all  of  the 
coefficients  are  weak  functions  of  ambient  temperature. 
However,  data  for  the  different  operating  modes  indicates 
that  the  anemometer  output  voltage  drift  is  independent  of 
the  fluid  velocity  for  measurements  in  air.  Therefore  the 
variations  of  B and  n with  temperature  may  be  ignored. 

The  analysis  can  be  modified  to  take  into  account  the 
variation  of  the  coefficient  A with  temperature.  However, 

' the  variation  of  A would  have  to  either  be  assumed  or 

experimentally  determined.  In  this  case,  the  compensating 
circuits  could  not  be  set  a priori  to  making  any  velocity 
measurements  and  the  experimenter  would  have  to  calibrate 
his  probe  for  both  velocity  and  temperature.  The  error 
introduced,  however,  by  making  this  assumption  should  be 
% small  and  will  be  experinientally  determined  in  future 

extensions  of  the  present  investigation. 

X 

t 
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CHAPTER  III 

RESPONSE  EQUATIONS  OF  INCLINED  CYLINDRICAL  SENSORS 
IN  STEADY  AND  FLUCTUATING  VELOCITY  FIELDS 


In  this  chapter,  the  X-probe  array  of  Figure  6 is  consid- 
ered to  be  Inunersed  in  a flow  field  described  at  any  instant 
of  tine  by  an  instantaneous  velocity  Q^.  This  flow  field  can 

be  deconposed  into  a mean  velocity  Q which  is  inclined  at 

s 

an  angle  y with  respect  to  the  x-axis,  and  the  tine  dependent 
velocities  g^  and  aligned  with  the  mean  streamline  and  the 
direction  normal  to  it,  respectively.  Throughout  this 
report  the  effect  of  fluctuations 'which  are  normal  to 
both  and  q^  is  not  considered.  The  response  of  the 
X-probe  array  to  this  low  turbulence  intensity  velocity  field 
will  be  determined  utilizing  various  relations  for  the  yaw 
sensitivity  to  determine  the  differences  between  them  in 
calculating  the  flow  field  characteristics. 


Time-Mean  Velocity  Field 

Only  the  effect  of  Q_  on  the  time  averaged  output  of 
the  X-probe  is  examined  here.  As  will  be  shown  later  this 
is  only  true  for  cases  of  low  turbulence  intensity.  Refer- 
ring to  Figure  6,  the  mean  effective  cooling  velocity  that 
each  sensor  is  subjected  to  can  be  related  to  the  mean 
streamwise  velocity  through  the  equation 


for  sensor  1 and 


Ogf^Oj+Y) 


(III-2) 


for  sensor  2.  The  function  f (6)  is  determined  from  the 


yaw  sensitivity  relation  utilized.  The  mean  effective 


cooling  velocity  is  related  to  the  mean  output  voltage  of 


the  sensor  by 


®e,l  " i « 1,  2 


(III-3) 


Where  H(E^)  can  take  many  forms.  Two  of  these  forms  are 


given  by 


H(E^) 


(III-4a) 


- A T 

fci 

• »i 


(III-4b) 


Where  is  the  calibration  constant  for  a linearized  sensor 


and  A^/  and  n^  are  the  calibration  constants  for  a 


non-linearized  sensor.  In  each  case  calibration  is  carried 


out  when  the  sensor  is  normal  to  the  free-streaon  velocity. 


Cosine  Law.  The  angle  of  inclination  Y and  the  mean 


U and  V velocity  components,  aligned  with  the  x and  y axes, 


respectively,  are  to  be  extracted  from  Equations  (III-l) 


and  (III-2) . If  one  is  fortunate  enough  to  have  sensors 


which  obey  the  cosine  law,  then  Equations  (II 1-1)  and  (II 1-2) 


are  linear  in  U and  V and  lead  to 


(cot  + cot  ©2) 


(tan  + tan 


.ifH^ 

I sin 

■1  r»<V 

^cos 


H(E2)  1 
sIHTTj 


(III-5) 


h(E2) 


cos  e. 


(111-6) 


where 


-irv 


Y - tan  I u 


(III-7) 
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When  the  hot-wire  signals  are  linearized.  Equations  (III-5) 
and  (111-6)  represent  a simple  weighted  sum  and  difference 
technique  and  are  therefore,  very  attractive  for  use  with 
analog  processing. 

As  Champagne  et.  al.  [25]  have  shown,  nearly  all  sensors 
deviate  from  the  cosine  law  and  the  magnitude  of  this 
deviation  is  a function  of  probe  characteristics.  In  order 
to  correct  for  these  deviations  more  sophisticated  yaw 
sensitivity  relations  have  been  proposed.  These  relations 
introduce  nonlinearities  in  such  a manner  that  a simple 
weighted  sum  and  difference  technique  as  in  Equations  j.lI-5) 
and  (1II-6)  is  usually  not  feasable.  To  determine  the  mean 
velocity  components,  the  angle  of  mean  flow,  y,  may  first 
be  determined  from  an  iteration  scheme  using 


- y)  H(E^) 

f^(e2  + y)  HCEj) 


(III-8) 


The  mean  8tre2unwise  velocity  can  then  be  determined  using 
Equation  (III-l)  or  (III-2)  with  Equation  (III-3) . A 
simple  decomposition  of  this  mean  velocity  then  yields  the 
mean  U and  V components  of  velocity. 

In  general,  proper  schemes  for  determining  the  mean 
components  when  talcing  into  account  deviations  from  the 
cosine  law  do  not  readily  lend  themselves  to  analog 
processing  techniques.  How  large  an  error  then  does  one  make 
in  the  determination  of  the  mean  velocity  if  one  assumes  the 
cosine  law  to  take  advantage  of  the  simplicity  previously 
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described?  For  ease  of  interpreting  the  results  we  assume 

that  “ ©2  “ that  both  sensors  have  identical  yaw 

sensitivity.  Let  y,  Q , U and  V be  quantities  determined 

utilizing  a yaw  sensitivity  function  taking  into  account 

the  deviations  from  the  cosine  law,  and  y , (Q  ) , U , V 

c s c c c 

be  those  quantities  calculated  from  the  seune  data  utilizing 
the  cosine  law.  It  can  then  be  shown  that 


. tan-l  - f^45.v)' 

° Lf'^(45«-Y)  + f'^(45+Y). 


(III-9) 


(IIl-lO) 


(III-ll) 


where 

- / If^(45-Y)]^  + [f^(45+Y)]^ 

^8  (III-12) 

Examining  Equations  (IZI-9)  through  (III-12) , one  observes 

that  as  the  response  of  the  sensors  deviates  from  the  cosine 

law  the  magnitude  of  the  errors  introdu-ed  increases.  These 

errors  will  be  significant  particularly  in  the  angle  of  the 

mean  flow  and  the  V component. 

It  was  pointed  out  in  Chapter  I that  there  are  many 
yaw  sensitivity  functions  in  the  literature  which  take  into 
account  the  deviations  from  the  cosine  law.  For  mean  velocity 
measurements  when  Y>0,  there  is  no  significant  difference 
in  accuracy  among  them.  In  the  next  few  chapters  it  will  be 
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shown  that  this  is  definitely  not  the  case  for  time  depen-* 
dent  velocity  measurments.  It  will  also  be  demonstrated  that 
the  present  relation  yields  the  highest  accuracy  in  all 
cases.  Thus,  for  the  prupose  of  examining  the  errors  intro- 
duced when  using  Equations  (III-5)  and  (III-6),  we  assume 
here  that  f(a)  is  given  either  by  the  present  relation  or 
the  Friehe  and  Schwarz  relation. 

Since  we  will  be  examining  the  difference  in  a quantity 
X as  calculated  by  the  present  relation  and  by  the  cosine 
law,  we  define  a relative  error  in  x by 


_E £ 


(III-13) 


where  p and  c denote  the  present  and  cosine  relations, 
respectively.  Similarly,  the  difference  in  the  saune  quantity 
X as  calculated  by  the  Friehe  and  Schwarz  relation  and  by 
the  cosine  law  can  be  written  as 


(III-14) 


lm*0. 5 

For  comparison,  the  difference  between  the  present  relation 
and  the  Friehe  and  Schwarz  relation  is  denoted  by 

An  ideal  X-probe  whose  sensors  have  identical  yaw 
response  and  are  inclined  at  ±45*  to  the  x-axis  will  now  be 
considered.  When  the  mean  streamwise  velocity  is  aligned 
with  the  x-axis,  using  Equations  (III-IO) , (111-12)  and 
(111-13)  it  can  be  shown  that 
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U 


- 1 


Y*0 


(III-15) 


and 


■U 


' Y»0 


rl-.159  r 

L l-b(l-2"”‘/^)  J 


- 1 


(III-16) 


For  all  cases  within  the  scope  of  this  work  it  was  found 
that  £u1y  0^  whenever  Y®0  the  cosine  law  is  quite 

adequate  in  determining  the  mean  velocity  irrespective  of  the 
sensor  being  used.  However,  the  magnitude  of  this  error  will 
be  shown  to  depend  on  a number  of  pareuneters. 

For  the  cases  where  Y^o,  sensors  whose  response  is 
close  to  the  cosine  law  are  first  considered.  Next,  sensors 
whose  response  is  far  from  the  cosine  law  will  be  discussed. 
As  an  example  of  the  former  case,  we  assume  that  b«1.046  and 
that  m=0.487  (details  of  these  results  can  be  found  in 
Chapter  V) . Using  Equations  (III-9)  through  (III-12)  one 
then  finds  that  for  Yjc25‘* 

Zy  z 2.5% 


ey  = "2% 


Therefore,  for  sensors  that  do  not  show  a marked  deviation 
from  the  cosine  law,  mean  velocity  measurements  are  not 
significantly  altered  by  employing  a "more  sophisticated" 
yaw  relation  and  the  simple  sum  and  difference  technique  of 
the  cosine  law  is  quite  adequate. 

For  the  latter  case  we  take  b“0.534  and  m“0.868.  Using 
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these  pareuneters  in  Equations  (III-9)  through  (III-12)  one 
obtains  for  y=5* 

e = 20% 

Y 

Eu  - -5% 

Ev  “ 20% 

As  the  angle  of  the  miean  flow  increases  these  errors  also 
increase.  Thus,  for  sensors  which  show  a noticeable  devia- 
tion from  the  cosine  law,  the  simple  sum  and  difference 
technique  does  not  apply  for  cases  where  Y^O,  making  analog 
processing  extremely  difficult  if  not  impossible. 

Furthur  Generalization  of  Cosine  Law.  For  sensors  which 
exhibit  small  deviations  from  the  cosine  law  it  was  shown 
that  utilizing  the  cosine  law  leads  to  quite  accurate  results. 
An  approximate  method  which  may  result  in  even  more  accurate 
results  is  based  on  the  use  of  a weighted  sum  and  difference 
technique  utilizing  a "more  sophisticated"  yaw  relation. 

Here  it  is  assumed  that  the  mean  effective  cooling  velocity 
can  be  decomposed  in  the  manner  similar  to  the  cosine  law. 

and 

«e,2  * "ap*^<®2>  ' ®2>  (III-17b) 

2 

The  function  f ( 6^)  once  again  refers  to  the  yaw  relation 
utilized  and  the  subscript  ap  refers  to  quantities  calculated 
using  this  approximate  method.  Since  the  angle  of  inclina- 
tion of  the  sensors  is  known,  then  once  the  sensors  are 


calibrated  for  yaw  sensitivity,  the  weightling  coefficients 

are  determined.  Solving  Equations  (III-17)  for  and 

ap 

V one  obtains 
ap 


f (Sj) 


f (90““ej^) 


H(E2) 

f^(90»-e2) 

f^(90»-e2) 


(III-18a) 


"h(Ej^)  H(E2) 

f^(e^)  ?(e2)_ 

‘f^(90“-ej^)  f^(90“-f)2)" 

f'^(e3^)  f (62) 


(III-18b) 


The  errors  introduced  using  this  scheme  can  be  calculat* 
ed  in  a manner  similar  to  that  used  for  Equations  (III-9) 


through  (111-12) . It  then  follows  that 


f (e^-Y) 
f^(e,) 


f‘(e,-Y) 


f^(e2+Y) 


f^(90*-e,)  f‘^(9o»-e,) 


f‘(0,+Y)  Ilf  (90*-e,)  f‘(9o*-e,) 


f^(eo 


f‘(90*-e,)  f‘(90*-e,) 


f‘(e,) 


f‘(90*-e,)  f‘{90*-e,) 
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This  approximate  technique  is  a function  of  the  sensors  angle 
of  inclination  and  its  yaw  sensitivity  response.  It  is  not 
possible  in  general  to  state  that  the  approximate  technique 
will  lead  to  more  accurate  results  than  those  based  on  the 
cosine  law.  The  errors  introduced  should  be  checked  in  both 
cases  to  determine  which  method  is  more  applicable. 

Lastly,  there  may  be  times  when  mean  velocity  data  are 
to  be  corrected  for  deviations  from  the  cosine  law  after  the 
data  has  been  taken.  Since  all  one  wants  to  do  is  to  "adjust" 
the  experimental  values  for  the  deviations  from  the  cosine 
law,  one  would  not  necessarily  want  to  go  through  the  itera- 
tion scheme  previously  mentioned.  Utilizing  Equations 
(III-IO)  and  (III-ll) , the  corrections  can  easily  be  calcu- 
lated once  the  true  angle  of  the  mean  flow  with  respect  to 
the  probe  axis  is  known.  From  Equation  (IZI-9)  one  can 
write  for  an  ideal  X-probe 

tan  Yc  - f*(Y)  (III-20) 


Expanding  f*(Y)  in  a Taylor  series  around  one  can  also  writs 
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tan  Yc  - + i7pMY)|  O-Y^)  * 

L J YEV 
' 'c 

The  mean  flow  angle  can  then  be  computed  from 

tan  Y-  - f*(Y  ) 

Y = Y^  + 


(III-22) 


(III-23) 


Normally  the  first  order  correction  should  lead  to  accurate 
results.  For  sensors  which  show  large  deviations  from  the 
cosine  law  a second  order  correction  term  must  be  added  to 


Equation  (III-22) . Utilizing  Equations  (III-IO) , (III-ll)  and 
(III-23)  corrections  for  mean  velocity  measurements  can  be 
made  after  the  data  has  been  taken  if  the  yaw  sensitivity  of 


the  wire  is  known. 


Time-Dependent  Velocity  Field 

In  this  section,  expressions  for  the  turbulence 
velocity  correlations  will  be  developed  for  an  ideal 
X -probe  in  a low-intensity  turbulent  flow.  Figure  6 shows 
an  x-probe  array  aligned  with  the  x-axis  and  the  mean 
velocity  inclined  at  an  angle  y with  respect  to  this  axis. 
Here  the  notation  of  references  [17]  and  [22]  have  been 
used.  Using  the  assumption  of  low-intensity  turbulence, 
second  and  higher  order  terms  in  the  velocity  correlations 
are  neglected  and  the  instantaneous  velocity  can  be  written 
in  the  form 

q- 

Qi  - Qs  1 + ^ (111-24) 

For  sensor  1,  we  obtain  the  expression  for  the  instantaneous 
angle  of  incidence,  assuming  an  ideal  X-probe,  6^^  « ©2  “ 
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and  = K2 


where 


K,  to  be 
cos  6 


X/? 


* 5-  0; 


1 .1  - Wg 

X B COS  Y + sin  Y 


(III-25) 

(III-26) 


and  ^ = cos  y ~ sin  y 

Similarly,  it  can  be  shown  that  for  sensor  2 


cos  6j  - (1  - ^ g2 


The  effective  instantaneous  velocity  is  given  by 


(III-27) 


Qe,j  =Qif‘(ej) 


j»l,2 


(III-29) 


where 


f^(6j)  = [l-b.  <l-cos"^  6j)]^ 

Assuming  that  the  sensor  output  voltage  is  linear  in  the 
effective  cooling  velocity,  i.e.. 


one  can  write 


E . = KQ  . 
3 ej 


Ej  - K(Q^  + qgJf  (Sj) 


(III-30) 


Expanding  f (6^)  in  a power  series  and  keeping  only  first 
order  terms,  one  obtains 

* Kg(b^,m^,X)Qg  + Kq^g (b^  ,m^ , X)  + K^q^h (b^ ,m^ , X) 

(III-31) 

and 

E2  » Kg (b2,m2,C)Qg  + Kq^g (b2  ,m2 » C)  - KXq^h (b2 ,m2 , C) 


(III-32) 


where 


m 


_ 1 m 

g(bj^,mj^,X)  « Il+(2  ^ X ^-Db^^] 


(III-33a) 


g(b^,m2,K)  = [l+(2 


(lll-33b) 


1"2~  "’1“^ 

h(bj^,mj,A)  = 2 A ^ 


- ^ m 

b^m^Il+(2  ^ A ^-Dbj^] 

(III-33c) 


^2  ”'2 

1 "5“  in_  “1  — y—  in« 

h(b2,in2,C)  = 2 C b2in2[l+(2  C -l)b2l 

(III-33d) 


The  sensor  output  voltage  can  be  decomposed  into  a mean 


and  a fluctuating  component^  that  is 

Ej  * Ej  + Cj  ; 


j«l,2 


(III-34) 


With  this,  the  response  of  the  x-probe  to  the  velocity 
fluctuations  becomes 

* Kqgg(bj^,mj^,A)  + K5qj^h(bj^,mj^, A)  (III-35) 

®2  “ Kqgg(b2,m2fC)  - KAq^h(b2,m2,5)  (III-36) 

In  the  remainder  of  this  section  the  turbulence  velocity 
correlations  are  determined  in  two  different  coordinate 
systems:  (1)  one  whose  axes  are  along  and  normal  to  the 

mean  velocity  vector  and  (2)  one  which  is  along  and  normal 
to  the  X and  y axes  of  Figure  6. 

Streamwise  Coordinate  System.  From  Equations  (III-35) 
and  (111-36)  and  defining  ^ 

A = Ah(b2,m2,C)g(bj^,mj^,A)  + Ch(bj^  ,mj^ , A)  g (b2  ,m2 » C)  (III-37) 


t 


we  obtain  the  following  expressions  for  the  instantaneous 
values  of  the  fluctuating  streamwise  and  normal  velocity 
components: 
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I 


,-1.-1 


- K '■A  *IXh(b2,m2,C)ej^  + Ch(b^,in^,A)e2]  (III-38) 


■ k“^A”^  [g  (b2,m2»C)ej^  - g (bj^ ,mj^ , X)  62!  (III-39) 

Taking  the  time  average  of  Equations  (III-31)  and  (III-32) 
one  arrives  at  the  following  expressions  for  the  mean 
velocity 
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<3  <3 


2 / 1_  _ r \ 9 ^ ^o » 9 f ^ ) 


t1 


9 ( ^2  '**'2  ^ ^ ^ 


T ~? 


r 2 n 

fTl 

g (bj^,mj^,X) 

^^s 

1 — 2 

g (b2,m2,5). 

A^{h2,m2>^^ 

— 2 

L 9 (b2.m2C) 


9 (b2fm2C) 


ol 


(III-44) 


In  the  special  case  v^en  y=0,  bj^»b2  and  mj^«»m2f  these 
expressions  reduce  to  the  forms  presented  by  Ch2unpa9ne  [17] 
and  Friehe  and  Schvarz  122],  i,e., 

1 


2 

„2 


1 ‘®l+®2> 


(III-45) 


“I 

q 


1L=  1 

_ T[h(b,m)_ 


2 (ej^-e2) 


q q 


(b,m) 
mT 


9(b, 

.mb; 


(III-46) 


(III-47) 


When  the  x-probe  is  non-ideal,  i.e.,  when  each  sensor 
has  a sli9htly  different  yaw  response,  there  is  a correction 
factor  in  the  streamwise  velocity  correlation,  even  in  the 
case  when  y>0.  In  this  case,  the  Reynolds  stress  cannot 
be  determined  usin9  the  above  scheme  without  prior  deter- 
mination of  the  streaunwise  and  normal  velocity  correlations. 
When  the  mean  flow  angle  is  not  zero,  in  order  to  determine 


46 
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the  velocity  correlations  using  analog  techniques,  a 
weighted  sum  and  difference  may  be  used.  The  magnitude  of 
the  weighting  functions  depends  on  the  yaw  sensitivity 
coefficients  of  both  sensors  of  the  probe  as  well  as  the 
mean  flow  angle. 

x-y  Coordinate  System.  In  most  cases  where  the  mean 
velocity  is  at  an  angle  y with  respect  to  an  x-y  coordinate 
system  the  fluid  mechanics  experimenter  is  not  interested 
in  determining  the  local  streeunwise  and  normal  fluctuations 
but  rather  he  is  interested  in  the  u & v fluctuations  along 
the  X and  y axes  respectively.  The  streamwise  and  normal 
fluctuating  velocities  can  be  decomposed  into  components 
along  the  x and  y axes  by  the  following  transformation, 

= u cos  Y + V sin  y (III-48) 

q^  ■ -u  sin  y + v cos  y (III-49) 

Substitution  into  Equations  (III-35)  and  (III-36) 
yields  the  response  for  the  fluctuating  output  voltage. 

This  can  be  written  as 

®1  “ KS^,^(y)u  + KS^^(y)v  (III-50) 

and 

®2  ” KS2y(Y)u  - KS2^(y)v  (III-51) 

where  the  velocity  sensitivity  coefficients  are  given  by 

Siu(y)  • g (bj^,mj^,X)cos  y - 5h(bj^,mj^,X)  sin  y (III-52) 

Siv(y)  “ g(bj^,mj^,X)  sin  y + Ch(bj^,mj^,X)cos  y (III-53) 

S2u(y)  ■ g(b2»m2»C)cos  y + Xh(b2»n\2,5)  sin  y (III-54) 
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and 

S2v(y)  = -g (b2,m2C) sin  y + Xh(b2,m2, C) cos  y (III-55) 

Solving  Equations  (III-50)  and  (III-51)  for  the  instantaneous 
velocity  components,  we  obtain 

u - K"^A"^lS2^(Y)e^  + S^^(Y)e2]  (111-56) 

V - K"^A"^lS2y(Y)e^  - Sj^^(Y)e2]  (III-57) 

which  again  is  a simple  weighted  sum  or  difference  of  the 
output  voltages.  It  should  be  pointed  out  that  the  mean 
velocity  components  must  be  determined,  in  order  to  evaluate 
the  angle  of  incidence  y,  before  Equations  (III-56)  and 
(III-57)  can  be  used. 

In  the  case  of  an  ideal  x-probe  and  Y=0»  the  deviations 
from  the  cosine  law  in  the  computed  turbulence  velocity 

2 

correlations  vould  appear  in  the  correction  terms  (g/h) 
and  (g/h)  for  the  normal  and  bi -normal  turbulence  intensities 
as  given  by  Equations  (III-45  through  47) . A plot  of  various 
constant  g/h's  as  a function  of  b and  m is  shown  in  Figure  7. 
With  this  plot,  one  can  determine  the  difference  between 
these  calculated  values  assuming  either  the  cosine  law,  the 
Friehe  and  Schwarz  relation  or  the  present  relation. 

Discrepancies  Between  Different  Yaw  Relations  in 
Calculated  Values  of  Turbulence  Velocity  Correlations . In 
many  complex  flows  one  is  faced  with  cases  in  vdiich  yj^O . 

How  then  do  these  three  relations  compare  in  the  computed 
velocity  correlations? 
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Squaring  and  time  averaging  Equation  (III -56)  we 
obtain 

(u^)p  = K"2A"^]s2^(Y)e^  + SjL^(Y)e2]^  (III-58) 

If  one  assumes  the  cosine  law  for  yaw  sensitivity,  (i.e., 
b=l  and  m=0.5)  Equations  (111-50)  and  (III-51)  become 


*1  • 

ej  . K2‘'‘(u^  - v^) 


(III-59) 

(III-60) 


where  u^  and  v^  are  the  time  dependent  velocities  determined 
using  the  cosine  law.  Note  that  employing  the  cosine  law 
results  in  no  explicit  dependence  on  y*  Substitution  into 
Equation  (III-58) , with  the  aid  of  Equation  (III-20a) , 
yields 


e— = c,  + c»  P + c,  P 
u*  lu  2u  u 3u 


where 


uv 


=1U  - - 1 


-2 


and 


'2u 


'3u 


IS2v(y)  -S,„(y)] 


Iv' 


(III-61) 

(III-62a) 

(III-62b) 

(III-62C) 


The  ratios  and  P^^,  which  represent  the  distribution  of 
turbulence  energy  in  the  different  directions  and  their 
interaction,  are  defined  as 

p 


(III-63a) 


ti 
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(uv)^ 

^UV  ' — K- 

(U^) 


(III-63b) 


Equation  (III -61)  represents  the  difference  in  the  mean 
square  of  the  u component  velocity  fluctuations  as  would  be 
computed  by  the  present  relation  and  by  the  cosine  law. 

It  is  a function  of  the  coefficients  of  the  present  relation, 
the  angle  of  incidence ,y ,and  the  distribution  of  the 
turbulence  as  wDuld  be  computed  with  the  cosine  law. 

In,  many  cases  the  difference  in  the  rms  fluctuations  is 
of  more  interest.  It  can  be  calculated  by 


u* 


(III-64) 


where  u'  is  the  root  mean  square  of  the  u component 
fluctuations. 

Following  a similar  procedure,  the  difference  between 
the  mean  square  of  the  v component  fluctuations  computed 
using  the  present  relation  and  that  obtained  by  the  cosine 
law  becomes 


?=c  +^+c  - 

V Pu  P 


uv 

u 


(III-65) 


where 


'=2v  ■ ^ -®lu 


(III-66b) 


and 
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C3^  = a“^IS2^(y)  - Siy(Y)]^  (III-66C) 

For  the  rms  fluctuations  we  obtain 

e„,  - /e  +n:  - 1 (III-67) 

^ ~~T 

V* 

Next,  we  exeunine  the  difference  in  the  computed  values 
of  the  Reynolds  stress  uv  as  obtained  by  tvo  different 
methods.  The  first  method  is  based  upon  the  subtraction  of 
the  mean  square  of  the  linearized  fluctuating  output 
voltages  of  the  two  sensors.  For  this  case  one  obtains 


e “ c,  + c- 
— ^luv  2uv 


where 


(G;rr  +1)  , 

^UV  ^ ^uv 


-1 


=luv  ■ - 1 (III-69a) 


[S^„(Y)  -SJ„(Y)1 


'2uv 


2«=luv 


-1 


(III-69b) 


and 


_ IS;',tY)  -S^^(Y)I 


(III-69C) 


Computation  of  the  Reynolds  stress  in  this  manner  can 
introduce  large  discrepancies  since  it  is  sensitive  to  the 
differences  in  the  computation  of  the  other  velocity 
correlations. 

Using  analog  or  digital  techniques,  the  instantaneous 


values  of  u and  v can  be  determined  from  Equations  (III-56) 
and  (III -57) . Since  for  analog  signal  processing  squaring 
circuits  are  more  readily  available  than  multipliers,  we 
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determine  the  instantaneous  value  of  uv  from  the  relation 


I (u+v)  ^ 


(u-v) 


(III-70) 


and  then  time  average  the  signal . 

Using  this  approach  the  difference  in  the  values  of 
uv,  computed  by  the  present  relation  and  by  the  cosine  law, 
is  given  by 


e 

uv 


°luv 


u 


3uv  ^ 

Puv  Puv 


(III-71) 


where 


e 

uv 


(uv)p  - (uv)^ 


Equation  (III-70) 


The  coefficients  can  be  calculated  from 


-2 


and 

o*u„  - 2-l4-2[S2^(Y)S2„(l,)  Y S^^(Y)S^^(Y) 

■ S2v<Y)Sj^^(y)]  (III-72b) 

In  each  of  the  previous  cases,  the  differences  between 
the  Friehe  and  Schwarz  relation  and  the  cosine  law  can  be 
computed  from 


“x’f 


b.b, 
m=0 . 5 


This  is  true  since  when  the  present  relation  is  evaluated 
using  an  exponent  of  0.5  and  the  best  coefficient  corre- 
sponding to  it, the  yaw  dependence  becomes  that  proposed 


by  Friehe  and  Schwarz.  The  difference  then  between  the 
Friehe  and  Schwarz  relation  and  the  present  relation  in 
computing  a quantity  x is 
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(III-73; 

Using  the  previous  definitions  the  differences  in  computing 
the  velocity  correlations  can  be  determined  as  will  be 
presented  in  Chapter  VII. 

Response  of  Sensor  To  Velocity  Fluctuations 

In  all  of  the  previous  derivations,  we  have  assumed 
that  the  instantaneous  effective  velocity,  Q^,  was 
related  to  the  instantaneous  velocity,  Q^,  through  the  yaw 
relation,  i.e..  Equation  (III-29)  or 

°e  = 

where  g is  the  instantaneous  angle  between  and  the 
direction  normal  to  the  sensor.  Although  this  relation 
is  assumed,  only  the  yaw  dependence  of  the  mean  velocity 
is  measured  v^en  calibrating  a probe  for  its  yaw 
characteristics.  The  calibration  does  not  give  any 
information  about  the  yaw  dependence  of  the  time -dependent 
velocities,  e.g.,  the  rms  of  the  fluctuations. 

Since  Q^,  and  6 are  instantaneous  values,  they 
can  be  decomposed  into  time -mean  and  fluctuating  components 
and  written  in  the  form 

Qe  * ^e  (III -74) 

A 

g = o + a 


X - 


(III-75) 


where  a is  the  mean  yaw  angle  and  is  the  time  average 


of  Q^.  Assuming  low  intensity  turbulence,  the  relation 


0 . **  0 + Cf  + 

'^1  '^s  ^S  2Q 


(III-76) 


may  be  utilized.  Expanding  f (6)  using  Equation  (III-16) 


and  noting  that 


(III-77) 


one  obtains 


f^(6)  = f^(a)  + d[f^(o)]  % 
' dot 


+ p(a)  ‘^n  - d [f^(a))  '^s'^n  + O.^n. 


da 

s 


where 


2 2 2in  2 2 in 

p(a)  = b^m^  cos^  a tan^  a + 2f (a) [bm(m-l)  tan^a  - bm] cos  a 

2 7" 


(III-78) 


Applying  this  relation,  the  time  average  of  Equation  (III-29) 


becomes 


Q = Q f^(a)  + (p(a)  + f^(a) 1 ^ 
® ® ~~T~  Q 


(III-79) 


Under  normal  calibration  conditions  the  second  term  will 


always  be  negligible  compared  to  the  first  so  that  the 


calibration  for  mean  velocity  becomes 


Qe  = Qsf‘(ct) 


(III-80) 
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Squaring  and  time  averaging  Equation  (III-29) , keeping 
up  to  second  order  terms,  and  utilizing  Equation  (III-79) 
leads  to 


qf  • ♦ (a§(f^a)l)2  ^ 


- 2f2(a)  (^[f2(o)l)- 


*s  ‘n 
(III-81) 

This  equation  relates  the  contributions  of  the  streamwise, 
normal  and  bi-normal  turbulence  intensities  to  the  effective 
fluctuating  velocity,  i.e.,  the  time  dependent  cooling  to 
which  the  sensor  responds.  Two  special  cases  of  Equation 
(III-81)  are  significant.  The  general  case  arises  for  flows 
in  which  the  correlation  between  the  streamwise  and  normal 
velocity  fluctuations  can  not  be  ignored.  In  this  case,  the 
effective  rms  velocity  which  the  sensor  indicates  has  con- 
tributions from  the  streamwise,  normal  and  bi-normal  turbulence 
stresses  and  is  described  by  Equation  (111-81)  . 

For  the  special  cases  where  the  streamwise  and  normal 


velocity  fluctuations  are  uncorrelated,  = 0,  the  sensor 

is  only  sensitive  to  the  streamwise  and  normal  intensities. 
The  first  case  arises  when  periodic  streamwise  fluctuations 
are  superimposed  on  the  calibration  flow,  i.e.,  a very  low 


turbulence  intensity  flow.  With  this  stipulation,  = 0 

“7  ”7 

and  *0  that  Equation  (III-81)  reduces  to 


q;  • q;f^a) 


(III-82) 


and  4^  arc  r»B  values.  Combining  this  with  Equation 


• • • . • 
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^ « 0^  = constant  (III-83) 

This  equality  states  that  the  effective  local  turbulence 
intensity,  measured  at  any  yaw  angle,  remains  constant. 

This  result  will  be  experimentally  examined  in  Chapter  V. 

The  last  special  case  of  Equation  (III-81)  also 
results  in  contributibns  to  the  effective  rms  velocity  from 
fluctuations  other  than  those  in  the  streamwise  direction. 
When  such  contributions  are  not  negligible  Equation  (III-82) 
will  not  hold,  i.e.. 


‘^e  2 

A t f (a)  (III-84) 

^s 

Contrasting  Equations  (III-80)  and  (III-84)  one  observes 
that  yaw  relations  of  hot-wires  and  hot-films  obtained 
from  calibration  in  a steady  velocity  field  may  only  be 
applied  to  mean  and  instantaneous  velocities  measured  in 
a turbulent  flow.  They  can  not  be  used  for  the  rms  of  the 
velocity  fluctuations.  A discussion  of  the  implications 
of  the  above  results  on  the  present  experimental  work,  as 
well  as  on  previous  ones,  e.g.,  Ilorrison  et  al^.  [26]  and 
Bruun  [27],  is  presented  in  Chapter  VI. 
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CHAPTER  IV 

EXPERIMENTAL  FACILITIES  AND  TECHNIQUES 
Experimental  Facilities 

Two  main  facilities  were  utilized  during  the  course 
of  the  experimental  part  of  the  present  study:  the  I.I.T. 
Environmental  Wind  Tunnel  and  the  I.I.T.  Calibration 
Tunnel. 

Environmental  Wind  Tunnel . The  I.I.T.  Environmental 
Wind  Tunnel  operates  in  a closed  return  mode  thus  permit- 
ting use  of  two  test  sections.  The  present  experiments 
utilized  the  high-speed  test  section.  The  dimensions  of 
this  section  are  2 ft.  deep,  3 ft.  high  and  10  ft.  long, 
and  the  free  stream  velocity  can  be  controlled  at  any 
speed  up  to  65  ft/sec.  Through  the  use  of  turbulence 
manipulators  (281,  a turbulence  level  much  less  than  0.1% 
is  maintained  throughout  the  test  section.  A schematic 
showing  the  side  view  of  the  tunnel  from  the  high-speed 
side  is  shown  in  Figure  8.  For  a more  detailed  descrip- 
tion of  the  wind  tunnel  the  reader  is  referred  to  the  re- 
port by  Tan-atichat  and  Nagib  [29]. 

Calibration  Tunnel . The  I.I.T.  Calibration  Tunnel 
was  utilized  in  two  modes  during  the  course  of  the  inves- 
tigation. In  the  experiments  of  Chapter  V it  was  used  in 
its  standard  configuration  for  velocity  calibration  of  the 
hot-wire  probes.  It  was  also  used  in  Chapter  VI  in 
another  configuration  to  produce  a flowfield  with  periodic 
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velocity  fluctuations  superimposed  on  the  mean  streamwise 
velocity.  In  both  tunnel  configurations,  compressed  air 
enters  an  acoustically  baffled  plenum  chamber  and  flows 
through  a bellmouth  into  a 3 in.  diameter  plexiglass  duct. 

In  the  calibration  mode,  a nozzle  is  attached  to  the 
end  of  the  duct  to  produce  a leuninar  free  jet  flowfield. 

The  probes  are  placed  in  the  core  of  the  one-inch  dieuneter 
jet  as  shown  in  the  bottom  photograph  of  Figure  9.  De- 
tailed information  regarding  the  construction  and  calibra- 
tion of  the  tunnel , and  on  the  turbulence  manipulators  used 
in  this  configuration  are- presented  by  Loehrke  and 
Nagib  [28] . 

A schematic  of  the  tunnel  in  the  second  flow  configu- 
ration is  shown  in  Figure  10.  To  produce  the  desired  flow 
conditions, a screen,  two  perforated  plates  and  two  honey- 
combs are  used  as  turbulence  manipulators  in  two  straight 
sections  of  the  duct.  The  straight  ducts  are  connected 
with  a plexiglass  "T"  section  as  shown  in  the  top  photo- 
graph of  Figure  11.  One  end  of  the  "T"  section  is  con- 
nected to  a conical  contraction  from  a 15  in.  speaker  which 
produces  the  periodic  forcing  of  the  flow. 

The  other  end  of  the  duct  is  open  to  the  laboratory.  The 
first  utilization  of  the  tunnel  in  this  mode  is  described 
by  Marcichow  et  al.  [30] . 

Angularity  Chuck . A special  angularity  chuck  de- 
signed and  built  by  Roberts  [23l  was  utilized  throughout 
the  course  of  this  investigation.  The  chuck  can  accurately 
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position  the  probes  at  various  angles  of  yaw  to  within 
+0.04".  The  chuck  is  shown  attached  to  the  side  of  the 
Environmental  Wind  Tunnel  in  the  top  photograph  of  Fig- 
ure 9,  and  positioned  in  the  traversing  mechanism  next  to 
the  Calibration  Tunnel  in  Figure  11.  In  conjunction  with 
the  angularity  chuck,  a specially  designed  probe  support 
system  was  used  to  attach  the  probes.  It  allowed  the 
velocity  probes  to  be  rotated  either  about  the  midpoint  of 
the  sensor  or  about  any  point  along  the  length  of  the 
probe  body.  In  acquiring  mean  yaw  calibration  data,  the 
probe  was  rotated  about  the  midpoint  of  the  probe  body 
(as  shown  in  Figure  9) , while  for  the  unsteady  yaw  cali- 
bration it  was  rotated  about  the  midpoint  of  the  sensor. 
The  arrangement  of  the  support  system  for  the  latter  case 
is  visible  in  the  photograph  of  Figure  11. 

Velocity  Sensors.  The  single-  and  X-wlre  probes  used 
were  designed  based  on  the  work  of  Comte-Bellot  ^ al. [31] 
and  Strohl  et  a]^.  [32]  to  minimize  the  aerodynamic  inter- 
ference on  the  sensors.  In  each  case,  the  probe  body  was 
constructed  from  0.25  inch  O.D.  stainless  steel  tubing 
18  inches  in  length.  The  prongs  of  the  single-wire  probes 
are  1 inch  long, 0.016  inch  in  mean  diameter  and  the 
spacing  between  their  tips  is  equal  to  0.20  in.  The 
prongs  of  the  X-wire  are  0.75  in.  long,  and  0.016  in.  in 
mean  diameter  with  a spacing  of  0.12  in.  between  prongs  at 
the  point  where  the  wires  are  soldered.  The  separation 
distance  between  the  wires  is  0.04  in. 
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The  sensors  used  in  all  measurements  are  made  of 
0.00015  in.  dieuneter  tungsten  wire.  The  sensor  l/d  was 
varied  from  75  to  1000  for  single  sensor  probes / while  a 
standard-length  sensor  (^/d=580)  was  used  for  all  X-wire 
measurements.  The  active,  i.e.,  heated,  length  of  the 
wire  was  controlled  by  changing  the  length  of  the  segments 
of  the  wire  coated  with  copper.  The  wires  are  soldered 
to  the  prongs  with  the  active  length  centered  between 
them.  Photographs  of  some  of  the  probes  are  shown  in 
Figure  12. 

Angularity  Measurements 

Steady  Velocity  Calibration.  Mean-velocity  yaw  cali- 
bration for  single-  and  X-wire  probes  was  performed  in  the 
Environmental  Wind  Tunnel.  However,  the  probes  were  first 
connected  to  a DISA  55D01  constant  temperature  anemometer 
and  velocity  calibration  data  were  obtained  by  placing 
them  in  the  jet  of  the  Calibration  Tunnel  positioned  next 
to  the  Environmental  Wind  Tunnel;  see  Figure  9.  Using  the 
angularity  jig  shown  in  the  bottom  photograph  of  Figure  9, 
X-probes  were  rotated  to  orient  each  of  the  wires  normal 
to  the  flow  direction  during  its  calibration.  Next,  the 
probes  were  connected  to  the  angularity  chuck  attached  to 
the  side  of  the  wind  tunnel.  The  probes  were  supported  at 
the  midpoint  of  the  body  and  aligned  such  that  the  prongs 
remained  in  the  same  vertical  plane  throughout  as  much  of 
the  rotation  as  possible.  The  probe  was  then  connected  to 
the  seme  DISA  55D01  constant  temperature  anemometer  using 
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the  same  overheat  ratio  utilized  in  the  velocity  cali- 
bration. 

At  any  constant  tunnel  velocity,  the  ambient  tempera- 
ture remained  constant  to  within  +0.5®F.  However,  the 
tunnel  air  temperature  varied  with  tunnel  velocity  so  that 
the  output  was  manually  adjusted  at  each  free  stream  velo- 
city to  yield  the  output  voltage  calculated  from  the  cali- 
bration curve.  The  adjustment  was  provided  by  the  variable 
bridge  resistance  in  the  anemometer.  The  accuracy  of  this 
scheme  of  compensation  of  the  anemometer  output  for  varia- 
tion in  the  tunnel  air  temperature  is  verified  in  Chap- 
ter II  and  Appendix  B.  The  anemometer  output  was  then  dis- 
played on  a Tektronix  564  oscilloscope,  time  averaged  for 
10  sec.  with  a Heath  multimeter  and  also  recorded  on  an 
H.P.  Mosley  strip  chart  recorder.  A photograph  and  a sche- 
matic of  the  instrumentation  for  these  measurements  are 
shown  in  Figures  9 and  13  respectively. 

To  measure  the  mean  tunnel  velocity  a Validyne  DP-45 
pressure  transducer  connected  to  a pitot  static  probe  was 
employed.  This  transducer  was  powered  by  a CD-15  carrier 
demodulator  and  its  gain  set  to  obtain  the  value  calibrated 
by  Corke  et  al.  [33]  of  405.28  volts/psld.  The  voltage 
proportional  to  the  differential  pressure  was  connected  to 
the  analog  averaging  circuit  of  Figure  13,  which  has  a 
time  constant  of  4 seconds.  The  averaged  output  was  then 
monitored  on  a Dana  Model  3300A  digital  multimeter  and 
displayed  on  Channel  2 of  the  oscilloscope.  For  further 
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details  of  the  pressure  transducer  system  and  its  cali- 
bration, the  reader  is  referred  to  Corke  et  [ 33]  . 

Fluctuating  Velocity  Calibration.  For  this  portion 
of  the  investigation,  the  modified  calibration  tunnel  was 
used  and  the  velocity  fluctuations  were  produced  via  a 
15  inch,  25  watt  Altec  Biflex  speaker.  The  probe  was 
attached  to  the  angularity  chuck  using  special  attachments 
so  that  the  probe  rotated  about  the  midpoint  of  the  sensor 
and  the  entire  sensor  remained  in  the  test  volume  shown 
in  Figure  10.  The  angularity  chuck  was  attached  to  a 
milling  machine  indexing  head  so  that  four  degrees  of  free- 
dom were  available  (three  translational  and  one  rotational) . 
Both  single  and  X-probe  measurements  were  made  and  the 
schematics  of  the  instrumentation  employed  are  shown  in 
Figures  14  and  15.  The  bottom  photograph  of  Figure  11  also 
depicts  the  setup  of  the  instrumentation.  In  either  case, 
the  output  of  a sensor  was  fed  to  a DISA  55D01  Constant 
Temperature  Anemometer,  utilizing  an  overheat  ratio  of  1.8, 
and  then  linearized  by  a DISA  55D10  Linearlzer.  This  signal 
which  is  proportional  to  the  velocity,  was  time 
averaged  with  a Heath  Digital  Instrument  model  EU805  for 
10  seconds.  The  accuracy  of  the  linearized  velocity,  i.e., 
+1%,  was  comparable  but  not  as  good  as  that  achieved  with 
the  non-linearized  operation  discussed  in  the  next  sec- 
tion. To  obtain  the  rms  of  the  signal,  the  linearized  out- 
put was  fed  into  a DISA  RMS  meter,  Model  55D35,  and  its 
output  read  using  an  H.P.  3440A  Digital  Voltiseter.  The 
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output  of  a Briiel  and  Kjaer  1022  Beat  Frequency  Oscillator 
powered  the  speaker  and  served  as  a reference  signal  to  the 
PAR  Waveform  Eductor  and  to  the  Acton  Labs  329B  Phasemeter. 
To  obtain  the  rms  of  the  periodic  fluctuations,  the  line- 
arized signal  was  connected  to  the  Eductor.  The  Waveform 
Eductor  serves  as  a periodic  averaging  device  triggered  by 
the  signal  producing  the  velocity  fluctuations.  The 
educted  output,  which  only  contained  the  part  of  the  sig- 
nal correlated  to  the  periodic  fluctuation,  was  then  fed 
to  the  RMS  meter  to  determine  its  level.  At  each  forcing 
frequency  exeunined,  the  phasemeter  was  used  to  measure  the 
phase  differences  between  the  reference  signal  and  the 
educted  sensor  signal.  The  phasemeter  output  was  monitor- 
ed while  traversing  the  probe  across  the  test  section  to 
insure  that  the  periodic  forcing  was  planar  over  the  test 
volume. 

Data  Processing 

Data  from  the  experiments  were  reduced  utilizing  a 
UNIVAC  1108  computer  and  a digital  data  acquisition  and 
processing  system  [34]  . Velocity  calibration  data,  with 
the  sensors  normal  to  the  flow,  were  fit  to  the  following 
equation: 

E^  « A*  + B*  u"  (lV-1) 


The  values  of  A',  B*  and  n were  obtained  by  minimizing 
the  standard  deviation  of  the  data  from  the  associated 
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fit.  The  equation  was  fit  for  a range  of  velocities  from 
10  to  50  ft/sec.  The  normalized  standard  deviation  of  the 
calculated  velocity  from  the  data  for  all  sensors  was 
approximately  0.4%.  This  normalized  standard  deviation  is 
defined  as 


^calculated 

^measured 


At  the  four  free-stream  velocities  of  12.5,  25,  37.5  and 
50  ft/sec  , calibration  data  were  obtained  for  the  yaw 
angles  -70*<a<70*.  At  each  position  for  the  yaw  angles 
the  anemometer  output  voltage  and  the  angle  of  yaw  with 
respect  to  the  flow  direction  were  recorded.  The  flow  di- 
rection was  initially  determined  by  maximizing  the  anemo- 
meter output.  The  wire  output  voltage  was  then  plotted 
versus  yaw  angle  for  a values  in  the  range  +30“  to  deter- 
mine the  angle  at  which  the  mean  flow  was  normal  to  the 
sensor,  i.e.,  0.  In  order  to  perform  accurate  compari- 

sons between  the  different  yaw  relations,  the  angle  was 
redetermined  to  a higher  degree  of  precision  in  the  data 
reduction  routine  as  will  be  discussed  in  Chapter  V.  The 
yaw  data  were  obtained  at  increments  of  three  degrees  in 
all  cases  and  every  one  degree  in  some  cases. 

Wire  output  voltages  were  then  converted  to  mean  velo- 
cities via  Equation  (IV-1)  and  normalized  by  the  velocity 
calculated  when  a * a^,  i.e.,  U^.  The  data  was  then  fit 
to  the  following  relation 
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m 

^ ■ a + b cos  a 

^00 


(IV-2) 


to  obtain  the  values  of  the  coefficients  a,  b and  m which 
minimize  the  standard  deviation  of  the  fit.  These  values 
are  denoted  by  b*  and  m*.  The  coefficient  a can  be 
replaced  in  Equation  (IV-2)  by  1-b,  since  a+b  must  equal 
to  one  to  satisfy  the  conditions  when  a » 0.  Therefore, 
we  can  rewrite  Equation  (IV-2)  as  a two  pareuneter  yaw  relation. 


^ = [ 1 - b(l.-  cos^a)]^  - f2(a)  (iv-3) 

For  the  rest  of  this  report  this  equation  will  be  referred 
to  as  the  present  yaw  relation. 

A grid  of  b and  m values  was  then  specified  and  at 
each  grid  point,  the  standard  deviation  of  U /U  from  the 
present  relation  was  determined.  The  standard  deviation, 
o^,  is  defined  as 

) - ^e,i  1 ^ (IV-4) 

J 

where  N is  the  number  of  data  points.  When  the  value  of 

the  exponent  m * 0.5,  then  Op»Oj,  where  is  the  standard 

deviation  of  the  data  from  the  Friehe  and  Schwarz  relation 

[22].  For  the  case  when  the  coefficient  b«l  and  m«0.5, 

a *0* , where  o*  is  the  standard  deviation  of  the  fit  using 
pc  c 

the  cosine  law. 


/i  ? r 

» ilA 
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These  values  were  then  written  onto  magnetic  tape  and 
processed  on  the  PDP  11-10  utilizing  a contour  plotting 
routine  in  the  Data  Acquisition  and  Processing  System 
Library.  This  contour  plotting  routine  was  written  by 
Richard  Wlezien  during  the  course  of  the  present  investi- 
gation. 
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[ CHAPTER  V 

I 

YAW  CALIBRATION  FOR  STEADY  VELOCITIES 

This  chapter  examines  the  influence  of  various  para- 
meters on  the  yaw  calibration  of  single  sensors  and  X-wire 
probes  utilizing  three  different  relations  for  the  yaw 
sensitivity;  the  cosine  law,  the  Friehe  and  Schwarz  rela- 
tion and  the  present  relation;  see  Chapter  I for  details. 

As  discussed  in  the  previous  chapter,  the  mean  yaw 
calibration  data  was  initially  centered  graphically  to 
determine  the  mean  flow  direction  which  corresponds  to 
0^=0.  An  example  of  this  is  shown  in  Figure  16,  where  the 
anemometer  output  voltage  is  plotted  versus  angle  of  yaw. 
The^  distance  between  the  curve  and  the  initial  line  of 
symmetry  was  measured  at  different  heights  and  used  to 
calculate  the  average  shift  in  the  center  line.  The  val- 
ues measured  in  the  example  of  Figure  16  are  shown  on  the 
graph.  This  process  was  carried  out  for  every  run  to  en- 
sure that  there  was  no  slippage  of  the  probe  in  the  angu- 
larity chuck.  Such  a procedure  is  sufficient  for  angular- 
ity calibration  for  normal  experimental  work.  However, 
for  the  purpose  of  doing  an  accurate  comparison  of  the 
standard  deviations  determined  from  the  three  yaw  rela- 
tions as  influenced  by  different  parameters  this  was  not 
sufficient.  The  reason  behind  this  is  found  in  Figure  17. 

Using  the  initially  shifted  data,  the  offset  angle, 
a^,  was  varied  from  -9*  to  't’9*  in  equal  increments  and 
the  siinimusi  standard  deviation,  o*.  corresponding  to  each 
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value  of  Oq  was  found  using  the  present  relation.  The  de- 
viation of  symmetric  data  from  an  even  yaw  relation,  as 
calculated  in  this  approach,  should  be  symmetric  about  the 
true  mean  flow  direction.  As  demonstrated  graphically  in 
Figure  17,  the  true  line  of  symmetry  appears  to  be  -0.18® 
away  from  that  determined  from  centering  the  data  by  the 
procedure  of  Figure  16.  Using  this  corrected  value,  the 
typical  variation  of  o*  with  offset  angle  (i.e.  probe  mis- 
alignment) is  shown  in  Figure  18.  At  a =0,  0*  achieves 

op 

itr,  minimum  value  of  0.15%  and  increases  linearly  with  off- 
set angle.  In  all  cases  examined,  the  variation  of  a* 
with  over  this  range  was  independent  of  sensor  t/d  or 
Re^  and  equal  to  0.8%  per  degree  of  misalignment.  The 
data  processing  program  included  a subroutine  which  auto- 
matically performed  this  accurate  alignment  of  the  data 
and  assigned  the  corrected  yaw  angles  to  all  data  points. 

Using  this  technique  to  determine  the  true  mean  flow 
angle,  the  variations  of  the  minimum  standard  deviations 
determined  using  the  cosine  law,  Friehe  and  Schwarz's  and 
the  present  relation  were  examined.  A typical  result  is 
shown  in  Figure  19  for  a sensor  t/d  of  250  and  a free 
stream  velocity  of  37.5  ft/sec.  If  the  probe  was  properly 
aligned,  an  error  of  1.8%  in  the  normalized  velocity, 
U^/U^^%fOuld  be  the  minimum  achievable  with  the  cosine  law. 
This  value,  which  indicates  the  goodness  of  the  fit.  was 
typically  equal  to  0.%%  if  the  Friehe  and  tchwari  relation 
wee  uee«)  eee  redu<'e<)  eci'T  <'«  iMtel>  0 1%  by  uoine 
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! the  present  relation.  The  values  of  a*  and  a|  increase 

by  0.8%  per  degree  of  misalignment  from  the  accurately  de- 
termined Oq.  If  the  probe  is  off  by  more  than  two  degrees 
o*  also  varies  linearly  by  0.85%  per  degree.  However, 

I a*  does  not  vary  linearly  with  for  values  of  offset 

i 

I less  than  two  degrees.  In  the  case  of  Figure  19,  a factor 

; of  four  in  accuracy,  as  measured  by  the  goodness  of  the 

!l 

ji  fit  of  the  relation  (i.e.,  o*)  , is  gained  by  using  the 

;j  Friehe  and  Schwarz  relation  over  the  cosine  law,  and  an- 

other factor  of  five  is  gained  in  going  from  the  Friehe 
and  Schwarz  relation  to  the  present  relation. 

I 

The  procedures  described  above  and  in  the  last  part 
of  the  previous  chapter  are  employed  for  single-wire 

! probes  in  the  following  sections  and  with  X-probes  later 

i in  this  chapter. 

Single-Wire  Probes 

The  best  fit  values  of  the  coefficients  which  mini- 
i mize  the  standard  deviation  for  the  present  relation  are 

! denoted  by  b*  and  m*.  For  the  Friehe  and  Schwarz  rela- 

I'  tion  ht  corresponds  to  the  minimum  value  of  standard  de- 

I viation,  oj.  One  method  of  determining  these  coefficients 

is  shown  in  Figure  20  where  the  deviations  away  from  the 
cosine  law  are  plotted  for  a standard-length  sensor 
I (f/d«580)  as  a function  of  Reynolds  number.  Since  the  de- 

via t ions  are  linear  in  cos  a,  the  value  of  b*  can  be 
I computed  from  the  slope  of  the  straight  line.  This 
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approach  is  similar  to  the  one  used  by  Friehe  and  Schwarz 
[22].  Both  positive  and  negative  deviations  from  the 
cosine  law  are  found  in  the  data  of  Figure  20.  This  be- 
havior will  have  important  implications  on  the  various 
yaw  relations  and  is  discussed  in  Chapter  VII . 

Another  method  to  indicate  the  regions  where  the  co- 
efficient b and  exponent  m minimize  the  standard  devi- 
ation is  based  on  examining  contours  of  constant  o in 
P p 

the  b-m  plane.  Using  this  technique  a comparison  be- 
tween all  three  yaw  relations  can  be  performed  at  the  same 
time,  since  when  m=0.5  the  present  yaw  relation  reduces 
to  the  Friehe  and  Schwarz  relation  and  when  b=l  and  m=0.5 
the  present  relation  reduces  to  the  cosine  lav , Typical 
examples  of  this  approach  are  shown  in  Figures  21  and  22. 
For  the  case  shown  in  Figure  21  {f/d=250)  , the  p>osition 

of  the  region  of  minimum  standard  deviation,  i.e.  less 
than  0.5%,  is  quite  removed  from  either  the  cosine  law  or 
the  Friehe  and  Schwarz  relation  as  depicted  graphically 
on  the  figure.  When  a longer  sensor  is  used,  as  in  Fig- 
ure 22,  this  region  will  lie  closer  to  the  cosine  law  and 
pass  through  the  Friehe  and  Schwarz  relation, although  the 
point  of  minimum  standard  deviation  is  still  removed  from 
these  yaw  relations. 

The  variations  of  these  regions  of  minimum  standard 
deviation  with  Reynolds  number  is  shown  in  Figure  23  for  a 
sensor  with  -t/d«250.  As  the  Reynolds  number  increases, 
these  regions  tend  to  move  upward  and  to  the  left  of  the 
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f figure,  although  remaining  far  removed  from  both  the  co- 

! 

sine  law  and  the  Friehe  and  Schwarz  relation  at  the  high- 
est velocity  used,  i.e.  50  ft/sec. 

Since  in  many  of  the  cases  examined  the  position  of 
the  minimum  standard  deviation  was  far  removed  from  either 
of  the  other  two  yaw  relations,  we  were  concerned  if  this 

i* 

type  of  behavior  was  due  to  the  manner  of  fabrication 
of  the  probes.  To  this  end,  an  independent 

[ set  of  yaw  calibration  data  was  obtained  from 

il 

Or.  John  Foss  of  Michigan  State  University  and  processed 
in  the  manner  described.  S^unple  results  for  a sensor 
f/d“200  are  shown  in  Figure  24  while  a more  complete  de- 
scription is  given  in  Appendix  D.  From  this  and  the  fig- 
^ ures  in  Appendix  D one  can  observe  that  similar  trends  are 

exhibited.  In  particular,  at  low  values  of  free-stream 
velocity,  the  region  of  minimum  standard  deviation  is  far 
removed  from  the  cosine  law  and  the  Friehe  and  Schwarz  re- 
lation. As  the  free-stream  velocity  increases,  these 
regions  shift  upward  and  to  the  left  and  in  the  case  of 
Figure  24,  they  display  regions  of  overlap.  Results  for 
a standard- length  wire  are  also  shown  in  the  same  figure. 
This  data  exhibits  more  similarity  to  the  MSU  results 

(shown  above  it  on  the  same  figure)  than  the  shorter 

il 

length  sensor  contours  of  Figure  23. 
j!  Effect  of  Probe  Misalignment.  We  next  examine  the 

variations  of  the  minimum  standard  deviations  obtained 
from  using  the  three  yaw  relations,  as  well  as  the  best 
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fit  coefficients  in  the  Friehe  and  Schwarz  and  present  re- 
lations/ with  offset  angle  (i.e.,  probe  misalignment)  and 
Reynolds  number.  Results  for  a sensor  with  an  f/d*>250  are 
presented  first  and  are  shown  in  Figures  25  through  29 
as  contours  cf  the  parameter  or  variable  examined  in  the 
Reynolds  number  versus  angle  of  misalignment  plane.  Fig- 
ure 25  indicates  that  o*  varies  by  approximately  0.85% 
per  degree  of  offset/  as  found  before/  and  that  the  mag- 
nitude of  a*  decreases  as  the  Reynolds  number  increases. 

The  variation  of  o|  is  shown  in  Figure  26  where  again  it 
varies  linearly  by  0.8%  per  degree  of  offset  and  decreases/ 
although  at  a smaller  rate  than  in  the  cosine  case/  as  the 
Reynolds  number  increases.  The  minimum  standard  deviation 
found  using  the  present  relation  is  shown  in  Figure  27. 

In  this  case/  o*  varies  linearly  by  0.8%  per  degree  of 
offset  as  found  from  Figures  18  and  19/  but  unlike  the 
previous  two  yaw  relations,  a*  is  independent  of  Reynolds 

number.  The  reason  that  o*  and  at  are  functions  of 

c z 

Reynolds  numbers  is  that  the  region  of  minimum  standard 
deviation  in  the  b-m  plane  shifts  toward  these  rela- 
tions as  the  Reynolds  nvimber  increases,  thereby  reducing 
the  difference  between  relations.  It  is  of  course  more  de- 
sirable to  have  the  accuracy  of  representing  the  calibra- 
tion independent  of  velocity,  as  achieved  by  the  present 
relation. 

Variations  in  bj  and  b*  with  and  Re^  are  shown  in 
Figures  28  and  29  respectively.  From  Figure  28  one 
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observes  that  b|  increases  with  Reynolds  number  and  tends 
toward  the  cosine  law.  It  is  important  to  note  that  small 
differences  between  the  mean  flow  direction  determined 
graphically  and  by  the  accurate  method  depicted  in  Fig- 
ure 17  will  have  a negligible  effect  on  the  computed  value 
of  bj.  Similar  trends  are  observed  for  b*  in  Figure  29, 
where  again  its  value  increases  with  Reynolds  number  and 
a small  (less  than  0.2")  change  in  offset  angle  will  usu- 
ally have  a small  effect  on  the  value  of  b*.  In  general, 
if  the  probe  was  misaligned  by  0.5",  the  effect  of  the 
best  fit  coefficients  is  equivalent  to  that  resulting  from 
calibration  at  a free  stream  velocity  3 ft/sec  less  than 
the  desired  one. 

The  variation  of  these  coefficients  with  Reynolds 
number  and  offset  angle  for  a standard-length  sensor  are 
shown  using  similar  plots  in  Figures  30  through  35.  The 
behavior  of  a|  and  a*  are  similar  to  that  of  the  previous 
case  and  is  shown  in  Figures  30  and  31 , respectively,  al- 
though for  this  standard  length  probe,  o|  depends  less  on 
Reynolds  number.  From  Figure  32,  we  observe  that  b|  has  a 
value  larger  than  1 at  low  Reynolds  numbers  and  decreases 
with  increasing  Reynolds  number.  This  behavior  is  oppo- 
site that  of  the  sensor  with  an  i/d=250.  However, the 
variations  of  b|  with  offset  angle  are  the  same  as  in  the 
previous  case.  Figures  33  and  34  examine  the  behavior  of 
the  best  fit  parameters  determined  from  the  present  rela- 
tion. We  note  that  b*  increases  with  Reynolds  number 
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while  m*  decreases r once  again  indicating  the  shifting  of 
the  region  of  minimum  standard  deviation  in  the  b-m 
plane.  Their  dependence  on  offset  angle  is  identical  to 
that  for  the  sensor  with  an  f/d>250. 

Associated  with  the  values  of  b*  and  m*  is  the  tur- 
bulence correction  factor  (g/h)*  described  in  Chapter  III. 
From  Figure  35,  (g/h) * increases  by  5%  over  the  indicated 
Reynolds  number  range  (i.e.,  values  of  free-stream  velo- 
city from  12.5  to  50  ft/sec).  In  this  case,  if  the  probe 
is  misaligned  by  0.5®,  the  value  of  (g/h)*  determined 
would  be  equivalent  to  recalibrating  the  probe,  properly 
centered,  in  a free-stream  with  a velocity  3 ft/sec  higher 
than  the  desired  value. 

Effect  of  Reynolds  Number . With  the  probe  properly 
centered  the  effects  of  changing  the  free-stream  velocity 
on  the  various  minimum  standard  deviations  and  best  fit 
coefficients  are  examined  next. 

If  the  values  of  o*  and  at  are  normalized  by  o*,  we 

c r P 

obtain  ratios  of  improvement  in  going  from  the  cosine  law 

or  the  Frlehe  and  Schwarz  relation  to  the  present  relation, 

respectively.  The  behavior  of  these  ratios  can  be  seen 

in  Figure  36  for  two  different  values  of  sensor  t/d.  For 

both  cases  the  ratio  oj/o*  decreases  with  increasing  Re^. 

The  magnitude,  however,  is  a function  of  sensor  t/d. 

With  the  long  sensor  (f/d«1000) , o*/o*  decreases  from  a 

c p 

value  of  25  at  a Reynolds  number  of  0.9  to  a value  of  3 at 
a Re^»3.6.  For  the  short  sensor  (t/d^lS) , over  the  same 


Re^  range,  the  ratio  decreases  from  9 to  3.  Surprisingly, 
this  indicates  that  in  some  cases  the  shorter  sensor  fits 
the  cosine  law  better  than  the  longer  sensor. 

Two  different  trends  are  observed  for  the  ratio 
o|/o*.  For  the  long  sensor  the  ratio  increases  from  1 at 
Re^“0.9  to  1.75  at  Re^*3.6.  This  is  opposite  to  the  trend 
observed  for  the  short  sensor  where  the  ratio  of  improve- 
ment decreases  from  a value  of  3.5  at  Re^=0.9  to  unity  at 
Re^*=3 . 6 . 

For  a sensor  with  f/d=1000,  we  next  examine  the 
Reynolds  number  dependency  of  the  best  fit  parameters  of 
the  yaw  relations.  These  are  shown  in  Figure  37.  As  the 
Reynolds  number,  Re^,  increases,  b|  decreases  by  10%, 
b*  decreases  by  15%  and  m*  decreases  by  10%.  In  this  par- 
ticular case,  the  yaw  sensitivity  coefficients  are  strong 
functions  of  the  Reynolds  number.  This  behavior  is  ex- 
pected and  will  be  discussed  in  Chapter  VII.  This  depend- 
ence can  also  lead  to  many  problems  when  using  analog 
techniques  to  measure  turbulence  velocity  correlations. 

Effect  of  Sensor  f/d.  By  changing  the  sensor  ■fc/d, 
previous  investigators  have  found  that  the  axial  tempera- 
ture distribution  along  the  wire  is  significantly  altered. 
This  may  lead  to  large  deviations  from  the  cosine  law. 
Therefore,  the  effect  of  sensor  t/d  on  the  minimum  stan- 
dard deviation  and  on  the  parameters  in  the  three  yaw  re- 
lations was  investigated. 


Five  different  values  of  f/d  were  tested:  75,  145, 
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75 
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250 f 580  and  1000  while  keeping  the  prong  spacing  con- 
stant. Two  identical  probes  were  used  in  collecting  the 
data,  and  one  of  them  (Probe  1)  was  used  with  the  f/d's 
of  75,  580  and  1000.  The  effect  of  sensor  f/d  on  t^e 
hot-wire  velocity  calibration  constants  of  Equation  (IV-1) 
is  shown  in  Figure  38.  The  coefficient  B'  remains  fairly 
constant  with  sensor  t/d  while  A'  is  a strong  function  of 
l/d.  The  exponent  n decreases  with  decreasing  values  of 
f/d.  The  exponent  of  King's  law,  n,  seems  to  tend  toward 
the  theoretical  value  of  one  half  as  £/d  approaches  very 
large  vlaues.  The  velocity  calibration  constants  obtained 
from  the  standard  length  sensors  of  the  X-wire  are  also 


I, 

shown.  Recalling  that  the  prong  spacing  was  different 
between  the  single-wire  and  the  X-wire  probes,  one  concludes 
that  probe  geometry,  including  prong  shape  and  spacing, 

I 

j greatly  affects  the  values  of  B'  and  n. 

To  separate  out  any  Reynolds  number  effect  we  next 
ex6unine  the  results  for  a free  stream  velocity  of  25  ft/sec, 
which  corresponds  to  a Reynolds  number  based  on  the  wire 
diameter  of  1.8.  The  ratio  of  improvement  over  the  cosine 

■ law,  a*/a*,  is  shown  in  Figure  39,  where  intermediate 

c p 

' value  points  were  determined  using  a three  point  Lagrangian 

I interpolation  scheme  (available  with  HP97  and  67  calculators) 

in  order  to  fit  a smooth  curve  through  the  data. 
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For  the  intermediate  range  of  sensor  t/d  (f/d=250) , 
this  ratio  reaches  a maximvim  value  of  over  20.  This  maxi- 
mum value  is  reduced  by  either  increasing  or  decreasing 
the  sensor  l/d.  Surprisingly,  the  ratio  of  improvement 
is  smaller  for  the  very  short  sensor  (Vd*75)  than  for  the 
long  sensor  (^d^lOOO)  . This  Indicates  that  the  cosine 
law  is  more  applicable  to  the  short  sensor,  rather  than  to 
the  long  sensor.' 

Figure  40  shows  the  variation  of  oj/o*  with  l/d. 

Once  again  this  ratio  reaches  large  values  for  the  inter- 
mediate sensor  lengths.  Since  the  present  relation  is  a 
more  complex  form  of  the  Friehe  and  Schwarz  relation. 


one  may  ask  ; when  does  the  benefits  of  greater  accuracy 
outweigh  the  additional  complexity?  If  we  arbitrarily 
assign  a ratio  of  improvement  of  two  as  a marginal  addi- 
tional accuracy,  remembering  that  in  all  cases  a**0.15%, 
then  regions  where  oJ/o*>2  would  be  regions  of  significant 
improvement.  Using  this  criterion,  one  finds  that  signi- 
ficant improvement  can  be  obtained  in  going  from  the 
Friehe  and  Schwarz  relation  to  the  present  one,  if  sensors 
with  -f/d's<550  are  used. 

The  dependence  of  b|  on  sensor  l/d  is  examined  in 
Figure  41.  The  variation  indicates  that  no  significant 
range  of  sensor  l/d*s  fits  the  cosine  law.  Note  that  b| 
can  take  on  values  greater  than  or  less  than  one  indi- 
cating negative  and  positive  deviations  from  the  cosine 
law,  respectively.  Figure  41  indicates  that  positive 


deviations  from  the  cosine  law  occur  for  the  intermediate 
range  of  sensor  f/d. 
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The  variation  of  b*  and  m*  with  sensor  f/d  is  shown  | 

i 

in  Figures  42  and  43 , respectively . The  value  of  b*  i 

reaches  a minimum  for  f/d»250  and  increases  while  either 
increasing  or  decreasing  the  sensor  f/d.  From  Figure  43, 
the  exponent  m*  has  the  largest  deviation  away  from  the 
Friehe  and  Schwarz  exponent  of  0.5  in  the  region  of  sig- 
nificant improvement  as  shown  in  Figure  40. 

The  turbulence  correction  factor,  as  determined  from  | 

Friehe  and  Schwarz  and  the  present  relation, plotted  as  a | 

function  of  f/d,  is  shown  in  Figure  44.  This  figure  indi- 
cates that  all  sensors  will  require  a turbulence  correc- 
tion factor.  However,  by  proper  selection,  this  can  be 
minimized.  We  note  again  as  in  Figures  39  through  43, 
that  large  deviations  from  the  cosine  law  exist  in  the 
intermediate  range  of  sensor  t/d.  Figure  44  also  reveals 
that  there  exists  only  a small  difference  between  (g/h)| 
and  (g/h)*.  However,  it  will  be  shown  in  Chapter  VII  that 
this  small  deviation  will  be  magnified  many  times  in  cases 
where  the  mean  flow  angle  y is  not  equal  to  zero  with 
respect  to  an  X-probe. 

Combined  Effects  of  Reynolds  Number  and  Sensor  t/d. 

The  combined  effects  of  Reynolds  number  and  sensor  t/d  are 
presented  in  Figures  45  through  50  where  the  various  para- 
meters Investigated  are  displayed  as  constant  value  con- 

i 

tours  in  the  Re^  versus  t/d  plane.  Contours  of  constant  j 
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ratio  of  improvement,  at/a*,  are  shown  in  Figure  44.  This 

^ It 

again  points  out  the  large  improvement  gained  in  going 
from  the  Friehe  and  Schwarz  relation  to  the  present  rela- 
tion for  the  intermediate  length  sensors,  in  particular 
at  low  Reynolds  numbers.  This  figure  also  depicts  the 
decrease  of  this  ratio  with  increasing  Reynolds  number, 
indicating  that  the  minimum  standard  deviation  found  using 
the  Friehe  and  Schwarz  relation  does  not  remain  constant, 
unlike  that  found  from  the  present  relation. 

The  variation  of  the  parameters  in  the  Friehe  and 
Schwarz  and  present  relations  is  shown  in  Figures  46 
through  48.  In  all  cases,  the  coefficients  are  strong 
functions  of  t/A.  In  general,  these  coefficients  are  also 
Reynolds  number  dependent.  However  it  appears  possible  to 
find  a probe  in  which  these  constants  are  weak  functions 
of  Reynolds  number.  It  is  most  significant  to  point  out 
that  the  trend  of  the  parameters  with  either  of  the  two 
variables  does  not  seem  to  approach  the  cosine  law,  but 
rather  it  seems  to  overshoot  it  without  any  observed  asymp' 
totic  values. 

The  same  can  also  be  said  for  the  turbulence  correc- 
tion factors  (g/h)*  and  (g/h)J,  namely,  that  they  are 
functions  of  both  sensor  l/d  and  Reynolds  nximber,  as  dis- 
played in  Figures  49  and  50. 

Application  of  Yaw  Relations  Over  Reynolds  Number 
Range.  As  a result  of  the  Reynolds  number  dependence  of 
the  yaw  sensitivity  coefficients,  one  is  lead  to  ask  the 
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following  questions  : does  a set  of  parameters  exist  for 
either  of  the  yaw  relations,  which  can  be  used  over  a 
Reynolds  number  range  and  which  would  give  the  best  accu- 
racy over  this  range?  If  such  a set  of  parameters  exists, 
how  good  is  the  accuracy  of  the  relation  in  representing 
the  data  over  this  range  of  Re^?  To  answer  these  ques- 
tions for  a particular  sensor,  the  values  of  the  standard 
deviation  calculated  in  the  b-m  plane  were  averaged  over 
the  four  Reynolds  numbers  at  each  grid  point  in  the  plane. 
This  value  is  denoted  Op.  The  standard  deviation  of  the 
four  values  from  this  value  is  denoted  by  Op.  The  sum  of 

and  will  then  indicate  an  overall  effective  standard 
P P 

deviation,  denoted  by  Z . This  effective  standard  devia- 

P 

tion  is  considered  a measure  of  the  quality  of  the  fit 
over  the  prescribed  Reynolds  number  range. 

Contours  of  constant  Z 's  for  an  f/d=250  are  shown  in 

P 

the  b-m  plane  of  Figure  51.  These  curves  have  a shape 
similar  to  those  at  a single  Re^,  however  the  magnitudes 
are  larger.  Once  again  we  see  that  the  minimum  region  is 
far  removed  from  the  cosine  law  although  there  is  negli- 
gible improvement  over  the  Friehe  and  Schwarz  relation 
when  the  present  relation  is  utilized.  Results  for  a sen- 
sor with  f/d=580  are  reproduced  in  Figure  52.  In  this 
case  there  is  a negligible  difference  between  all  three 
yaw  relations  including  the  simple  cosine  law. 

This  is  an  expected  result,  since  based  on  the  ex- 
perimental data  the  variation  of  U^/U^  with  yaw  angle 
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strongly  depends  on  Reynolds  number.  At  any  prescribed 
value  of  Re^f  the  present  relations  will  most  accurately  de- 
scribe the  yaw  dependence,  but  with  different  values  of  the 
parameters  b*  and  m*.  Small  changes  in  these  parameters  may 
lead  to  very  large  changes  in  the  errors  introduced  by  the 
relation  in  describing  the  data,  as  Indicated  by  the  sharp 
gradients  in  the  contours  of  in  the  b-n  plane.  Aver- 
aging over  any  number , n , of  Reynolds  numbers  will  then  lead 
to  significant  contributions  to  the  errors  from  at  least 
n-1  of  the  Reynolds  numbers.  Therefore,  if  the  turbulence 
velocities  and  correlations  are  calculated  using  these 
averaged  values,  significant  errors  will  result.  This  leads 
one  to  the  conclusion  that  a different  set  of  parameters 
must  be  utilized  in  the  yaw  relation  for  each  mean  velocity 
the  X-probe  is  used  at;  a most  severe  restriction. 

Effect  of  Overheat  Ratio.  The  effect  of  the  operating 
overheat  ratio  on  the  yaw  coefficients  was  also  examined 
using  Probe  2,  equipped  with  a standard  length  hot-wire.  Yaw 
calibration  was  carried  out  at  two  values  of  freestream 
velocity  and  three  different  values  of  overheat  ratio:  1.4, 

1.6  and  1.8,  the  latter  being  the  most  commonly  used  by  ex- 
perimenters. (Note  that  this  resulted  in  repeating  two 
runs  at  the  1.8  overheat  ratio  with  two  identical  probes 
using  similar  t/d  sensors) . These  results  are  described  in 
Figures  53  and  54.  In  general,  the  ratio  of  improvement 
a|/a*  decreases  with  decreasing  the  overheat  ratio,  while 
the  turbulence  correction  factor  remains  constant  within  2%. 
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Figure  54  proves  that  the  parameters  b*,  m*  and  b|  vary 
slightly  with  overheat  ratio;  b*  decreasing  with  decreasing 
overheat  ratio;  m*  increasing  with  decreasing  overheat 
ratio;  and  b|  increasing  with  decreasing  overheat  ratio. 

All  of  these  results  indicate  that  the  yaw  sensitivity  of 
wires  is  not  significantly  influenced  by  the  operating  over- 
heat ratio. 

The  values  of  the  parameters  of  the  yaw  relations  ob- 
tained using  Probe  2 with  a standard  length  sensor  operated 
using  an  overheat  ratio  of  1.8  can  be  compared  to  the  re- 
sults for  a standard  length  sensor  using  Probe  1 presented 
in  earlier  parts  of  this  chapter  in  connection  with  Fig- 
ures 46  through  48.  This  comparison  can  serve  as  a repeat- 
ability check  on  the  experimental  procedure  and  may  shed 
light  on  the  generality  of  the  parameters.  The  values  of 
bj  and  m*  from  Figure  54  are  nearly  identical  to  those  in 
Figure  46  and  48  respectively.  The  variation  of  b*  with 
Re^  in  Figures  47  and  54  is  also  identical,  however  the 
values  are  about  15%  higher  in  Figure  54 . A comparison  of 
the  turbulence  correction  factor  (g/h) * in  Figures  50  and 
53  reveals  that  the  values  are  identical.  Therefore,  the 
two  sets  of  experiments  demonstrate  that  the  results  are 
repeatable  using  the  same  probe  and  sensor  design.  How- 
ever, at  this  point  one  cannot  reach  any  conclusions  re- 
garding the  generality  of  the  yaw  dependence  parameters. 

Effect  of  Angular  Range  of  Data.  The  mean  yaw  cali- 
bration data  discussed  in  this  chapter  were  taken  in  the 
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range  of  yaw  angles  from  -70®  to  +70®.  The  effect  of 
changing  this  range  on  the  yaw  coefficients  is  demonstrated 
in  this  section.  Results  for  a sensor  with  an  t/d  of  145 
are  shown  in  Figures  55  and  56.  For  this  particular  probe 
the  value  of  the  turbulence  correction  factor  (g/h) * is 
independent  of  the  angular  range  of  fit,  while  the  ratio  of 
improvement,  o|/a*,  decreases  to  unity  for  all  Reynolds 
numbers  when  the  angular  range  of  data  is  +45®.  This 
means  that  the  improvement  in  using  the  present  relation 
over  the  Friehe  and  Schwarz  relation  comes  at  relatively 
large  yaw  angles  and  particularly  at  low  Reynolds  nxunbers. 

The  variation  of  b*  and  m*  is  shown  in  Figure  56. 

These  coefficients  remain  constant  from  an  angular  range  of 

+50®  to  +70®.  However,  the  value  of  b|  increases  linearly  | 

as  the  angular  range  of  the  data  is  decreased. 

A similar  variation  in  these  quantities  is  observed 
for  a sensor  of  f/d=580  as  depicted  in  Figures  57  and  58. 

Here,  (g/h)*  remains  constant  and  CT|/a*  decreases  to  unity 
as  Aa  decreases.  However,  the  variations  of  b*  and  m*  are 
not  as  well  defined  as  in  the  previous  case:  e.g.,  b| 
remains  constant  with  Aa  to  within  1% . 


X-Wire  Probe 

In  measuring  turbulence  velocity  correlations  using 
X-probe,  the  time  dependent  component  of  the  sensor  output 
voltages  must  be  related  to  the  fluctuating  velocity  com- 
ponents. One  method  has  been  examined  In  Chapter  III,  where 
expressions  accounting  for  the  deviations  from  the  cosine 
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law  are  derived.  Thus,  once  the  probe  has  been  calibrated 
for  yaw  sensitivity,  these  relations  can  be  analytically 
evaluated.  Many  investigators  do  not  employ  this  method 
but  rather  calibrate  the  X-probe  in  a manner  such  that  the 
deviations  from  the  cosine  law  are  accounted  for  in  a cali- 
bration constant.  This  section  will  compare  these  two 
methods  of  determining  the  longitudinal  and  transverse  ve- 
locity sensitivity  coefficients.  A brief  description  of 
the  direct  calibration  method  follows. 

Determination  of  Transverse  and  Longitudinal  Sensi- 
tivity Coefficients . If  the  anemometer  output  voltage  is 
linearized,  then  one  can  write,  [22,35] 


E = E(Q^,e) 


so  that 


(V-1) 


3^  dQi  + ^ dB 


(V-2) 


By  use  of  Equations  (III-24) , (111-75)  and  (III-78) , this 
becomes 


• -[isjs  - 5;[||]  «n 

Only  a comparison  between  the  calculated  and  calibrated 
values  of  the  transverse  sensitivity,  OE/3a)/Q_,  is  attempt- 
ed  in  the  following.  However,  similar  results  for  the 
longitudinal  sensitivity  coefficients  can  be  obtained. 
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From  the  mean  yaw  calibration  data  of  the  X-probe,  the 
velocity  as  a function  of  yaw  angle  was  computed  so  that  the 
coefficient  in  Equation  (V-3)  could  be  evaluated.  In  order 
to  obtain  this  coefficient^  methods  of  determining  local  de~ 
rivatives  of  the  data  had  to  be  used.  In  order  to  minimize 
any  effects  of  the  differentiation  scheme  on  the  results, 
four  different  methods  of  determining  the  local  derivatives 
were  employed.  The  first  method  utilized  a two  point  local 
derivative  scheme  where  3E/9a  is  approximated  by 

E(a.+h)  - E(a.-h) 

— jiT ^ <V-4) 

where  h is  the  a increment  between  the  data  points.  This 
method  should  be  the  most  susceptible  to  the  scatter  in  the 
data.  The  second  method  utilized  a four  point  local  deri- 
vative where 


3E 


3E 

7K 


la«ot 

I 


8 [E (a^+h) -E (a^-h) ] -E (a^+2h) +E (a^-2h) 

12h 


(V-5) 


Another  method  involved  performing  a linear  least  squared 
fit  on  five  adjacent  data  points,  1*  apart,  to  determine 
the  slope  of  the  line,  thereby  obtaining  the  derivative  at 
the  mid-point.  The  last  method  consisted  of  fitting  the 
same  five  points  to  a second  order  polynomial  of  Chebyshev 
type,  from  which  the  derivative  at  the  midpoint  could  be 
determined  analytically. 
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A complete  yaw  calibration,  similar  to  those  presented 
for  the  single  sensor,  was  carried  out  at  four  different 
mean  velocities  for  each  of  the  sensors  of  the  X-probe.  A 
finer  grid  of  data  was  collected  around  the  axis  of  sym- 
metry of  the  X-wires.  While  data  were  obtained  every  three 
degrees  for  -70*<a<70*,  where  a is  measured  from  the  axis  of 
symmetry,  the  calibration  was  performed  every  one  degree 
over  the  range  -10<a<10*.  Results  for  the  transverse  sensi- 
tivity are  shown  in  Figures  59  through  63.  The  results  for 
the  parameters  of  the  yaw  relation  are  listed  in  tables  in 
the  same  figures  and  discussed  in  the  last  part  of  this 
chapter  and  in  Chapter  VII.  Figure  59  compares  the  varia- 
tion of  the  four  differentiation  schemes  with  the  values 
calculated  from  Equation  (111-35) . The  method  which  pro- 
duces the  largest  scatter  in  the  results  is  the  two  point 
local  differentiation  method.  The  five  point  least  squares 
method  averages  out  most  of  the  scatter  in  the  data.  The 
values  computed  from  the  Friehe  and  Schwarz  and  the  present 
relations  lie  through  the  values  obtained  by  direct  cali- 
bration. This  indicates  that  computing  the  velocity  sensi- 
tivity coefficients  through  the  Fpiehe  and  Schwarz  relation 
or  the  present  relation,  using  the  present  or  a similar 
technique,  yields  more  accurate  results.  The  utilization 
of  any  "sophisticated"  yaw  relation  will  lead  to  the  same 
result  since  it  will  not  be  sensitive  to  small  scatter  in 
the  data  which  is  magnified  when  local  differentiation  is 
applied  or  when  the  experimenter's  judgement  is  required 
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for  a different  differentiation  scheme,  e.g.,  graphical 
methods.  We  also  observe  that  the  cosine  law  does  not 
accurately  describe  the  velocity  sensitivity  coefficients. 
It  should  be  noted  that  the  sensor  used  in  the  example  of 
Figure  59  lies  in  the  region  of  minimum  improvement  over 
the  Friehe  and  Schwarz  relation  thereby  giving  negligible 
differences  in  the  computed  sensitivities  between  the  pres- 
ent relation  and  their 's.  If  an  intermediate  length  sensor 
was  utilized,  the  differences  would  be  more  significant. 
This  will  be  discussed  in  Chapter  VII. 

Examples  using  different  mean  free-stream  velocities 
and  the  other  wire  in  the  X-probe  are  shown  in  Figures  60 
through  63.  Only  the  two  point  differentiation  method  is 
used  in  these  figures  since  based  on  the  comparison  of 
Figure  59  the  difference  between  the  methods  of  obtaining 
local  derivatives  from  the  present  data  is  considered  in- 
significant. In  each  case  the  values  calculated  from  the 
yaw  sensitivity  function,  taking  into  account  deviations 
from  the  cosine  law,  yield  a higher  degree  of  accuracy  than 
those  obtained  from  the  direct  calibration  method.  When 
the  best  fit  parameters  of  the  yaw  relation  are  near  those 
corresponding  to  the  cosine  law,  the  cosine  derivatives 
appear  to  be  just  as  good;  see  Figures  60  and  61. 

Recalling  that  the  design  and  spacing  of  the  prongs 
utilized  in  the  X-probe  were  different  than  those  for  the 
single-wire  probe,  a comparison  of  the  yaw  coefficients  for 
a standard  length  sensor  (Figures  46  through  48)  can  be 
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made  with  the  values  presented  in  the  tables  of  Figures  59 
through  63.  The  Reynolds  number  dependency  is  identical  . 
However, the  values  of  the  coefficients  are  substantially 
different.  These  results  indicate  that  the  yaw  relation 
parameters  are  strongly  influenced  by  the  probe  design  and 
construction.  Therefore,  each  probe  must  be  calibrated  to 
determine  the  yaw  parameters  associated  with  it  at  the 
Reynolds  numbers  of  interest. 


t 


88 


CHAPTER  VI 

YAW  CALIBRATION  FOR  FLUCTUATING  VELOCITIES 

When  a sensor  is  calibrated  to  determine  its  yaw  sensi- 
tivity/ the  effective  cooling  of  the  wire,  as  represented  by 
the  anemometer  output  voltage#  is  recorded  as  a function  of 
angle  during  a static  calibration  procedure  utilizing  a low 
turbulence#  calibration- type  flowfleld.  The  results  of  the 
calibration  are  represented  by  a relation  like  Equation 
(III-80) . However#  when  the  sensor  is  used  for  measuring 
turbulence  velocity  correlations#  it  is  exposed  to  fluc- 
tuating velocities  of  varying  degrees  of  intensity  and  fre- 
quency distributions.  This  chapter  examtines  the  effect  of 
such  fluctuations#  as  represented  in  particular  by  control- 
led periodic  velocities#  on  the  yaw  sensitivity  of  the  sen- 
sor, which  is  given  by  the  relation  between  the  time- 
averaged  effective  cooling  and  the  mean  flowfleld.  It  also 
examines  the  validity  of  Equation  (111-83) # which  was  de- 
rived to  relate  the  yaw  dependence  of  the  mean  and  unsteady 
fluctuations  and  was  based  on  a quasi-steady  approximation 
of  the  yaw  relation. 

Establishment  of  Fluctuating  Velocity  Field 

The  calibration  tunnel  was  initially  modified  to  super- 
impose streamwise  periodic  fluctuations  on  the  mean  flow  in 
the  manner  described  by  Marcichow  and  Way  [ 30] . In  this 
configuration#  a linear  contraction  leading  to  a one-inch 


diameter  nozzle  was  attached  through  a short  duct  segment 
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to  the  open  end  of  the  plexiglass  "T"  section  (see  Chap- 
ter IV  and  Figures  9 and  10)  . Driven  by  the  spea)cer  at  dif- 
ferent frequencies,  the  resulting  free- jet  flowfield  pro- 
vided a fluctuating  velocity  field  with  a controllable  air 
speed  and  amplitude  of  fluctuations.  In  particular,  high 
levels  of  intensity  of  the  fluctuations  could  be  achieved. 

A standard  hot-wire  attached  to  the  angularity  jig  and 
indexing  head  was  used  to  carefully  map  the  flowfield  pro- 
duced by  this  configuration  at  different  values  of  free 
stream  velocity,  forcing  frequency  and  fluctuations  inten- 
sity. Vertical  and  lateral  profiles  at  various  positions 
downstream  of  the  jet  were  taken  utilizing  the  instrumenta- 
tion of  Figure  14.  These  profiles  were  of  the  mean  velo- 
city, the  rms  of  the  educted  fluctuations  and  the  phase  dif- 
ference between  the  reference  signal,  from  the  beat  frequen- 
cy oscillator,  and  the  educted  wire  signal.  For  this  con- 
figuration, it  was  determined  that  the  phase  changed  by  as 
much  as  40**  from  the  jet  centerline  to  the  edge  of  the  jet 
at  a position  0.25  inch  downstream  of  the  jet  exit.  This 
phase  difference  increased  with  downstream  distance  and  was 
found  to  be  asymmetric.  Initially  this  was  thought  to 
be  due  to  misalignment  of  the  speaker-duct  connection.  The 
speaker  was  then  inclined  at  various  angles  to  the  duct  and 
the  measurements  repeated.  The  second  set  of  measurements 
indicated  that  the  angle  of  the  speaker  had  little  influence 
on  the  flow  field  out  of  the  nozzle,  and  that  the  contrac- 
tion and  nozzle  were  the  likely  component  contributing  to 
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the  skewness  of  the  unsteady  flowfield. 

Considering  the  size  of  the  probe  with  respect  to  the 
nozzle  exit,  one  can  surmise  that  such  phase  differences 
would  lead  to  fluctuations  incident  upon  the  sensor  which 
are  not  planar.  This  in  turn  would  lead  to  inaccuracies 
in  the  measurements,  in  particular  at  large  yaw  angles. 

In  an  attempt  to  minimize  these  effects,  the  linear  con- 
traction was  replaced  by  a smooth  contraction  section 
leading  to  an  identical  size  free  jet.  (The  same  contrac- 
tion section  is  utilized  when  the  tunnel  is  used  for  velocity 
calibration  of  probes;  see  Chapter  IV  and  Reference  [28]). 
After  careful  mapping  of  the  flowfields  produced,  it  was 
found  that  the  phase  difference  was  reduced  to  approxi- 
mately 20*.  This  was  still  not  acceptable,  in  particular 
since  the  rms  of  the  fluctuations  was  not  constant  over  the 
test  volume  in  which  the  sensor  would  be  positioned  during 
yaw  calibration. 

The  main  advantages  of  the  above  two  configurations  are 
the  low  background  turbulence  level  in  the  free  jets,  and 
the  wide  range  of  velocities  and  fluctuation  intensities 
achievable.  However,  the  two-dimensionality  of  the  fluc- 
tuating velocity  field  was  considered  the  most  important 
aspect  and  the  configuration  of  the  Calibration  Tunnel  was 
modified  to  achieve  it.  The  final  modification  consisted 
of  the  removal  of  the  contraction  section  leading  to  the 
nozzle,  thereby  having  the  duct  open  to  the  laboratory.  By 
mapping  the  flow  in  this  configuration,  it  was  found  that 
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as  long  as  the  sensor  remained  within  the  duct  ,the  velocity 
fluctuations  were  planar  over  the  entire  cross  section. 

Assured  that  the  fluctuations  in  the  duct  were  two 
dimensional,  turbulence  manipulators  were  added  to  reduce 
the  level  of  the  turbulence  to  an  acceptable  level  but  with 
minimum  interference  to  the  unsteady  flowfield.  The  final 
arrangement  of  these  is  shown  in  Figure  10.  Without  peri- 
odic forcing  of  the  flow,  at  a freestream  velocity  of 
12  ft/sec,  the  turbulence  intensity  near  the  exit  of  the 
duct  was  found  to  be  2%.  For  more  details  on  the  flow  mani- 
pulators and  the  flow  characteristics  at  the  test  position, 
including  velocity  and  turbulence  profiles  and  spectra, 
the  reader  is  directed  to  the  report  of  Loehrke  and 
Nagib  [29l . 

Additional  experiments  were  conducted  to  select  the 
combinations  of  mean  velocity,  fluctuations  amplitude  and 
frequency  suitable  for  calibration.  All  the  selected  con- 
ditions exhibited  two  dimensional  fluctuations  across  the 
test  section  and  simple  harmonic  character  of  the  unsteady 
j velocity  component.  The  only  limitations  were  the  maximum 

I steady  velocity  achievable  and  the  maximum  intensity  of 

fluctuations  which  could  be  obtained  in  this  configuration 
of  the  tunnel.  However,  as  will  be  evident  from  the  fol- 
lowing sections,  the  available  conditions  were  adequate  for 
the  purposes  of  the  experiments. 

In  the  remainder  of  this  chapter , single  and  X-probes 
are  tested  in  the  fluctuating  velocity  field  described 
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above.  It  should  be  emphasized  at  this  point  that  the  sen- 
sors are  exposed  simultaneously  to  periodic  and  random 
fluctuations  and  that  their  cooling  is  affected  by  both 
unsteady  velocity  fields.  However,  with  the  aid  of  the 
periodic  sampling  techniques  utilizing  the  Waveform  Eductor 
(for  details  see  Chapter  IV) , one  can  separate  the  periodic 
part  of  the  hot-wire  output  signal.  This  periodic  part 
will  be  useful  in  the  comparison  of  the  yaw  relations  with 
the  aid  of  equations  given  in  Chapter  III. 

Single-Wire  Probe 

A linearized  standard-length  hot  wire  was  used  in  the 
modified  Calibration  Tunnel  as  described  above  and  in  Chap- 
ter IV.  The  velocity  fluctuations  presented  here  refer  to 
the  educted  signal  from  the  wire,  i.e.,  the  linearized  hot- 
wire output  was  periodically  averaged  to  provide  the  rms  of 
the  periodic  component  of  the  measured  velocity.  With  the 
free-stream  velocity  set  at  12  ft/sec,  the  forcing  frequency 
set  at  lOOOHz  and  the  fluctuations  intensity  equal  to  3.8%, 
the  effective  mean  and  rms  velocities  were  recorded  as  the 
yaw  angle  was  varied.  The  results  shown  in  Figure  64, 
for  various  angles  of  yaw  in  the  range  +48®,  indicate  that 
the  effective  turbulence  intensity,  v^/Ug,  remains  constant. 
The  behavior  of  the  data  in  this  figure  is  represented  by 
Equation  (III-83)  which  testifies  to  its  validity. 

To  examine  this  in  greater  detail,  forcing  frequencies 


of  different  intensities  were  used.  The  effective  turbu- 
lence intensity,  normalized  by  its  value  when  the  sensor  is 
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normal  to  the  f low  ,is  shown  in  Figure  65  as  a function  of 
yaw  angle  for  different  frequencies  and  intensities.  A wide 
range  of  frequencies  encompassing  a major  portion  of  the 
hot-wire  spectrum  is  used.  In  all  cases,  the  turbulence 
intensity  remains  constant  and  within  2%  or  less  of  its 
value  at  zero  yaw  angle.  With  these  results  the  validity  of 
Equation  (III-83)  is  confirmed. 

Having  established  that  at  any  forcing  frequency,  the 
yaw  sensitivity  for  the  mean  and  rms  fluctuations  of  the 
flow  are  the  same,  the  effect  of  the  velocity  fluctuations 
on  the  parameters  of  the  yaw  relation  based  on  mean  calibra- 
tion is  examined  next.  Maintaining  the  free-stream  velocity 
constant  while  varying  the  frequency,  the  yaw  sensitivities 
for  both  mean  and  rms  velocities  were  obtained.  For  ease  of 
comparing  the  results,  only  the  variation  in  b|  will  be  pre- 
sented here.  This  is  done  because  the  Friehe  and  Schwarz 
relation  has  only  one  variable  coefficient  and  its  variation 
in  this  case  should  be  easier  to  follow.  These  results  are 
shown  in  Figure  66.  In  absence  of  the  velocity  fluctuations, 
the  value  of  b|  was  determined  to  be  equal  to  1.0.  When 
periodic  fluctuations  were  superimposed  on  the  flow,  keeping 
the  time-mean  velocity  constant,  the  values  of  b|  obtained 
from  both  mean  and  rms  yaw  calibration  changed  by  as  much  as 
3%  at  the  lower  frequencies  and  by  as  little  as  1%  at  the 
higher  frequencies.  The  figure  displays  b|,  obtained  in 
presence  of  the  fluctuation,  decreasing  with  increasing 
frequency  and  approaching  the  value  obtained  from  static 
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calibration , it  should  be  pointed  out  that  the  same 
amplitudes  of  fluctuations  were  used  in  Figures  65  and 
66. 

This  frequency  dependence  of  the  yaw  coefficients  will 
cause  a discrepancy  between  static  and  dynamic  calibration 
for  the  X-probe  sensitivity  coefficients  described  in  Equa- 
tion (III-35) . In  addition,  Figure  €6  brings  out  other  im- 
portant problems.  If  mean  velocities  are  to  be  measured  in 
an  essentially  laminar  flow  using  inclined  wires,  only  a 
static  calibration  will  describe  its  true  yaw  behavior.  In 
order  to  accurately  measure  mean  velocities  with  an  X-wire 
in  a turbulent  flow  the  probe  should  be  calibrated  in  a fluc- 
tuating flowfield  of  similar  intensities  and  frequencies. 

In  cases  where  fluctuating  velocities  are  to  be  measured 
using  inclined  wires,  neither  static  nor  dynamic  calibration 
will  yield  accurate  results  since  the  sensitivities  for  in- 
clined wires  will  be  frequency  dependent.  This  frequency 
dependence  can  also  lead  to  substantial  errors  in  spectral 
measurements  utilizing  inclined  sensors. 

Although  the  data  of  Figure  66  are  based  on  the  period- 
ic part  of  the  hot-wire  signal,  the  results  may  be  influ- 
enced by  the  "presence  of  the  random  2%  turbulence  fluctua- 
tions. Based  on  Equation  (111-79)  and  in  presence  of  turbu- 
lence Reynolds  stresses,  uv,  one  would  expect  the  difference, 
depicted  in  Figure  66,  between  the  mean-velocity  yaw  cali- 
bration with  and  without  fluctuations  present.  A similar 
conclusion  can  be  made  regarding  the  unsteady  velocity  from 
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Equation  (111-81) . In  the  calibration  condition  used  in  ■ 

> this  chapter  such  Reynolds  stresses  could  be  present  in  the  j 

! 

random  fluctuations  or  due  to  their  interaction  with  the 
periodic  velocity  field.  However  it  remains  somewhat  puz- 
zling that  Equation  (III-83)  was  verified  earlier,  while  in 
view  of  the  results  of  Figure  66  it  may  not  necessarily 
hold.  It  could  be  that  the  additional  terms  in  Equation 
(111-79)  and  (III-81)  balanced  each  other  leading  to  the 
behavior  described  by  Equation  (III-83)  without  having  to 
satisfy  all  the  assumptions  made  to  derive  it. 

X-Wire  Probe 

The  applicability  of  Equation  (III-83)  was  also  exam-, 
ined  through  measurements  utilizing  a linearized  X-wire 
probe.  Results  for  the  normalized  turbulence  intensity  for 
different  ranges  of  angular  position  are  shown  in  Figures  67 
and  68  for  frequencies  of  different  intensities.  Once 
again,  these  show  that  the  effective  turbulence  intensity 
remains  essentially  constant.  The  scatter  in  the  data  in 
Figure  68  arises  from  measurements  with  forcing  frequencies 
having  intensities  of  1 to  2%,  where  a small  error  in  mea- 
suring the  rms  velocity  leads  to  large  errors  in  the  mea- 
sured turbulence  intensity. 

The  transverse  sensitivity  coefficient,  numerically 
calculated  by  a two  point  local  derivative  method  (for  fur- 
ther details  see  the  section  of  "X-Wire  Probe"  in  Chapter  V) , 
for  mean  and  fluctuating  velocities,  using  a forcing  frequen- 
cy of  40  Hz,  is  shown  in  Figure  69.  These  should  be  equal 
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through  use  of  Equation  (III-83)  . The  transverse  sensitivi- 
ty coefficient  determined  from  Equation  (III-35)  using  the 
yaw  coefficients  obtained  from  yaw  calibration  in  the  ab- 
sence of  periodic  velocity  variations,  is  also  shown  on  the 
same  figure.  A significant  difference  in  the  calculated 
sensitivity  coefficients  exists  between  static  and  dynamic 
calibration.  The  reason  for  this  lies  in  the  frequency 
dependence  of  the  yaw  coefficients  as  described  in  Figure  66. 
Although  reasonable  agreement  is  demonstrated  between  the 
three  curves  of  Figure  69,  other  data  did  not  yield  as 
good  an  agreement.  However,  in  all  cases  larger  discrepan- 
cies were  found  between  the  static  and  dynamic  calibration 
for  the  yaw  dependence  of  mean  velocities. 

Morrison  et  al.  [26]  found  that  large  differences  in 
both  sensitivity  coefficients  of  the  X-wire  existed  between 
static  and  dynamic  calibration.  Bruun  [27]  argues  that 
part  of  this  difference  is  caused  by  poor  approximation  of 
the  static  calibration  data  over  a large  velocity  range  by 
use  of  a constant  power  exponent  law,  and  seems  to  attribute 
all  of  it  to  this  approximation.  However,  his  results  con- 
sistently show  the  statically  determined  longitudinal  sen- 
sitivity coefficient  less  than  that  determined  by  dynamic 
calibration  by  approximately  2%.  This  is  consistent  with 
the  present  results  and  hence  could  not  be  attributed  to 
experimental  error  as  Bruun  states. 
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CHAPTER  VII 
DISCUSSION 

Temperature  Dependence 

A scheme  utilizing  a temperature  probe  immersed  in  the 
working  fluid  to  compensate  for  the  dependence  of  the  hot- 
wire velocity  calibration  on  ambient  temperature  variations 
was  exaunined  in  Chapter  II.  Requiring  only  the  knowledge  of 
the  resistances  of  the  velocity  sensor  and  the  temperature 
probe,  and  their  temperature  dependence,  this  scheme  was 
found  capable  of  compensation  over  as  much  as  40"C  with  as 
small  as  a 1%  error  in  the  indicated  velocity  (The  error  is 
essentially  proportional  to  the  square  of  the  maximum  tem- 
perature variations) . This  error  may  be  reduced  even  fur- 
ther if  the  optimum  circuit  design  is  employed.  The  anal- 
ysis for  this  optimum  design  yields  a range  of  solutions 
from  which  a minimum  error  condition  can  be  derived,  in 
principle.  This  derivation  requires  an  optimization  ap- 
proach and  has  not  been  completed  yet. 

The  analysis  and  error  estimates  were  based  on  the  pre- 
mise that  the  coefficient  A in  King's  law  is  independent  of 
temperature.  However,  it  was  determined  in  Appendix  B that 
this  coefficient  exhibits  a slight  variation  with  ambient 
temperature.  This  variation  will  add  a small  contribution 
to  the  errors  in  the  compensation  attained. 

If  one  wishes  to  utilize  this  scheme  without  velocity 
calibration  at  various  temperatures,  one  may  use  the 


t 


f 


J 


98 


I 


analysis  based  on  the  constant  A assumption.  However,  if 
one  uses  the  scheme  after  such  velocity  calibration  has  been 
performed,  the  analysis  can  be  modified  slightly  to  take 
into  account  the  variation  of  A with  temperature.  The 
errors  can  be  determined  from  the  data  presented  in 
Appendix  B.  They  will  be  calculated  by  comparing  the  values 
based  on  the  experimental  data  of  dR^/dT,  i.e.,  the  values 
required  for  a constant  anemometer  output  at  a fixed  velo- 
city with  varying  temperature  , with  those  derived  by  the 
analysis  of  Chapter  II. 

Yaw  Dependence 

In  regard  to  the  yaw  dependence  of  hot-wires,  results 
for  yaw  calibration  in  steady  velocities  are  presented  in 
Chapter  V.  As  a check  on  these  results,  they  are  compared 
to  those  available  in  the  literature. 

Since  the  Friehe  and  Schwarz  relation  was  one  of  the 
yaw  relations  examined  in  the  course  of  the  present  investi- 
gation, the  yaw  coefficients  determined  when  using  it  can  be 
compared  to  the  results  of  Friehe  and  Schwarz  [22] . This 
comparison  is  summarized  in  Figure  70.  Although  there  is  a 
fair  amount  of  scatter,  due  to  the  Reynolds  number  depend- 
ence, a general  trend  is  drawn.  The  data  of  Friehe  and 
Schwarz  [22]  fall  along  this  trend  indicating  that  the  re- 
sults are  consistent.  However,  the  present  data,  having 
cases  when  b|>l,  i.e.  negative  deviations  from  the  cosine 
law,  reveals  one  aspect  which  they  did  not  report. 

Another  check  on  the  data  is  performed  by  comparing  the 
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results  to  the  data  of  Champagne  [17]  . Champagne  used  the 
following  yaw  relation 
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(VII-1) 


where  k is  a small  parameter  accounting  for  the  axial 

cooling  of  the  wire.  The  variation  of  k with  sensor  t/d 

at  velocities  approximately  equal  to  110  ft/sec  was  reported 

in  his  work  [17,25].  Data  from  the  present  relation  were 

2 

converted  to  an  equivalent  average  value  of  k through  the 
following  relation 
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where  Aa  is  the  angular  range  over  which  the  yaw  calibration 

data  were  obtained.  These  results,  plotted  in  Figure  71, 

indicate  that  k^  is  a strong  function  of  Reynolds  Number, 

i.e.,  decreasing  with  increasing  Re^;  a result  which  is  also 

supported  by  Kjellstrom  and  Hedberg  [19]  . In  spite  of  this 

Reynolds  number  dependence,  a general  trend  with  sensor  l/d 

can  be  also  drawn  in  this  figure.  For  intermediate  sensor 
2 

lengths,  k is  positive.  However,  it  becomes  negative  for 

shorter  sensors  (i.e.  f./d<150)  and  for  longer  sensors 

(i.e.  i/d>600) , a result  not  determined  by  Champagne.  In 

all  cases  one  observes  that  Champagne's  data  follows  the 

trend  of  the  present  results.  Again,  as  for  b|>l,  the 

2 

negative  values  of  k represent  effective  cooling  of  the 
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sensor  less  than  that  described  by  the  cosine  law.  Fig- 

2 

ure  71  is  plotted  using  k instead  of  k to  avoid  intro- 
ducing imaginary  values  in  the  coefficients  of  the  yaw  rela- 
tions. However,  if  one  is  to  use  Equation  (VII-1)  to  repre- 
sent experimental  data,  the  present  results  indicate  that 

2 

one  should  either  use  k as  the  parameter  or  should  allow 
k to  be  imaginary  if  it  is  to  be  utilized  directly.  In  Fig- 
ure 71  each  of  the  present  data  sets  is  represented  twice: 
once  based  on  its  fit  to  the  present  relation,  and  the  other 
based  on  the  best  correlation  of  the  data  by  the  Friehe  and 
Schwarz  relation.  The  standard  deviation  of  the  two  fits  is 
different  and  is  always  smaller  in  the  former  case.  The 
trends  remain  the  same  in  both  representations. 

Based  on  the  results  of  Chapter  V and  in  Figures  70  and 
71,  it  is  evident  that  as  the  velocity  over  the  sensor 
changes,  its  yaw  sensitivity  is  altered.  At  any  velocity, 
for  a constant  temperature  anemometer,  the  conduction  from 
the  wire  to  the  prongs  is  constant.  However,  as  the  velo- 
city increases,  the  convection  increases  thereby  changing 
the  ratio  between  convection  and  conduction,  which  in  turn 
will  alter  the  axial  temperature  profile  of  the  wire. 

Horvatin  [36]  concludes  that  as  this  ratio  increases  (in 
his  case  by  decreasing  the  conduction  effects)  the  tempera- 
ture profile  becomes  more  uniform  along  the  entire  length 
of  the  wire.  This  can  be  also  calculated  directly  through 
use  of  the  analysis  presented  by  Sandborn  et  al.  (37) . 
Therefore,  one  would  expect  that  at  higher  velocities,  where  th 
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temperature  profile  becomes  more  uniform,  that  the  sensors' 
yaw  response  should  be  closer  to  the  cosine  law  approxima- 
tion. 

The  variation  of  the  yaw  coefficients  obtained  from  the 
present  relation  are  shown  in  Figures  72  and  73.  At  any 
l/d,  there  is  a wide  range  of  scatter  of  the  coefficients 
with  Reynolds  number  and  type  of  probe  used.  The  calibra- 
tion data  obtained  by  an  independent  source  over  a wider 
range  of  Reynolds  numbers  and  data  from  hot-films  cali- 
brated in  water  are  all  included  in  these  and  the  following 
figure.  Figures  72  and  73  indicate  that  the  yaw  coeffi- 
cients obtained  for  any  probe  are  strong  functions  of  probe 
design  and  construction.  Therefore,  each  probe  must  be  cali- 
brated separately  for  yaw  sensitivity.  Similar  results  for 
the  turbulence  correction  factor  are  presented  in  Figure  74, 
indicating  again  that  no  typical  values  of  the  yaw  para- 
meters can  be  found  for  any  of  the  f/d's  shown.  In  particu- 
lar, the  figure  displays  the  wide  range  of  correction  fac- 
tors of  the  turbulence  intensity  correlations  which  may  be 
encountered  for  typical  probes. 

Included  in  Figures  72,  73  and  74  are  the  data  from  the 
different  overheat  ratios,  the  different  angular  ranges  of 
calibration,  and  the  X-wire  calibration.  The  X-probe  data 
brings  out  a most  significant  result.  Comparing  the  values 
of  b|,  b*  and  m*  for  Wire  1,  listed  in  the  tables  of  Figures 
59,  60  and  61,  to  those  for  Wire  2,  given  in  Figures  62  and 
63,  one  notes  the  substantial  difference  between  them.  As 
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illustrated  by  the  photographs  of  Figure  12,  the  two  wires 
of  the  X-probe  were  made  as  identical  as  one  can  expect  them 
to  be.  The  photograph  proves  that  they  appear  identical  and 
symmetric  to  the  naked  eye  even  with  the  enlargement  of  the 
pictures.  In  addition,  the  two  wires  were  made  from  the 
same  stock  using  the  identical  technique.  After  careful 
examination  under  the  microscope  and  while  using  an  optical 
comparator,  it  was  found  that  Wire  2,  which  exhibited  the 
larger  deviation  from  the  cosine  law,  had  a slight  double 
curvature  along  its  length. 

The  above  results  lead  us  to  believe  that  not  only 
each  probe  but  also  each  new  sensor  mounted  on  the  probe 
must  be  calibrated  even  if  the  new  sensor  is  identical 
to  the  one  being  replaced. 

An  interesting  correlation  between  b*  and  m*  is  shown 
in  Figure  75  where  b*  varies  approximately  as  (2m*)  If 

this  correlation  is  used,  the  present  yaw  relation  can  be 
reduced  to  a one  parameter  relation.  This  would  still  give 
a more  accurate  description  of  the  data  than  the  Friehe  and 
Schwarz  relation  but  will  no  longer  give  the  minimum  errors 
which  can  be  achieved  with  the  two  parameter  relation. 

Again, data  from  MSU  and  from  hot-films  in  water  are  included 
in  this  figure.  The  small  scatter  about  the  correlation  is 
surprising  in  view  of  the  variety  of  the  data  included. 

Considering  the  added  complexity  introduced  by  a two 
parameter  relation  over  the  various  one  coefficient  rela- 
tions, it  is  important  to  find  out  when  a more  complex 


relation  is  most  needed.  The  most  obvious  comparison  is  be- 
tween the  present  and  the  Friehe  and  Schwarz  relations . 

Taking  into  account  the  effect  of  f/d  and  the  Reynolds  num- 
ber, the  regions  where  a substantial  improvement  in  the 
accuracy  is  gained  by  using  the  present  relation  are  shown 
in  Figure  76.  "Significant"  improvement  is  indicated  where 
the  ratio  ot/a*  is  larger  than  two  as  described  in  Chap- 
ter  V.  This  figure  can  then  be  used  as  an  operating  chart 
to  determine  which  of  the  two  yaw  relations  is  more  "suit- 
able" for  a particular  t/d. 

Considering  5%  deviations  from  unity  in  the  value  of 
the  correction  factor  of  the  turbulence  correlations  to  be 
small,  sensors  which  behave  like  the  cosine  law  can  be  iden- 
tified (Note  that  unity  corresponds  to  the  cosine  relation) . 
The  cross  hatched  regions  of  Figure  77  encompass  all  the 
conditions  found  to  fall  under  this  definition.  The  results 
are  somewhat  surprising  in  that  the  short  sensors  should  have 
a non-uniform  axial  temperature  distribution  which  some  in- 
vestigators suggest  causes  deviation  from  the  cosine  law. 
However  the  reason  that  the  longer  sensors  show  some  devia- 
tion from  the  cosine  law  may  be  explained  by  the  method  used 
to  change  t/d.  In  all  t/d  cases,  the  spacing  between  the 
hot-wire  prongs  is  kept  constant  so  that  as  the  sensor  t/d 
increases^  the  active  portion  of  the  sensor  will  be  closer 
to  the  prongs  thereby  increasing  the  conduction  and  aero- 
dynamic interference  effects.  These  increased  effects  are 
probably  only  significant  for  4/d*1000. 
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I If  one  uses  a probe  similar  to  that  used  in  this  inves- 

1 

tigation.  Figures  45  and  76  can  be  utilized  as  operating 
charts  to  determine  which  sensor  would  be  most  applicable. 
From  these,  one  observes  that  going  to  a very  short  sensor, 
e.g.,  £/d=100,  would  be  most  advantageous.  With  this  sen- 
' sor  , the  effects  of  gradients  in  mean  velocity  on  the  wire 

I sensitivity  as  described  by  Gessner  and  Moller  [38]  would 

! be  minimized.  The  short  sensor  would  also  be  more  suitable 

I 

j for  measurement  of  turbulence  quantities,  as  described  by 

Frenkiel  [391 . 

In  Chapter  V,  the  idea  of  finding  the  best  values  of 
the  yaw  relation  parameters  over  a range  of  Reynolds  num- 
bers, i.e.,  velocities,  is  introduced.  Such  values  would  be 
* useful  when  X-probes  are  utilized  to  obtain  the  various  pro- 

I files  of  the  velocity  components  and  their  correlations. 

This  is  particularly  true  in  connection  with  analog  signal 
processing.  The  selection  of  the  best  values  of  the  para- 
meters is  based  on  an  effective  standard  deviation  of  the 
data  from  the  yaw  relation,  over  the  desired  range  of 
Reynolds  numbers.  As  demonstrated  in  Figures  51  and  52, 
employing  the  various  yaw  relations  will  lead  to  different 
effective  standard  deviations.  The  effective  overall  ratio 
of  improvement  gained  by  changing  from  one  relation  to 

I another  can  be  calculated  for  the  range  of  Reynolds  numbers 

! 

I under  investigation  by  dividing  the  values  of  the  best 

effective  standard  deviations  found  from  the  different  yaw 
relations.  This  improvement  ratio  is  given  in  Figure  78  for 


i 


a comparison  between  the  present  relation  and  both  the  co- 


sine law  and  the  Friehe  and  Schwarz  relation.  A negligible 


improvement  is  achieved  by  using  the  present  instead  of  the 


Friehe  and  Schwarz  relation.  A small  improvement  is  made 


by  substituting  the  present  relation  for  the  cosine  law. 


However,  these  results  must  be  interpreted  with  some 


care.  The  magnitude  of  is  typically  an  order  of  magni- 


tude larger  than  that  of  o*  , which  indicates  the  large 

hr 


errors  introduced  by  trying  to  represent  a range  of  Reynolds 


numbers  by  a single  set  of  yaw  relation  parameters.  The 


present  exercise,  therefore,  indicates  that  one  should  use 


different  yaw  calibration  coefficients  for  different  Rey- 


nolds numbers,  or  velocities.  Based  on  the  results  of  the 


exercise  one  can  also  surmise  that  if  a single  set  of  para- 


meters is  to  be  used,  the  accuracy  achieved  by  a simple 


relation  like  that  of  Friehe  and  Schwarz  is  as  good  as  one 


may  expect.  In  fact  for  most  sensor  f/d's, using  the  cosine 


law  is  just  as  good.  Therefore,  it  can  be  stated  that  for 


analog  signal  processing  one  may  just  as  well  use  the  sim- 


plest yaw  relation,  i.e.,  the  cosine  law.  However,  after 


the  data  has  been  gathered,  the  results  must  be  corrected 


for  the  large  errors  introduced  by  the  cosine  law.  These 


errors  are  present  in  all  of  the  mean  and  rms  velocity  com- 


ponents measured  by  the  X-probe  and,  of  course,  in  their 


correlations.  The  correction  factors  must  be  obtained  with 


the  aid  of  the  yaw  calibration  data  of  the  wires  and  most 


likely  they  are  functions  of  the  velocity,  i.e.,  Reynolds 
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number . 

Such  correction  factors  have  been  derived  by  Champagne 
et  a^.  [25]  based  on  his  yaw  relation  for  the  case  when 
the  mean  flow  is  aligned  with  the  axis  of  symmetry  of  the 
X-wires.  They  also  demonstrate  that  a simple  summing  and 
differencing  technique  can  be  used  to  determine  the  various 
turbulence  intensities  and  Reynolds  stresses.  Bruun[40] 
presents  an  approximate  method  which  he  argues  can  be  used 
to  determine  these  quantities  when  the  mean  flow  angle  is 
not  aligned  with  the  probe.  His  approximate  method  again 
relies  on  simple  sum  and  difference  techniques.  However, 
based  on  the  present  analysis,  his  approximation  will  lead 
to  large  errors  in  some  cases.  In  Chapter  III,  it  was 
shown  that  to  first  order,  an  X-probe  can  be  used  with  the 
same  degree  of  accuracy  in  flows  where  the  mean  flow  angle 
is  not  aligned  with  the  probe  as  when  the  mean  flow  angle 
is  zero.  The  approach  of  Chapter  III  is  based  on  simple 
weighted  summing  and  differencing  of  the  signals  from  the 
sensors.  The  weighting  functions  depend  on  the  yaw  coef- 
ficients and  the  mean  flow  angle.  Having  shown  in  Chapter 
VI  that  dynamic  calibration,  as  for  example  used  by  Perry 
et  a^.  [41]  , or  that  calibration  in  a turbulent  flow  may  be 
necessary,  one  ascertains  that  the  method  proposed  in  Chap- 
ter III  is  much  more  suitable.  In  particular,  since  it 
has  no  restrictions  on  the  type  of  flowfield  it  is  cali- 
brated in.  This  method  will  also  express  the  accuracy  of 
the  relation  representing  the  calibration  data  by  the  , 
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quantity  a^. 

To  demonstrate  the  importance  of  knowing  the  standard 
deviation  of  the  yaw  data,  and  how  even  small  differ- 

ences in  can  cause  substantial  discrepancies  in  the  com- 
puted turbulence  intensities  for  different  mean  flow  angles, 
Figures  79  through  81  are  presented  . In  order  to 
determine  these  discrepancies,  the  distribution  of  energy  be- 
tween the  different  components  of  the  turbulence,  as  calcu- 
lated by  the  cosine  law  must  be  known.  Sample  results  are 
shown  for  two  different  flow  conditions,  one  typical  of  the 
outer  region  of  a turbulent  boundary  layer,  the  other  typi- 
cal of  the  inner  region.  Data  for  these  flow  conditions 
were  based  on  the  measurements  of  Klebanoff  [42]  . While 
only  very  small  mean  flow  angles  with  respect  to  the  X-probe, 
y,  are  present  in  the  boundary  layer,  Klebanoff 's  data  is 
used  here  for  0“<  y £ 25° to  illustrate  the  trends  and  to 
give  the  reader  a feeling  for  the  magnitudes  of  the  errors 
involved . 

The  yaw  parameters  used  in  Figures  79  and  80  are  for 
the  case  where  the  largest  deviations  were  observed  between 
the  three  yaw  relations,  i.e.,  the  worst  case  expected  based 
on  the  probes  tested  here. 

In  Figures  79,  80  and  81  the  differences  between  values 
computed  by  the  cosine  law  and  those  computed  by  either  the 
present  relation  or  the  Friehe  and  Schwarz  relation  are  re- 
presented by  and  respectively.  The  subscript  i 

denotes  the  turbulence  velocity  correlation  computed.  For 
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example,  an  e of  0.1  indicates  a 10%  difference  between  the 
values  computed  by  the  two  relations  in  question  based  on 
the  same  output  signals  from  the  X-wire  and  the  same  yaw 
calibration  data.  The  differences  between  the  present 
and  the  Friehe  and  Schwarz  relation  can  be  obtained  by  the 
difference  between  and  (e^);  see  Equation  (III-73)  . 

In  Figure  79,  when  the  mean  flow  angle,  y,  is  zero,  a 
4.5%  difference  in  computed  values  of  v*  exist  between  the 
Friehe  and  Schwarz  relation  and  the  present  relation  even 
though  only  less  than  0.5%  difference  in  fitting  the  mean 
velocity  is  observed  in  Figure  §0.  This  difference  in  v' 
increases  to  10%  as  y increases.  The  major  effect  appears 
in  the  Reynolds  stress  uv  where  a difference  of  over  70% 
exists  between  the  Friehe  and  Schwarz  and  the  present  rela- 
tions for  increasing  values  of  y,  while  the  difference  in 
uv  between  the  cosine  law  and  Friehe  and  Schwarz  varies 
between  9 and  16%.  When  y = 0,  the  minimum  difference  of 
5.5%  in  the  calculated  uv  is  found  between  the  present  and 
the  Iriehe  and  Schwarz  relations. 

The  distribution  of  the  turbulence  does  have  some  ef- 
fect on  the  differences  in  the  computed  values  of  v' , how- 
ever it  plays  a larger  role  in  the  determination  of  the 
Reynolds  stress.  As  the  magnitude  of  the  v component  of  the 
fluctuations  approaches  that  of  the  u f lucutations,  the 
difference  in  uv  computed  by  the  present  relation  and  Friehe 
and  Schwarz'  relation  substantially  increases.  At  the  same 
time  the  difference  between  the  Friehe  and  Schwarz  relation 


and  the  cosine  law  increases  only  by  approximately  7%  at 
the  large  y's. 
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Differences  in  mean  and  streamwise  fluctuating  velo- 
cities for  the  same  case  appear  in  Figure  80.  Although 
small  differences  in  u*  are  found  when  y - 0,  these  dif- 
ferences grow  to  12%  as  y increases.  For  y = 0 the  differ- 
ence in  the  computed  value  of  u'  is  identical  to  the  dif- 
ference in  mean  velocity.  Therefore  the  ratio  of  these, 
i.e.,  the  streamwise  turbulence  intensity,  remains  the  same 
in  all  three  cases.  This  is  true  in  the  case  of  an  ideal  X- 
probe.  However  once  a typical  X-probe  in  which  each  sen- 
sor has  a slightly  difference  yaw  sensitivity  is  used, 
there  will  be  differences  in  the  computed  turbulence  inten- 
sities. The  results  of  the  analysis  demonstrate  that  the 
longitudinal  turbulence  intensity  is  the  only  quantity  that 
can  be  measured  exactly,  even  then  only  with  an  ideal  probe 
with  matched  sensors  at  equal  angles  to  the  mean  flow  dir- 
ection. It  is  ironic  that  this  quantity  can  be  measured 
just  as  accurately  by  the  much  simpler  method  of  one  hot- 
wire perpendicular  to  the  flow  direction. 

It  should  be  noted  that  in  all  cases,  when  Y = 0,  the 
differences  between  the  computed  values  reduce  to  those 
given  by  Equations  (III-45)  through  (III-47)  and  are  in 
agreement  with  the  results  of  Friehe  and  Schwarz  [22]  and 
of  Roberts  et  al.  [23,  24] . While  the  present  analysis  is 
limited  to  small  turbulence  intensities,  by  extending  it  to 
the  cases  y 0,  one  is  able  to  get  an  indication  of  the 
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trends  at  high  intensities.  Under  large  turbulence  condi- 
tions the  instantaneous  value  of  the  angle,! .e. , 3 ,can  achieve 
very  large  values.  If  the  heat  transfer  from  the  wire  be- 
haves in  a somewhat  quasi-steady  fashion,  the  y 0 results 
may  give  good  indications  of  the  errors  introduced  by  the 
simpler  yaw  relations. 

The  case  for  an  X-probe  constructed  of  standard-length 
hot-wires  corresponding  to  the  region  of  negligible  im- 
provement between  the  three  relations,  is  shown  in  Figure  81 
for  the  scune  flow  conditions  as  those  of  Figures  79  and  80. 
The  differences  in  the  computed  values  of  u'  and  v'  are 
negligible  between  the  Friehe  and  Schwarz  and  the  present 
relations  while  the  differences  in  uv  are  only  as  large 
as  4%.  The  typical  differences  for  such  an  ideal  probe  are 
of  the  order  of  1%.  The  difference  between  for  Friehe 
and  Schwarz  and  the  present  data  in  this  case  is  0.01%. 

The  results  of  Figure  81  demonstrate  that  large  dif- 
ferences in  the  computed  turbulence  velocity  correlations 
can  exist  even  for  insignificant  differences  (i.e.,  0.01%) 
in  fitting  the  mean  yaw  calibration  data.  In  addition,  the 
magnitude  of  the  standard  deviation  of  the  yaw  calibration 
fit  must  be  very  small.  The  data  of  Figures  79  and  80  show 
that  a value  of  of  approximately  0.7%  is  certainly  not 
small  enough,  a most  surprising  conclusion.  Reducing  this 
value  by  a factor  of  about  4.8,  i.e.,  to  a ; 0.15%,  leads 
to  very  significant  reductions  in  the  errors  in  measuring 


the  various  turbulence  quantitites.  The  errors  in  the  most 
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critical  of  these  quantities,  namely  uv,  are  reduced  by  6 to 
70%  depending  on  the  inclination  of  the  X-wires  with  respect 
to  the  mean  flow  direction  (or  depending  on  the  intensity 
of  the  turbulence) . Therefore,  a very  high  degree  of  accu- 
racy in  fitting  the  mean  yaw  calibration  data  is  needed  to 
attain  accurate  results. 

The  sample  results  presented  in  Figures  79,  80  and  81 
utilized  information  on  the  distribution  of  the  turbulence 
energy,  i.e.,  and  , based  on  the  classical  experi- 

ments of  Klebanoff  [42].  In  practice,  the  correction  fac- 
tors must  be  derived  based  on  the  data  obtained  from  the 
experiment  in  addition  to  the  yaw  calibration  parameters. 
Such  corrections  can  be  made  a posteriori  if  analog  signal 
processing  is  used  or  can  be  incorporated  in  a subroutine 
to  the  digital  data  processing  program.  In  the  former  case 
the  analog  system  would  be  set  based  on  the  cosine  relation 
because  of  the  discussion  on  the  accuracies  achieved  over  a 
Reynolds  number  range  presented  earlier  in  this  section.  If 
the  Reynolds  nximber  of  the  probe  remains  unchanged  during 
the  experiment,  the  corrections  can  be  included  a priori  to 
the  experiments.  They  would  be  automatically  incorporated 
in  the  weighted  sums  and  differences  of  the  output  signal 
based  on  the  yaw  calibration  parameters  corresponding  to  the 
operating  Reynolds  number.  The  details  of  this  approach  are 
given  in  Chapter  III. 

In  regard  to  the  effects  of  velocity  fluctuations  on 
the  yaw  calibration  of  hot  wires,  the  results  obtained  in 
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this  study  are  presented  and  discussed  in  Chapter  VI . How- 
ever the  main  thrust  of  the  findings  will  now  be 
summarized . 

Various  sensors  from  single-  and  X-probes  were  cali- 
brated under  steady  conditions  as  well  as  in  the  presence 
of  velocity  fluctuations.  The  results  indicate  that  while 
dynamic  yaw  calibration  of  the  probes  is  not  necessary,  the 
mean  yaw  sensitivity  of  the  probe  should  be  obtained  in  the 
presence  of  some  velocity  fluctuations,  e.g.,  in  the  pres- 
ence of  some  background  turbulence.  These  results  also 
suggest  that  the  velocity  calibration  of  the  probe  may  fol- 
low the  same  trends.  This  of  course  is  contradictory  to  the 
common  practice  of  using  calibration  flow  conditions  with 
the  minimum  possible  turbulence  intensity.  However,  they 
are  not  necessarily  surprising  in  view  of  the  changes  which 
may  be  induced  by  unsteady  velocities  in  the  boundary  lay- 
ers ai.d  in  the  flow  around  the  sensor  and  prongs  as  well  as 
in  other  aspects  of  the  convective  heat  transfer  from  the 
wire . 
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CHAPTER  VIII 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  conclusions  drawn  from  the  various  parts  of  the  in- 
vestigation are  presented  in  the  following  order:  first, 
conclusions  obtained  from  the  analysis  of  the  temperature 
compensating  scheme,  excunined  in  Chapter  II  and  Appendix  A, 
are  presented.  Next,  the  main  findings  of  the 
experiments  on  temperature  compensation,  reported  in  Appen- 
dix B,are  summarized.  The  major  conclusions  drawn  from  the 
analysis  describing  the  yaw  response  of  hot-wires  in  steady 
and  fluctuating  velocity  fields  are  presented  next.  This 
analysis  is  given  in  detail  in  Chapter  III.  Finally,  con- 
clusions based  on  the  results  of  the  experiments  on  the  yaw 
sensitivity  of  hot-wires  are  listed.  The  discussion  of 
these  results  is  included  in  Chapters  VII,  VI  and  V. 

Conclusions 

Based  on  the  analysis  of  the  scheme  proposed  for  com- 
pensation of  the  anemometer  output  for  ambient  temperature 
variations,  it  was  found  that: 

1.  Temperature  compensation  does  not  depend  on 
whether  the  anemometer  output  is  linearized  or  not. 

2.  One  needs  only  to  know  the  properties  of  the  ane- 
mometer bridge  and  the  velocity  and  temperature  probes 
(i.e.,  their  resistance  and  temperature  dependence  co- 
efficients) in  order  to  achieve  the  compensation,  without 
any  need  for  temperature  calibration.  Hence,  the  scheme 
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provides  means  for  incorporating  the  temperature  compensa- 
tion a priori  to  conducting  the  experiments. 

3.  Constraints  on  the  parameters  of  the  temperature 
compensating  probe,  valid  for  all  three  circuit  configura- 
tions, were  derived.  With  this  and  the  knowledge  of  the 
properties  of  available  velocity  sensors,  operating  over- 
heat ratios,  anemometer  bridge  ratio  and  operating  mode,  a 
suitable  temperature  compensating  probe  may  either  be  fab- 
ricated or  purchased. 

4 . Exact  solutions  for  temperature  compensating  cir- 
cuits 1 and  2 were  obtained  for  three  different  operating 
modes;  i.e.,  constant  overheat  ratio,  constant  overheat  dif- 
ference and  constant  output  voltage.  Operating  charts  are 
provided  (Figures  4 and  5)  for  the  selection  and  adjustment 
of  the  resistances  in  the  circuits. 

5.  When  no  constraints  are  imposed,  a family  of  solu- 
tions exist  for  compensating  circuit  3.  Constraints  to 
minimize  the  non-ideal  behavior  are  being  derived  for  this 
circuit  making  it  the  optimum  circuit  configuration. 

6.  Estimates  for  the  errors  introduced  through  the 
non-ideal  behavior  of  the  compensating  circuits  have  been 
derived  and  indicate  that  the  error  is  proportional  to  the 
square  of  the  temperature  variation. 

7.  For  a typical  hot-wire,  temperature  compensator 
pair,  it  was  determined  that  circuit  1 with  an  ambient  tem- 
perature increase  of  40"C  produces  a -2%  change  in  the  in- 
dicated velocity,  while  utilizing  circuit  2 reduces  the 
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error  to  -1%.  This  error  can  be  further  reduced  if  the  op- 
timized configuration  of  circuit  3 is  employed. 

8.  The  temperature  dependence  of  the  coefficient  A 
in  King's  Law,  which  was  found  from  the  experiments,  was 
ignored  in  the  analysis  and  error  estimates.  Its  variation 
with  temperature  will  add  a small  contribution  to  the  errors 
in  the  compensation  attained.  In  order  to  account  for  the 
temperature  dependence  of  this  or  the  other  coefficients 
in  King's  Law,  i.e.,  n and  B,  the  a priori  setting  of  the 
scheme  and  estimating  of  the  errors  will  not  be  possible 
without  velocity  calibration  at  different  temperatures. 

Based  on  the  experiments  on  temperature  compensation 
conducted  in  air,  the  following  conclusions  were  drawn: 

1.  For  a constant  velocity,  if  the  ambient  temperature 
changes,  operating  in  a constant  overheat  mode  overcompen- 
sates but  seems  to  approach  the  ideal  compensation  to  a 
greater  degree  of  accuracy  as  compared  to  the  constant  dif- 
ference mode  of  operation,  which  undercompensates. 

2.  The  effective  temperature  compensation  of  all  oper- 
ating modes  was  found  to  be  independent  of  fluid  velocity 
and  initial  overheat  ratio. 

In  the  derivation  of  the  response  equations  of  inclined 
cylindrical  sensors  in  steady  and  fluctuating  velocity 
fields,  previous  analyses  have  been  extended  for  the  pres- 
ent yaw  relation  to  include  X-probes  with  dissimilar  senscrs 
as  well  as  mean  velocities  at  an  angle  y to  the  line  of 
symmetry  of  the  X-wires.  The  following  conclusions  are 
based  on  the  results  of  the  analysis. 
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1.  For  sensors  with  yaw  dependence  that  does  not 
closely  follow  the  cosine  law,  neither  the  simple  nor  an 
approximate  weighted  sum  and  difference  technique  applies 

when  yf^O . This  makes  analog  processing  to  obtain  accurate  mean 
velocities  extremely  difficult  if  not  impossible.  However, 
for  cases  where  the  mean  streamline  is  aligned  with  the  probe 
axis  of  symmetry,  the  errors  introduced  are  usually  less 
than  5%. 

2.  In  low-level  turbulence,  time  dependent  velocity 
components  can  be  obtained  using  a weighted  sum  and  dif- 
ference technique  once  the  mean  flow  angle  is  known. 

3.  By  extending  the  analysis  to  cases  where  y ^ 0, 
while  using  low  intensity  asstunptions , one  is  able  to  get  an 
indication  of  the  trends  at  high  intensities  for  which  the 
instantaneous  flow  angle  can  assume  very  large  values. 

4 . The  correction  factors  derived  for  the  measured  mean 
and  turbulence  quantities  must  be  based  on  the  data  obtained 
from  the  experiment  in  addition  to  the  yaw  calibration  para- 
meters. Such  corrections  can  be  made  a posteriori  if  analog 
signal  processing  is  used  or  can  be  incorporated  in  a sub- 
routine to  the  digital  data  processing  program. 

5.  Discrepancies  exist  between  different  yaw  relations 
in  calculated  values  of  turbulence  velocity  correlations 
which  are  derived  from  the  same  data.  These  discrepancies 
depend  on  the  mean  flow  angle,  the  distribution  of  the  tur- 
bulence energy  between  the  velocity  components  and  the  ac- 
curacy to  which  the  mean  flow  calibration  data  fits  the  par- 
ticular yaw  function.  For  differences  in  this  accuracy  of 
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0.5%,  discrepancies  as  large  as  10%  and  70%  may  exist  in  the 
computation  of  v'  and  uv,  respectively.  Differences  in  fit- 
ting mean  calibration  data  of  0.01%  produce  discrepancies  in 
turbulence  quantities  of  the  order  of  1%  when  y * 0 and  as 
large  as  4%,  in  the  Reynolds  stress,  at  large  y's. 

6.  The  longitudinal  turbulence  intensity  is  the  only 
quantity  that  can  be  measured  exactly,  even  then  only  with 
an  ideal  probe  having  matched  sensors  at  equal  angles  to  the 
mean  flow  direction. 

7.  Measurements  in  many  complex  flowfields  such  as 
"swirling  flows"  or  near  bluff  bodies,  in  particular  in  the 
presence  of  turbulence,  require  that  the  yaw  sensitivity  of 
both  sensors  of  an  X-probe  must  be  known  to  a high  degree  of 
accuracy . 

Conclusions  based  on  the  experiments  on  yaw  sensitivity 
in  steady  and  fluctuating  velocities  are  listed  in  the 
following: 

1.  Unlike  the  cosine  law  or  the  Priehe  and  Schwarz  re- 
lation, the  standard  deviation  of  fitting  the  yaw  calibra- 
tion data  by  the  present  relation  is  independent  of  the 
Reynolds  number  and  the  length  of  the  sensor,  i.e.,  l/d. 

2.  In  all  cases  examined,  including  hot-films  in  water 
and  data  from  an  independent  source,  the  present  relation 
fit  the  yaw  calibration  data  to  a higher  degree  of  accuracy 
than  the  cosine  law  or  the  Friehe  and  Schwarz  relation.  The 
improvement  in  the  accuracy  over  the  cosine  law  can  be  as 
large  as  a factor  of  25,  while  that  over  the  Friehe  and 
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Schwarz  relation  can  only  be  as  large  as  5.  This  improve- 
ment, however,  is  dependent  on  probe  design  and  Reynolds 
number . 

3.  The  ratio  of  improvement  over  either  the  cosine  law 
or  the  Friche  and  Schwarz  relation  is  smaller  for  the  very 
short  sensors  (f/d  ~ 100)  than  for  the  long  sensor  (£/d  = 
1000) . This  indicates  that  the  cosine  law  is  more  appli- 
cable to  the  short  sensors,  rather  than  to  the  long  ones! 

The  variation  of  the  yaw  parameters  was  not  found  to  be 
monotonic  with  sensor  f/d  and  did  agree  with  previously 
reported  results. 

4.  The  improvement  in  using  the  present  relation  over 
the  Friehe  and  Schwarz  relation  comes  at  relatively  large 
yaw  angles  (i.e.,  larger  than  50“)  and  particularly  at  low 
Reynolds  niombers.  The  yaw  coefficients  determined  from  the 
present  relation  are  independent  of  the  angular  range  of 
data  between  ±50“  and  +70“. 

5.  A good  yaw  relation  must  account  for  both  higher 
. d lower  effective  cooling  rates  than  that  determined  by 
the  cosine  law. 

6.  The  yaw  sensitivity  of  hot  wires  is  not  significant- 
ly Influenced  by  the  operating  overheat  ratio  in  the  range 
between  1.4  and  1.8. 

7.  The  trend  of  the  yaw  parameters  with  either  sensor 
l/d  or  Reynolds  number  does  not  seem  to  approach  the  cosine 
law,  but  rather  it  seems  to  overshoot  it  without  any  observ- 
ed asymptotic  values. 
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8.  In  general,  if  an  X-probe  is  misaligned  by  0.5®, 
the  effect  on  the  yaw  coefficients  is  equivalent  to  that 
resulting  from  recalibrating  the  probe,  properly  centered. 


at  a freestream  velocity  3 ft/sec  lower  than  the  desired 
one . 

9.  An  approximate  correlation  can  be  used  to  reduce 
the  present  two  parameter  yaw  function  to  a single  para- 
meter relation  and  is  given  by  b*  = l/2m*.  Employing  this 
relation  would  still  yield  a more  accurate  description  of 
the  yaw  calibration  data  but  will  no  longer  give  the  mini- 
mum possible  errors. 

10.  For  all  cases  examined,  values  of  the  transverse 
sensitivity  coefficient  calculated  from  a yaw  sensitivity 
function,  taking  into  account  deviations  from  the  cosine 
law,  yield  a higher  degree  of  accuracy  than  those  obtained 
directly  from  the  data. 

11.  All  sensors  require  a turbulence  correction  factor 
to  take  into  account  deviations  from  the  cosine  law.  How- 
ever, by  proper  design  and  selection,  the  magnitude  of  this 
correction  factor  can  be  minimized. 

12.  One  should  use  different  yaw  calibration  coef- 
ficients for  different  velocities.  If  a single  set  of  para- 
meters is  to  be  used,  the  accuracy  achieved  by  a simple  re- 
lation like  the  cosine  law  is  as  good  as  one  may  expect. 
However,  after  the  data  has  been  gathered,  the  results  must 
be  corrected  for  the  large  errors  introduced  into  the  mean 
and  rms  velocity  components. 
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13.  The  yaw  parameters  are  influenced  by  the  probe 
design  and  construction.  Therefore,  each  probe  must  be  cali- 
brated to  determine  the  yaw  parameters  associated  with  it  at 
the  Reynolds  numbers  of  interest.  The  results  lead  us  to 
believe  that  not  only  each  probe  but  also  each  new  sensor 
mounted  on  the  probe,  must  be  calibrated,  even  if  the  new 
sensor  is  identical  to  the  one  being  replaced. 

14.  The  parameters  of  the  yaw  relation  are  affected  by 
the  presence  of  unsteady  velocities  in  the  flowfield  and 
depend  on  the  frequency  of  these  fluctuations. 

15.  While  dynamic  yaw  calibration  of  the  probes  is  not 
necessary,  the  mean  yaw  sensitivity  of  the  probe  should  be 

obtained  in  the  presence  of  some  velocity  fluctuations, 
e.g.,  in  the  presence  of  some  background  turbulence  of  simi- 
lar spectral  content  to  the  flow  under  consideration. 

Recommendations 

The  following  items  are  recommended  for  future  investi- 
gations in  connection  with  the  proposed  temperature  compen- 
sation scheme: 

1.  The  optimization  analysis  of  compensating  circuit 
3 should  be  completed. 

2.  The  analysis  of  the  compensating  scheme  should  be 
extended  to  incorporate  the  variations  of  the  coefficients 
A and  B in  King's  Law  with  ambient  temperature;  this  should 
include  estimates  of  the  errors  involved.  The  variation  of 
A with  temperature  is  important  in  air  while  the  variation 
of  B with  temperature  will  probably  be  important  in  water. 
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3.  The  errors  can  be  calculated  exactly  from  the  ex- 
perimental data  presented  in  Appendix  B.  They  should  be 
found  by  comparing  the  values  based  on  the  experimental  data 
of  dR3/dT,  i.e.,  the  values  required  for  a constant  anemom- 
eter output  at  a fixed  velocity  with  varying  temperature, 
with  those  derived  by  the  analysis  of  Chapter  II. 

4.  Experiments  using  real  compensating  probes  should 
be  performed  to  verify  the  analysis  of  Chapter  II. 

The  following  items  are  recommended  to  develop  further 
insight  into  the  yaw  dependence  of  hot-wires. 

1.  The  analysis  can  be  extended  to  account  for  the 
effect  of  small  lateral  velocity  fluctuations,  i.e.,  w,  on 
the  various  turbulence  correlations. 

2.  The  analysis  can  also  be  extended  to  determine  the 
variation  of  the  measured  u'  and  v'  as  functions  of  y. 

3.  Experiments  on  yaw  calibration  in  fluctuating  velo- 
city fields  should  be  modified  to  also  include  controlled 
periodic  velocity  fluctuations  in  the  transverse  direction. 
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INITIAL  ERROR  ESTIMATES  FOR  TEMPERATURE 
COMPENSATING  CIRCUITS 

Instead  of  the  approach  presented  in  Chapter  II  for 
determining  the  errors  introduced  when  using  the  tempera- 
ture compensating  circuits,  the  following  approximate  meth- 
od was  initially  used.  If  one  assumes  that  over  a reason- 
able range  of  ambient  temperature  variations  (e.g.,  for  a 
temperature  change  of  40“C  in  air)  the  equivalent  resist- 
ance Rj  varies  linearly,  the  operating  velocity  probe  re- 
sistance Rjj  will  also  be  linear  due  to  the  bridge  balance 
condition  . The  resistance  R^  can  then  be  written  as: 

■Si  > "Sio  “h  - ’’o*  1 

Operating  the  anemometer  and  compensating  circuit  in 
the  constant  bridge  output  voltage  mode,  we  assume  that  the 
temperature  function  t (T) , which  is  defined  in  Chapter  II, 
can  be  approximated  by 

t(T)  • t(T^)  - m(T  - T^)  (A-2) 

where  t (T  ) is  a constant  evaluated  at  the  reference  tem- 
o 

perature  and  m(T  - T^)  is  a small  variation  away  from  this 
constant  value.  Here,  the  coefficient  m can  be  deter- 
mined once  R^(T)  and  R^CT)  are  known. 

Using  Equation  (A-2)  with  Equation  (11-21)  one  ob- 
serves that  3Ej^/3T  ^ 0.  The  resulting  variation  of  bridge 
output  voltage  with  euobient  temperature  can  be  given  as 
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(A-3) 


If  the  anemometer  output  voltage  is  linearized,  fol- 
lowing a similar  procedure  one  finds  that 
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1 + ag(T  - T^)_ 

(A-4) 


With  available  information  on  the  temperature  depend- 
ence of  the  compensating  probe  and  the  velocity  probe. 
Equation  (A-3)  or  (A-4)  may  be  used  to  estimate  the  drift 
in  output  voltage.  The  accuracy  of  this  method  solely  de- 
pends upon  the  accuracy  of  Equation  (A-1)  which  can  be 
experimentally  determined. 
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EXPERIMENTS  ON  TEMPERATURE 
COMPENSATION 

In  parallel  with  the  analysis  of  the  temperature  compen- 
sating circuits  presented  in  Chapter  II  and  Appendix  A,  ex- 
periments were  conducted  in  air  to  determine  the  validity 
of  several  of  the  assumptions  made  and  to  determine  the 
temperature  dependence  of  some  of  the  parameters  derived. 

Measurements  were  carried  out  utilizing  the  I.I.T.  com- 
pressed-air driven  hot-wire  calibration  tunnel.  Dry  air  was 
heated  by  a 4.5  kw  "Chromalox"  Model  GCH3405  heater  prior  to 
entering  the  plenum  chamber  of  the  calibration  tunnel.  The 
air  temperature  was  varied  by  regulating  the  current  to  the 
heater  via  a 3-phase  voltage  control  unit.  Air  temperature 
was  monitored  in  the  test  section  of  the  calibration  tunnel 
by  a specially  constructed  circuit  utilizing  a thermistor. 

A digital  voltmeter  indicated  the  temperature  in  °C  directly 
and  the  accuracy  and  repeatability  of  the  readings  were 
within  0.5°C  over  the  range  from  20  to  60°C.  At  high  velo- 
cities a thermistor  near  the  upstream  end  of  the  test  sec- 
tion was  used,  while  at  low  velocities,  the  temperature  was 
monitored  approximately  2 cm.  downstre£un  of  the  hot-wire  to 
reduce  errors  due  to  heat  loss  to  the  surroundings.  The 
flow  rate  was  controlled  by  a Watts  Model  110  pressure  regu- 
lator in  line  with  and  ahead  of  the  air  heater.  A schematic 
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of  the  entire  setup  is  shown  in  Figure  82. 

The  temperature  coefficient  of  resistance  of  two  compen- 
sating probes  and  a standard  length  hot-wire  (i/ <1=580)  were 
obtained  by  exposing  the  probes  to  heated  air  streams  of 
known  temperature.  Probe  resistance  was  measured  using  the 
bridge  of  a Disa  55D01  anemometer  unit  to  within  +0.005  fl. 

The  temperature  was  varied  and  measurements  taken  through- 
out the  range  from  20  to  60*'C.  Graphs  of  probe  resistance 
versus  temperature  were  plotted  and  the  temperature  coeffi- 
cients of  resistance  determined  from  the  slope  of  the  curve. 
The  temperature  dependence  of  all  probe  resistances  was 
found  to  be  linear  for  the  range  of  temperature  variation. 

Guided  by  the  analysis#  three  different  operating  modes 
were  examined.  The  first  kept  the  operating  overheat  ratio 
constant  by  manually  adjusting  the  decade  resistance  on  the 
anemometer  at  each  temperature  where  data  was  taken.  The 
second  mode  kept  the  operating  resistance  difference, 

constant  by  manually  adjusting  the  decade  resistance 
at  each  temperature.  The  third  mode  of  operation  involved 
adjusting  the  decade  resistance  on  the  anemometer  to  achieve 
constant  output  voltage  at  each  fixed  velocity  as  the  air 
temperature  was  varied.  The  operating  resistance  as  a func- 
tion of  ambient  temperature,  required  to  maintain  the 
bridge  output  constant,  was  found  from  the  latter  approach. 

Results  from  this  portion  of  the  experiment  are  shown 
in  Figures  83  through  85.  Figure  83  demonstrates  that  with- 
out compensation,  the  non-dimensional  bridge  voltage 
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decreases  linearly  with  increasing  temperature  while  main- 
taining the  velocity  constant.  This  result  was  found  to  be 
independent  of  flow  velocity  from  3.5  to  30  m/s.  This  drop 
in  the  bridge  output  voltage  is  dependent  on  the  overheat 
ratio  used  and  is  largest  for  low  overheats.  Three  over- 
heats were  used  : 1.4/  1.6  and  1.8,  the  latter 

being  the  most  common  value  employed  by  various  ex- 
perimenters. 

Holding  the  resistance  difference  constant,  again  while 
maintaining  a constant  velocity,  the  drop  in  output  voltage 
was  found  to  decrease.  Ideally,'  the  bridge  output  voltage 
for  a velocity  sensor  should  remain  constant  for  a constant 
velocity  if  the  ambient  temperature  changes.  The  constant 
overheat  mode  over-compensates  but  seems  to  approach  the 
ideal  compensation  to  a greater  degree  of  accuracy  than  the 
constant  difference  mode.  Both  of  these  operating  modes  are 
found  to  be  independent  of  fluid  velocity  and  initial  over- 
heat ratio. 

If  the  bridge  output  is  forced  to  remain  constant  at  a 
given  velocity,  i.e.,  the  third  mode  of  operation  described 
above,  the  variation  in  the  overheat  ratio  and  resistance 
difference  can  be  determined  as  shown  in  Figure  84.  The 
resistance  difference  has  to  increase  substantially  while 
the  overheat  ratio  must  drop  slightly  as  the  temperature  in- 
creases in  order  to  achieve  ideal  compensation.  This  tends 
to  suggest  that  compensation  using  a constant  overheat  mode 
would  give  better  results  as  compared  to  the  constant 
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difference  mode,  a result  which  is  in  agreement  with  Fig- 
ure 83.  Each  of  these  cases,  normalized  in  the  proper  way, 
is  independent  of  fluid  velocity  and  initial  overheat  ratio. 

The  variation  of  the  function  t with  ambient  tempera- 
ture is  shown  in  Figure  85.  Use  of  this  figure  permits  the 
evaluation  of  the  coefficient  m,  discussed  in  Appendix  A, 
which  is  needed  for  the  approximate  estimate  of  the  errors 
of  the  compensating  circuits.  The  linear  variation  of  x 
with  temperature  substantiates  the  assumptions  made  in 
Appendix  A. 

The  hot-wire  used  in  these  measurements  was  calibrated 
for  velocity  over  a range  of  overheat  ratios  and  ambient 
temperatures.  The  summary  of  the  results  is  shown  in  Fig- 
ure 86.  Both  of  the  coefficients  A'  and  B'  are  found  to  be 
dependent  on  the  eunbient  temperature  and  initial  overheat 
ratio.  A*  increases  with  temperature  and  with  overheat 
ratio  while  the  coefficient  B'  decreased  with  overheat 
ratio  and  with  ambient  temperature.  The  exponent  n was 
found  to  be  the  least  sensitive  to  temperature,  varying  only 
within  4%  of  its  mean  value  for  all  cases. 

The  parameters  A and  B in  King's  Law  can  be  extracted 
from  these  results  with  the  aid  of  Equation  (11-19) . The 
calculated  values  are  plotted  in  Figure  87.  The  curves  in- 
dicate that  the  coefficient  A has  a slight  temperature  de- 
pendence while  the  coefficient  B remains  almost  the  same  for 
all  cases. 

The  analysis  of  the  constant  output  voltage  mode  in 
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Chapter  II  was  based  upon  the  assumption  that  the  parameters 
in  King's  law,  i.e.,  A,  B and  n,  are  independent  of  temper- 
ature. Based  on  Figures  83  and  87  one  can  conclude  that  the 
variation  in  B and  n may  be  neglected.  However,  the  tem- 
I perature  dependence  of  A will  introduce  a small  error  in 

! the  resistances  computed  for  proper  compensation.  The  mag- 

nitude of  this  error  can  be  calculated  from  the  experimental 
data . 
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YAW  CALIBRATION  RESULTS  FROM  HOT-FILMS 
IN  WATER 

Initial  work  on  the  present  yaw  relation  was  done  at 
l.I.T.  by  F.  A.  Roberts  [23]  using  hot-films  in  water  in 
an  attempt  to  explain  some  of  the  trends  shown  in  Fig- 
ure 88.  Data  presented  by  Friehe  and  Schwarz  [22]  has  been 
reproduced  in  this  figure  using  their  method  of  presenta- 
tion. The  data  indicates  that  hot-films,  which  have  small 
f/d's,  seem  to  deviate  less  from  the  cosine  law  as  compared 
to  the  hot-wires  which  have  much  larger  f/d's.  Although 
for  each  of  the  films  or  the  wires  increasing  t/d  decreases 
the  deviation  from  the  cosine  relation,  the  behavior  of  the 
data  in  Figure  88  remained  unexplained  after  the  studies  of 
Friehe  and  Schwarz  [22]  and  of  Roberts  [23] . 

The  summary  of  Roberts'  data  is  presented  in  this  ap- 
pendix and  in  Figures  72  through  75  since  his  original 
work  is  not  readily  available  to  the  reader.  His  data,  for 
a hot-film  with  f/d=18,  are  plotted  in  a similar  fashion  to 
those  of  Friehe  and  Schwarz  in  Figure  89.  The  Reynolds 
number  dependence  of  the  best  fit  yaw  coefficients  (in  the 
least  squares  sense)  as  determined  from  the  present  rela- 
tion demonstrates  the  importance  of  studying  yaw  rela- 
tions over  a range  of  Reynolds  numbers.  To  date  we  have 


not  been  able  to  find  many  such  studies  in  the  literature. 

Contours  in  the  b-m  plane  of  constant  standard  de- 
viation of  the  data  from  the  present  relation  are  shown  in 
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Figures  90  through  92  for  the  same  sensor.  It  should  be 
noted  here  that  the  definition  of  a'  differs  from  that  of 
Op  in  Chapter  V.  Here  o'  is  the  standard  deviation  of 
f(a),  namely. 


o ' 


f(a^) 


U 


OQ 


2 


(C-1) 


where  la|  varied  from  0®  to  70*.  Figures  90  and  91  show 
trends  of  constant  o'  similar  to  those  described  in  Chap- 
ter V. 

Figure  92  shows  the  dependence  of  the  minimum  regions 
of  constant  standard  deviation  on  Reynolds  number  for  a 
standard  commercial  hot-film  in  a typical  range  of  velo- 
cities. The  movement  of  these  regions  is  not  consistent 
with  the  trends  observed  in  Chapter  V for  hot-wires  in  air. 
This  may  suggest  that  even  though  both  types  of  velocity 
sensors  obey  the  present  relation,  their  variation  with 
pareuneters  examined  in  the  present  work  may  be  vastly 


different. 
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YAW  CALIBRATION  RESULTS  FROM 
X-WIRE  PROBES;  DATA  COURTESY  OF  MSU 

In  an  attempt  to  verify  the  behavior  of  the  hot  wires 
used  at  I.I.T.  as  typical,  two  independent  sets  of  yaw  cali- 
bration data  were  made  available  by  Dr.  John  Foss  of  Michi- 
gan State  University. 

Both  sets  of  calibration  data  were  obtained  using  an 
X-wire  probe  with  sensors  having  an  £/d=200.  Each  set  was 
obtained  through  different  calibration  procedures.  The 
initial  X-wire  data,  called  Set  1,  consisted  of  a free- 
stream  velocity  calibration  of  each  sensor  at  various  angles 
of  yaw. 

In  his  analysis  of  similar  data,  Foss  [43]  used  the 
Friehe  and  Schwarz  yaw  relation  and  assumed  that  the  velo- 
city calibration  equation  at  any  angle  of  yaw  can  be  written 
as 

eJ  = A'  + B'  u”  (D-1) 

b e 

2 

Since  U^  and  U^  are  related  by  f (a) , Equation  (D-1)  can  be 
rewritten  as 


eJ  = A'  + B'  [f  (a)]^"  Uj; 
or 

eJ  = A'  + B'  u'' 
b a * 


(D-2) 


where  the  coefficient  B^  will  depend  on  the  angle  of  yaw. 
Both  A'  and  n are  assumed  to  remain  constant  for  all  a's. 
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At  each  angle  of  yaw,  the  hot-wire  velocity  calibration 
constants  were  determined  by  fitting  the  data  from  Set  1 to 
Equation  (D-2)  while  minimizing  the  deviation  of  the  data 
from  the  equation.  The  results,  shown  in  Figure  93,  indi- 
cate that  the  coefficient  B'  decreases  with  increasing 
angle  of  yaw  as  was  expected  from  Equation  (D-2) . On  the 
other  hand  the  coefficients  A'  and  exponent  n increase 
with  increasing  angle  of  yaw.  The  effect  of  forcing  A'  and 
n to  remain  constant  on  the  standard  deviation  of  fitting 
the  data  has  not  as  yet  been  examined.  However,  based  on 
the  available  results  it  is  likely  that  the  errors  intro- 
duced by  this  scheme  would  be  larger  than  those  achieved 
by  the  present  approach  and  discussed  in  the  following  part 
and  in  earlier  chapters. 

A second  set  of  calibration  data  was  requested  from 
MSU  and  Dr.  Foss  obtained  the  data  using  the  same  calibra- 
tion procedures  employed  in  the  present  work.  First,  he 
performed  a velocity  calibration  of  each  sensor  when  it  was 
normal  to  the  flow.  Next,  at  each  free  stream  velocity,  he 
varied  the  yaw  angle  of  the  probe  to  obtain  readings  of 
effective  velocities  through  the  use  of  Equation  (IV-1) . 

The  velocity  range  of  Set  1 was  between  10  ft/sec.  and 
110  ft/sec.  while  the  second  set  varied  from  15  ft/sec.  to 
55  ft/sec.  For  Set  1,  once  a free  stream  velocity  was 
chosen,  the  corresponding  output  voltages  were  calculated 
using  Eqtiation  (D-2)  at  each  yaw  angle.  Using  the  coeffi- 
cients retained  when  the  probe  was  nearly  normal  to  the 
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flow,  the  same  equation  was  utilized  to  calculate  the  effec- 
tive velocities.  Both  sets  of  data  were  then  analyzed  in 
the  manner  described  in  Chapter  IV  and  some  of  the  resulting 
correlations  are  shown  in  Figures  94  through  96. 

Sample  contours  of  constant  standard  deviation  in  the 
b-m  plane  for  one  velocity  of  the  yaw  calibration  data 
taken  from  Set  2 are  shown  in  Figure  94 . A comparison 
with  Figures  21  and  22  reveals  the  similar  trends  with  near- 
ly the  same  values  for  the  standard  deviation  contours.  The 
variations  with  Reynolds  number  of  the  ratio  of  improvement 
in  the  fit  of  the  data,  a|/a*,  and  the  yaw  coefficients  from 

Friehe  and  Schwarz  and  the  present  relation  for  both  sets  of 

* 

data  are  shown  in  Figures  95  and  96.  The  coefficients  b^ 
and  b*  increase  with  Reynolds  number  while  m*  and  the  ratio 
of  improvement  decrease  with  Reynolds  number.  Considering 
the  t/d  of  the  sensor  (£/d=200) , these  trends  can  be  com- 
pared to  those  presented  in  Figures  45  through  48.  Such 
comparison  indicates  that  in  each  case,  the  trends  are  iden- 
tical although  there  is  a fair  amount  of  scatter  in  the  mag- 
nitudes of  the  coefficients. 

These  results  indicate  that  the  behavior  of  the  wires 
used  at  I.I.T.  are  typical.  The  MSU  data  also  reaffirms 
the  fact  that  the  yaw  coefficients  depend  strongly  on  the 
probe  design  and  construction. 


Figure  1.  Effect  on  Present  Yaw  Relation  of  Changing 
Coefficient  b While  Holding  Exponent  m Constant 


f (a) 


Figure  2.  Effect  on  Present  Yaw  Relation  of  Changing 
Exponent  m While  Holding  Coefficient  b Constant 


Figure  7.  Turbulence  Correction  Factor,  g/h,  as  a Function  of  b and  m for 
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Instrumentation  for  Mean  Yaw  Calibration 
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Figure  11.  Photographs  of  Calibration  Tunnel,  Modified  to 
Produce  Streamwise  Periodic  Fluctuations,  with  Single 
Sensor  Probe  Mounted  in  Angularity  Chuck,  and  Instru- 
mentation for  Unsteady  Yaw  Calibration 


^ irr:  IJ.  Photograph*  of  Single-Sensor  Probe,  and 
% •«i4  Top  Vlaws  of  X-Probe;  All  Wires  have  ^/d>580 


Figure  13.  Instrumentation  Schematic  for  Mean  Velocity  Yaw  Calibration 
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Figure  18.  Typical  Variation  with  Offset  Angle  of  Stan- 
dard Deviation,  from  Present  Relation,  of  Mean  Yaw  Cali- 
bration Data  for  Properly  Shifted  Data 
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Figure  19.  Typical  Variation  of  Standard  Deviation  of  Mean 
Yaw  Calibration  Data  with  Probe  Misalignment  Using 
Three  Different  Yaw  Relations 


gure  21.  Contours  of  Constant  Standard  Deviation  in  b-m  Plane  for  a Single  Wire 
Probe  at  a Fixed  Reynolds  Number;  f/d=250  and  Re. =1.8 
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Figure  24.  Variation  with  Free  Stream  Velocity  of  Region 
of  Minimum  Standard  Deviation  for  X-Wire  Data  from  MSU 
and  a Standard-length  Single  Wire  Probe 
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25.  Effect  of  Probe  Misalignment  and  Reynolds  Number  on  Standard  Deviation 
Cosine  Law,  of  Mean  Yaw  Calibration  Data  for  a Single  Wire  Probe;  f/d=250 
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Figure  27.  Effect  of  Probe  Misalignment  and  Reynolds  Number  on  Standard  Deviation 
from  Present  Relation,  of  Ne2m  Yaw  Calibration  Data  for  a Single  Hire  Probe; 
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Figure  36.  Effect  of  Reynolds  Number  on  Ratio  of  Standard 
Deviation  Improvement  Between  Cosine,  Friehe  and 
Schwarz  and  Present  Relations  for  f/d*75  and  1000 
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Figure  37.  Variation  with  Reynolds  Number  of  Best  Fit  Co 
efficients  Determined  from  Friehe  and  Schwarz  and  Pres- 
ent Relations  for  f/d=1000 
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Figure  38.  Variation  of  Hot-Wire  Velocity  Calibration 
Constants  with  t/d 
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Figure  39.  Variation  of  Ratio  of  Standard  Deviation  Im- 
provement Between  Cosine  and  Present  Relations  with 
Sensor  f/d  at  a Fixed  Reynolds  Number;  Re^-1.8 


Figure  40.  Variation  of  Ratio  of  Standard  Deviation 
Improvement  Between  Frlehe  and  Schwarz  and  Present 
Relations  with  Sensor  t/d  at  a Fixed  Reynolds  Number 
Re,*! .8 


Figure  41.  Variation  of  Best  Fit  Coefficient  Determined 
from  Friehe  and  Schwarz  Relation  with  Sensor  //d  at  a 
Fixed  Reynolds  Number;  Re. *1.8 
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Figure  44.  Variation  of  Turbulence  Correction  Factors,  as 
Determined  from  Friehe  and  Schwarz  and  Present  Relations 
with  Sensor  t/d  at  a Fixed  Reynolds  Number;  Re.^1.8 


ire  45.  Combined  Effect  of  Reynolds  Number  and 
iviation  Improvement  Between  Friehe  and  Schwarz 
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Figure  53,  Variation  with  Overheat  Ratio  of  Turbulence 
Correction  Factor  and  Ratio  of  Standard  Deviation 
Improvement  Between  Friehe  and  Schwarz  and  Present 
Relations  for  a Standard-length  Probe;  f/d»580 
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Figure  54.  Variation  of  Best  Pit  Coefficients  and  Exponent 
Determined  from  Friehe  and  Schwarz  and  Present  Relations 
with  Overheat  Ratios  for  a Standard-length  Probe;  f/d«580 


1.025 


0.975 


• Re^j-S.b 
A Rej*  1.8 
A Re^  *0.9 

I/d*  145 


-A  a^a  ^ Aa 


Aa- DEGREES 


Figure  55.  Variation  of  Turbulence  Correction  Factor  and 
Ratio  of  Standard  Deviation  Improvement  Between  Friehe 
and  Schwarz  and  Present  Relations  with  Angular  Range  of 
Calibration;  .C/d«145 
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Figure  56.  Variation  of  Best  Fit  Coefficients  and  Expo 
nent  Determined  from  Friehe  and  Schwarz  and  Present 
Relations  with  Angular  Range  of  Calibration;  t/d»145 
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Figure  57.  Variation  of  Turbulence  Correction  Factor  and 
Ratio  of  Standard  Deviation  Improvement,  Between  Friehe 
and  Schwarz  and  Present  Relations,  with  Angular  Range  of 
Calibration  for  a Standard-length  Single  Wire  Probe;  i/d«580 
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Figure  58.  Variation  of  Best  Fit  Coefficients  and  Expo- 
nent Determined  from  Friehe  and  Schwarz  and  Present 
Relations  with  Angular  Range  of  Calibration  for  a Stan 
dard-length  Single  Wire  Probe;  i/d«580 
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Effective  Velocity  for  a Single  Wire  Probe 


SINGLE  WIRE  PROBE  ' 

OPEN  SYMBOLS-MEAN  CALIBRATION 
CLOSED  SYMBOLS -rms  CALIBRATION 
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Figure  72.  Variation  of  Best  Fit  Coefficient  Determined 
Using  Present  Relation  with  Sensor  t/A  for  Single  Wires 
X-Wires  and  Hot-Film  Probes 
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Figure  73.  Variation  of  Best  Fit  Exponent  Determined 
Using  Present  Relation  with  Sensor  for  Single  Wires, 
X-Wires  and  Hot-Film  Probes 
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Figure  74,  Variation  of  Turbulence  Correction  Factor 
Determined  Using  Present  Relation  with  Sensor  t/d  for 
Single  Wires,  X-Wires  and  Hot-Film  Probes 
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Figure  75.  Correlation  Between  Best  Fit  Coefficients  and 
Exponents  Determined  Using  Present  Relation  for  Single 
Wires,  X-Wires  and  Hot-Film  Probes 
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Figure  79.  Typical  Variation  with  Mean  Flow  Angle  of  Com 
puted  Differences  in  v'  and  uv  Obtained  Using  Different 
Yaw  Relations  for  a Sensor  with  Large  Deviations  from 
Cosine  Law 
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Figure  80.  Typical  Variation  with  Mean  Flow  Angle  of  Com- 
puted Differences  in  U and  u'  Intensity  Obtained  Using 
Different  Yaw  Relations  for  a Sensor  with  Large  Devia- 
tion from  Cosine  Law 
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Figure  81.  Typical  Variation  with  Mean  Flow  Angle  of  Com- 
puted Turbulence  Velocities  Obtained  Using  Different  Yaw 
Relations  for  a Standard-length  Single  Wire  Probe 
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Variation  with  Air  Temperature  of  Overheat  Ratio  and  Resistance  Differ 
Constant  Output  Voltage  Operation  at  Several  Fixed  Velocities 


Figure  86.  Temperature  Dependence  of  Hot-Wire  Velocity 
Calibration  Coefficients  for  Different  Overheat  Ratios 
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Figure  92.  Dependence  of  Region  of  Minimum  Standard  Devi 
ation  in  b-m  Plane  on  Reynolds  Number  for  Hot-Films 
in  Water 
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Figure  95.  Variation  with  Reynolds  Number  of  Best  Fit 
Coefficient  Determined  from  Friehe  and  Schwarz  Relation 
and  Ratio  of  Standard  Deviation  Improvement  Between 
Friehe  and  Schwarz  and  Present  Relations  for  MSU  Data 
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Figure  96.  Variation  with  Reynolds  Number  of  Best  Fit 
Coefficient  and  Exponent  Determined  From  Present  Rela 
tion  for  MSU  Data 
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